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We consider novel unusual effects in superconductor-ferromagnet (S/F) structures, fn particular 
we analyze the triplet component (TC) of the condensate generated in those systems. This compo- 
nent is odd in frequency and even in the momentum, which makes it insensitive to non-magnetic 
impurities. If the exchange field is not homogeneous in the system the triplet component is not 
destroyed even by a strong exchange field and can penetrate the ferromagnet over long distances. 
Some other effects considered here and caused by the proximity effect are: enhancement of the 
Josephson current due to the presence of the ferromagnet, induction of a magnetic moment in 
superconductors resulting in a screening of the magnetic moment, formation of periodic magnetic 
structures due to the influence of the superconductor, etc. We compare the theoretical predictions 
with existing experiments. 
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I. INTRODUCTION 

Although superconductivity has been discovered by H. Kammcrlingh Onnes almost one century ago (1911), the 
interest in studying this phenomenon is far from declining. The great attention to super conductivity within the las t 
15 years is partly due to the discovery of the high temperature superconductors (HTSC) IjBednorz and Mullen . Il986). 
which promises important technological applications. It is clear that issues such as the origin of the high critical 
temperature superconductivity, effects of external fields and impurities on HTCS, etc, will remain fields of interest 
for years to come. 

Due to the successful investigations of the HTSC and its possible technological applications, the interest in studying 
properties of traditional (low T c ) superconductors was not as broad. Nevertheless this field has also undergone 
a tremendous development. Technologically, the traditional superconductors are often easier to manipulate than 
high T c cuprates. One of the main achievements of the last decade is the making of high quality contacts between 
superconductors and normal metals (S/N), superconductors and ferromagnets (S/F), superconductors and insulators 
(S/I), etc. All these heterostructures can be very small with the characteristic sizes of submicrometers. 

This has opened a new field of research. The small size of these structures provides the coherence of superconducting 
correlations over the full length of the N region. The length of the condensate penetration into the N region £jy is 
restricted by decoherence processes (inelastic or spin-flip scattering). At low temperatures the characteristic length 
over which these decoherence processes occur may be quite long (a few microns) . Superconducting coherent effects in 
S/N nanostructures, such as conductance oscillations in an external magnetic field, were studied inten sively during 
the last decade (see for example the review articles bv lBeenakkerl l)l997j) : lLambert and Raimondil (jX908^) ^ . 

The interplay between a superconductor (S) and a normal metal (N) in simpler types of S/N structures (for 
example, S/N bilayers) has been unde r study for a long time and the main physics of th i s so c alled proximity effect is 
well described in the review articles bv lde Gennesl l)l964|) and lDeutscher and de Penned lTl969). In these works it was 
noticed that not only the superconductor changes the properties of the normal metal but also the normal metal has 
a strong effect on the superconductor. It was shown that near the S/N interface the superconductivity is suppressed 
over the correlation length £g, which means that the order parameter A is reduced at the interface in comparison with 
its bulk value far away from the interface. At the same time, the superconducting condensate penetrates the normal 
metal over the length £jv, which at low temperatures may be much larger than £<j. Due to the penetration of the 
condensate into the normal metal over large distances the Josephson effect is possible in S/N/S junctions with the 
thicknesses of the N regions of the order of a few hundreds nanometers. The Josephson effects i n S/N/S junctions 
were s t udied in many p apers and a good overview , both experimental and theoretical, is given bv lKulik and Yanson.1 
lll97ffl . lLikha"revl lll979l) . andU arone and Paternol l)l982() . 

The situation described above is quite different if an insulating layer / is placed between two superconductors. 
The thickness of the insulator in S/I/S structures cannot be as large as of the normal metals because electron wave 
functions decay in the insulator on atomic distances. As a consequence, the Josephson current is extremely small in 
S/I/S structures with a thick insulating layer. 

But what about S/F/S heterojunctions, where F denotes a ferromagnetic metal? In principle, the electron wave 
function can extend in the ferromagnet over a rather large distance without a considerable decay. However, it is 
well known that electrons with different spins belong to different energy bands. The energy shift of the two bands 
can be considered as an effective exchange field acting on the spin of the electrons. The condensate of conventional 
superconductors is strongly influenced by this exchange field of the ferromagnets and usually this reduces drastically 
the superconducting correlations. 

The suppression of the superconducting correlations is a consequence of the Pauli principle. In most superconductors 
the wave function of the Cooper pairs is singlet so that the electrons of a pair have opposite spins. In other words, 
both the electrons cannot be in the same state, which would happen if they had the same spin. If the exchange field of 
the ferromagnet is sufficiently strong, it tries to align the spins of the electrons of a Cooper pair parallel to each other, 
thus destroying the superconductivity. Regarding the S/F interfaces and the penetration of the condensate into the 
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ferromagnet, these effects mean that the superconducting condensate decays fast in the region of the ferromagnet. 
A rough estimate leads to the conclusion that the ratio of the condensate penetration depth in ferromagnets to the 
one in non-magnetic metals with a high impurity concentration is of the order of y^T c /h, where h is the exchange 
energy and T c is the critical temperature of the superconducting transition. The exchange energy in conventional 
ferromagnets like Fe or Co is several orders of magnitude higher than T c and therefore the penetration depth in the 
ferromagnets is much smaller than that in the normal metals. 

Study of the proximity effect in the S/ F str uctures s tarted n ot long ago but it h a s already evolv e d int o 
a very active field of rese a rch ( for a review see iBuzdinl l)2005a|) : iGolubov etall i|2004D : llzvumov et~al\ (|2002ft : 
iLvuksvutov and Pokrovskvl <)2004|) '). The effect of the suppression of superconductivity by the ferromagnetism is 
clearly seen experimentally and this corresponds to the simple picture of the destruction of the singlet superconduc- 
tivity by the exchange field as discussed above. 

At first glance, it seems that due to the strong suppression of the superconductivity the proximity effect in S/F 
structures is less interesting than in the S/N systems. However, this is not so because the physics of the proximity 
effect in the S/F structures is not exhausted by the suppression of the superconductivity and new very interesting 
effects come into play. Moreover, under some circumstances superconductivity is not necess arily suppressed by the 
ferromagnets because th e presence of the latter may lead to a triplet superconducting pairing l)Bergeret et all . l2001al 
iKadierobov et al\ . 12001). In some cases not only the ferro magnetism tends to destroy the superconduct i vity b ut also 
the superconductivity may suppress the ferromagnetism l(Bergeret et al .1 1200(1 iBuzdin and Bula,evskuL M 988ft . This 
may concern "real" strong ferromagnets like iron or nickel with a Curie temperature much larger than the transition 
temperature of the superconductor. 

In all, it is becoming more and more evident from recent experimental and theoretical studies that the variety 
of non-trivial effects in the S/F structures exceeds considerably what one would have expected before. Taking into 
account possible technological applications, there is no wonder that S/F systems attract nowadays a lot of attention. 

This review article is devoted to the study of new "exotic" phenomena in the S/F heterojunctions. By the word 
"exotic" we mean phenomena that could not be expected from the simple picture of a superconductor in contact 
with a homogeneous ferromagnet. Indeed, the most interesting effects should occur when the exchange field is not 
homogeneous. These non-homogeneities can be either intrinsic for the ferromagnetic material, like e.g. domain walls, 
or arise as a result of experimental manipulations, such as multilayered structures with different directions of the 
magnetization, which can also be spoken of as a non-homogeneous alignment of the magnetic moments. 

Of course, we are far from saying that there is nothing interesting to be seen when the exchange field is homogeneous. 
Although it is true that in this case the penetration depth of the superconducting condensate into the ferromagnet 
is short, the exponential decay of the condensate function into ferromagnets is accompanied by oscillations in space. 
These oscillations lead, for example, to oscillations of the critical superconducting temper ature T r and the critical 
Josephson current I r in S/F stru ctures as a function of the thickness (If- Being predicted bv lBuzd in and^ Kupriya noyl 
(1990) and lRadovic et q?J l|l99l|) . the observation of such oscillatory behavior was first reported bv lJiang"e7anifl995[) 
on Gd/Nb structures . Indications to a non-monotonic behav i or of T r as a funct i on of dp was also reported by 
iMercaldo etall l|l996fl : lMuhge et all l)l996j) : l0biand et all l[l999fl : IStrunk et all l|l994j) : IVelez et aD l)l999fl : IWong et all 
1119861) . 

However, in other expe riments the dependence of T c on dp was monotonic. For example in Ref. 
ijBourgeois and Dvnesl 12002^) the critical temperat ure of the bilayer Pb /Ni decreased by increasing the F layer thick- 
ness d p in a monotonic way. In the experiments by Miihge et al. (1998) on Fe/Nb/Fe structures and bv lAarts et all 
l)l997l) on V/Fe systems both a monotonic and non-monotonic behavior of T c has been observed. This different 
behavior was attributed to changes of the transmittance of the S/F interface. A compr ehensive analysis taking into 
account the samples quality was made for different materials bv lChien and Reichl |)1999|) . 

More convincing results were found by measuring the Josephson critical current in a S/ F/S junction. Due to 
the oscillatory behavior of the superconducting condensate in the F region the critical J osephson current should 
change its sign in a S/F/S junction {it— junction). This phenomenon predicted long ago bv iBulaevskii et al\ 
has been confirmed experimentally only recently ijBauer et all . 120041 iBlum et all 120021: iKontos et all 12002 
iRvazanov et aLLl200lHSellier et alll2004|) . 

Experiments o n tra nspor t properties of S/F structures were also performed in the last years. For example, 
iPetrashov et all 1^999) and iGiroud et all l)l998j) observed an unexpected decrease of the resistance of a ferromag- 
netic wire attached to a superconductor when the temperature is lowered below T c . In both of the experiments strong 
ferromagnets Ni and Co, respectively, were used. One would expect that the change of the resistance must be very 
small due to the destruction of the superconductivity by the ferromagnets. However, the observed drop was about 
10% and this can only be explained by a long-range proximity effect. 

This raises a natural question: how can such long range superconducting effects occur in a ferromagnet with a 
strong exchange field? We will see in the subsequent chapters that provided the exchange field is not homogenous a 
long-range component of the condensate may be induced in the ferromagnet. This component is in a triplet state and 
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can penetrate the F region over distances comparable with £jv, as in the case of a normal metal. 
We outline now the structure of the present review. 

In Chapter [n] we discuss the proximity effects in S/N structures and S/F structures with a homogeneous magneti- 
zation. Chapter[n]may serve as an introduction into the field. The main results illustrated there have been presented 
in other reviews and we discuss them here in order to give the reader an impression about works done previously. 
Chapter [H] can also help in getting the basic knowledge about calculational methods used in subsequent chapters. 
One can see from this discussion that already homogeneous ferromagnets in contact to superconductors lead to new 
and interesting physics. 

Nevertheless, the non-homogeneities bring even more. We review below several different effects arising in the non- 
homogeneous situation. It turns out that a non-homogeneous alignment of the exchange field leads to a complicated 
spin structure of the superconducting condensate. As a result, not only the singlet component of the condensate exists 
but also a triplet one with all possible projections of the total spin of the Cooper pair (S z — 0, ±1). In contrast to the 
singlet component, the spins of the electrons in the triplet one with S z = ±1 are parallel to each other. The condensate 
(Gor'kov) function f tr of the triplet state is an odd function of the Matsubara frequency 1 . The singlet part f sng is, as 
usual, an even function of a; but it changes sign when interchanging the spin indices. This is why the anticommutation 
relations for the equal-time functions f tr {t,t) and f sn g(t,t) remain valid; in particular, f tr (t,t) = and f S ng(t,t) =/= 0. 
Therefore the superconductivity in the S/F structures can be very unusual: alongside with the usual BCS singlet 
part it may contain also the triplet part which is symmetric in the momentum space (in the diffusive case) and odd 
in frequency. Both components are insensitive to the scattering by non-magnetic impurities and hence survive in the 
S/F structures even if the mean free path I is short. When generated, the triplet component is not destroyed by the 
exchange field and can penetrate the ferromagnet over long distances of the order of £n = \J Dp /litT . 

In Chapter IIIII we analyze properties of this new type of superconductivity that may arise in S/F structures. We 
emphasize that this triplet superconductivity is generated by the exchange field and, in the absence of the field, one 
would have the conventional singlet pairing. 

The superconductor-ferromagnet multilayers are a very interesting and natural object for observation of Josephson 
effects. The thickness of both the superconductor and ferromagnetic layers, as well as the transparency of the interface, 
can be varied experimentally. This makes possible a detailed study of many interesting physical quantities. As we have 
mentioned, an interesting manifestation of the role played by the ferromagnetism is the possibility of a 7r-junction. 

However, this is not the only interesting effect and several new ones have been recently proposed theoretically. 
As not so much time has been passed, they have not been confirmed experimentally unambiguously but there is no 
doubt that proper experiments will have been performed soon. In Chapter II VI we discuss new Josephson effects in 
multilayered S/F structures taking into account a possible change of the mutual direction of the magnetization in 
the ferromagnetic layers. We discuss a simple situation when the directions of the magnetic moments in a SF/I/FS 
structure are collinear and the Josephson current flows through an insulator (J) but not through the ferromagnets. 
Naively, one could expect that the presence of the ferromagnets leads to a reduction of the value of the critical 
current. However, the situation is more interesting. The critical current is larger when the magnetic moments of 
the F-layers are antiparallel than when they are parallel. Moreover, it turns out that the critical current for the 
antiparallel configuration is even larger than the one in the absence of any ferromagnetic layer. In other words, the 
ferromagnetism can enhance the critical current ijBergeret et all l200Ib|) 

Another setup is suggested in order to observe the odd triplet superconductivity discussed in Chapter IIIII Here 
the current should flow through the ferromagnetic layers. Usually, one could think that the critical current would 
just decay very fast with increasing the thickness of the ferromagnetic layer. However, another effect is possible. 
Changing the mutual direction of the additional ferromagnetic layers one can generate the odd triplet component of 
the superconducting condensate. This component can penetrate the ferromagnetic layer as if it were a normal metal, 
leading to large values of the critical current. 

Such structures can be of use for detecting and manipulating the triplet component of the condensate in experiments. 
In particular, we will see that in some S/F structures the type of superconductivity is different in different directions: 
in the longitudinal direction (in-plane superconductivity) it is caused mainly by the singlet component, whereas in the 
transversal direction the triplet component mainly contributes to the superconductivity. We discuss also possibilities 
of an experimental observation of the triplet component. 

Although the most pronounced effect of the interaction between the superconductivity and ferromagnetism is the 
suppression of the former by the latter, the opposite is also possible and this is discussed in Chapter Ivl Of course, a 
weak ferromagnetism shoul d be strongly affected b y the superconductivity and this situation is realized in so called 
magnetic superconductors ijBulaevskii et all Il985|) . Less trivial is that the conventional strong ferromagnets in the 
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S/F systems may also be considerably affected by the superconductivity. This can happen provided the thickness 
of the ferromagnetic layer is small enough. Then, it can be energetically more profitable to enforce the magnetic 
moment to rotate in space than to destroy the superconductivity. If the period of such oscillations is smaller than 
the size of the Cooper pairs £s, the influence of the magnetism on the superconductor becomes very small and the 
superconductivity is preserved. In thick layers such an oscillating structure (cryptoferromagnetic state) would cost 
much energy and the de struction of the superconductivity is more f avorable. Results of several experiments have been 
interpreted in this way l)Garifullin et al ., 200i lMuhge et a/11199^ . 

Another unexpected phenomenon, namely, the inverse proximity effect is also presented in Chapter It turns out 
that not only the superconducting condensate can penetrate the ferromagnets but also a magnetic moment can be 
induced in a superconductor that is in contact with a ferromagnet. This effect has a very simple explanation. There 
is a probability that some of the electrons of Cooper pairs enter the ferromagnet and its spin tends to be parallel to 
the magnetic moment. At the same time, the spin of the second electron of the Cooper pair should be opposite to 
the first one (the singlet pairing or the triplet one with S z — is assumed). As a result, a magnetic moment with the 
direction opposite to the magnetic moment in the ferromagnet is induced in the superconductor over distances of the 
superconducting coherence length £s. 

In principle, the total magnetic moment can be completely screened by the superconductor. Formally, the appear- 
ance of the magnetic moment in the superconductor is due the triplet component of the condensate that is induced 
in the ferromagnet F and penetrates into the superconductor S. It is important to notice that this effect should 
disappear if the superconductivity is destroyed by, e.g. heating, and this gives a possibility of an observation of 
the effect. In addition to the Meissner effect, this is one more mechanism of the screening of the magnetic field by 
superconductivity. In contrast to the Meissner effect where the screening is due to the orbital electron motion, this is 
a kind of spin screening. 

Finally, in Chapter VI we discuss the results presented in the review and try to anticipate future directions of 
the research. The Appendix [S] contains necessary information about the quasiclassical approach in the theory of 
superconductivity. 

We should mention that s everal review articles on S/F related topics have been published recently ijBuzdinl 
l2005al iGolubov et all . 12004 llzvumov et aH 120021: iLvuksvutov and Pokrovskvl |2004|) . In these reviews various 
properties of the S/F structures are discussed for the case of a homogeneous magnetization. In the review by 
ILvuksvutov and Pokrovskvl l)2004|) the main attention is paid to effects caused by a magnetic interaction between 
the ferromagnet and superconductor (for example, a spontaneous creation of vortices in the superconductor due to 
the magnetic interaction between the magnetic moment of vortices and the magnetization in the ferromagnet). We 
emphasize that, in contrast to these reviews, we focus on the discussion of the triplet component with all possible 
projections of the magnetic moment (S z = 0,±1) arising only in the case of a nonhomogeneous magnetization. In 
addition, we discuss the inverse proximity effect, that is, the influence of superconductivity on the magnetization M 
of S/F structures and some other effects. Since the experimental study of the proximity effects in the S/F structures 
still remains in its infancy, we hope that this review will help in understanding the conditions under which one can 
observe the new type of superconductivity and other interesting effects and hereby will stimulate experimental activity 
in this hot area. 



II. THE PROXIMITY EFFECT 

In this section we will review the basic features of the proximity effect in different heterostructures. The first part 
is devoted to superconductors-normal metals structures, while in the second part superconductors in contact with 
homogeneous ferromagnets are considered. 

A. Superconductor-normal metal structures 

If a superconductor is brought in contact with a non-superconducting material the physical properties of both ma- 
terials may change. This phenomenon called the proximity effect has been studied for many years. Both experiments 
and theory show that the properties of superconducting layers in contact with insulating (I) materials remain almost 
unchanged. For example, for superconducting films evaporated on glass substrates, the critical temperature T c is very 
close to the bulk value. However, physical properties of both metals of a normal metal/superconductor (N/S, see 
FIG.^) heterojunction with a high N/S interface conductance can change drastically. 

Study of the prox i mity effect goes back to the beginning of the 1960's and was reviewed in many publications 
(see, e.g. Ide Gennesl ()l964fl and iDeutscher and de Gennesl l)l969(l '). It was found that the critical temperature of the 
superconductor in a S/N system decreased with increasing TV layer thickness. This behavior can be interpreted as 
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FIG. 1 S/N bilayer. 



the breaking down of some Cooper pairs due to the penetration of one of the electrons of the pairs into the normal 
metal where they are no longer attracted by the other electrons of the pairs. 

At the same time, penetrating into the normal metal the Cooper pairs induce superconducting correlations. For 
example, the influence of the superconductivity on the physical properties of the TV metal manifests itself in the sup- 
pression of the density of states. Experiments dete rmining the density of states of S/N bilayers with the help of tunnel- 
ing spectroscopy were performed many years ago ijAdkins a ndJKm gtonLll969tlTot)licar and Finnemore L Il977l). While 
spatia lly resolved density of states were later measured bv lAnthore^ra^J l)2003() : TGue ? roriT^ni^9q ): iGupta et all 
l)2004j) (see FIG.EJ). The simplest way to describe the proximity effect is to use the Ginzburg-Landau (G-L) equation 




V(mV) 

FIG. 2 Tunnelling density o f states measured at 60 mK at the Au surface of dif f erent Nb/Au bilayers samples with varying Au 
thickness djv- Adapted from lGupta. Crtinon. Moussv. Pannetier. and Courtoisl (120041) . 



for the order parameter A ijOinzburg and LandaiJ . ll95flL This equation is valid if the temperature is close to the 
critical temperature of the superconducting transition T c . In this case all quantities can be expanded in the small 
parameter A/T c and slow variations of the order parameter A in space. 
Using the G-L equation written as 



d 2 A (r) 
<9r 2 



A (r) (sgn(T c , w ,5 — T) — A 2 (r) / A 2 ,) = 



(2.1) 



one can describe the spatial distribution of the order parameter in any N/S structure. Here £,gl is the coherence 
length in the N and S regions at temperatures close to the critical temperatures T c n,s- in the diffusive limit this 
length is equal to 



£gl - JnD NjS /8\T - T c 



cN,S\ 



(2.2) 



where Dn^s is the diffusion coefficient in the N and S regions. The quantity Ao is the bulk value of the order 
parameter in the superconductor S. It vanishes when T reaches the transition temperature T c . 

It should be noticed though, that the region of the applicability of Eq. i|2.1|> for the description of the S/N 
contacts is rather restricted. Of course, the temperature must be close to the transition temperature T c but this is 
not sufficient. The G-L equation describes variations of the order parameters correctly only if they are slow on the 
scales vf/Tc for the clean case or WDn\s/T c in the diffusive "dirty" case. This can be achieved if the normal metal 
is a superconducting material taken at a temperature exceeding its transition temperature T c n and the transition 
temperatures T c g and T c n are close to each other. If this condition is not satisfied (e.g. T c n = 0) one should use 
more complicated equations even at temperatures close to T c s, as we show below. 



7 



It follows from Eq. I|2.1(l that in the S region, far from the N/S interface, the order parameter A (r) equals the 
bulk value A , whereas in the N region A (r) decays exponentially to zero on the length £jy 
The order parameter A (r) is related to the condensate function (or Gor'kov function) 

/(M')=<^T(*W)> (2-3) 

via the self-consistency equation 

A N , s (t) = X N ,sf{t,t) , (2.4) 

where Xn,s is the electron-electron coupling constant leading to the formation of the superconducting condensate. 

Ea. l|2.1|l describes actually a contact between two superconductors with different critical temperatures T c n^s, when 
the temperature is chosen between T c $ and T c n. In the case of a real normal metal the coupling constant A at is equal 
to zero and therefore An = 0. However, this does not imply that the normal metal does not possess superconducting 
properties in this case. The point is that many important physical quantities are related not to the order parameter 
A but to the condensate function /, Eq. i|2.3[l . For example, the non-dissipative condensate current j$ is expressed in 
terms of the function / but not of A. If the contact between the N and S regions is good, the condensate penetrates 
the normal metal leading to a finite value of jg ^ in this region. 

In the general case of an arbitrary Ajv it is convenient to desc ribe the penet ration of the condensate (Cooper pairs) 
into the N region in the diffusive limit by the Usadel equation ijUsadell Il970|) which is valid for all temperatures and 
for distances exceeding the mean free path I. This equation determines the so called quasiclassical Green's functions 
(see Appendix 0) which can be conveniently used in problems involving length scales larger than the Fermi wave 
length Xf and energies much smaller than the Fermi energy. Alternatively, one could try to find an exact solution 
(the normal and anomalous electron Green's functions) for the Gor'kov equations, but this is in most of the cases a 
difficult task. 

In order to illustrate the convenience of using the quasiclassical method we calculate now the change of the tunnelling 
density of states (DOS) in the normal metal due to the proximity effect with the help of the Usadel equation. The 
DOS is a very important quantity that can be measured experimentally and, at the same time, can be calculated 
without difficulties. 

We consider the S/N structure shown in FIG.H an d assume that the system is diffusive (i.e. the condition er << 1 
is assumed to be fulfilled, where r is the momentum relaxation time and e is the energy) and that the transparency 
of the S/N is low enough. In this case the condensate Green's function /(e) = J dtf(t — t') exp(ie(t — i')) is small in 
the N region and the Usadel equation can be linearized (see Appendix |A"|) . 

Assuming that the boundary between the superconductor and normal metal is flat and choosing the coordinate x 
perpendicular to the boundary we reduce the Usadel equation in the N region to the form 

D N d 2 f/dx 2 + 2ief = 0, (2.5) 

where Dn = vfI/3 is the classical diffusion coefficient. 

The solution of this equation can be found easily and we write it as 

/ = / cxp(-aV-2ie/ J D JV ) , (2.6) 

where /q is a constant that is to be found from the boundary conditions. 

We see that the solution for the condensate function / decays in the N region exponentially at distances inversely 
proportional to y/e. In many cases the main contribution to physical quantities comes from the energies e of the order 
of the temperature, e ~ T . This means that the superconducting condensate penetrates the N region over distances 
of the order of £jy = \J Dm /2nT. At low temperatures this distance becomes very large, and if the thickness of the 
normal metal layer is smaller than the inelastic relaxation length, the condensate spreads throughout the entire N 
region. 

In order to calculate the DOS it is necessary to know the normal Green's function g which is related to the condensate 
function / via the normalization condition (see Appendix lAf 

g 2 -f 2 = l (2.7) 

Eqs. Ij2.5[l and 1|2.7|1 are written for the retarded Green's function (/ = f R , see Appendix lA"|l . They are also valid 
for the advanced Green's functions provided (e + iO) is replaced by (e — iO). The normalized density-of-states (we 
normalize the DOS to the DOS of non-interacting electrons) v(e) is given by the expression 



v{e) = Reg(e) 



(2.8) 
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As the condensate function / is small, a correction 8v to the DOS due to the proximity effect is also small. In the 
main approximation the DOS v is very close to its value in the absence of the superconductor, vsjI. Corrections to 
the DOS 8v are determined by the condensate function /. From Eq. 1|2.7(1 one gets 

6u « f/2 . 

Now we consider another case when the function / is not small and the correction 5v is of the order of unity. Then 
the linearized Eq. Q2.5JI may no longer be used and we should write a more general one. For a S/N system the general 
equation can be written as (see Appendix 0) 

- iDs,Nd{gdg/dx)s,N/dx + e[f 3 ,g s , N ] + [A s ,g s , N } = . (2.9) 

This non-linear equation contains the quasiclassical matrix Green's function g. Both normal g and anomalous Green's 
functions / enter as elements of this matrix through the following relation (the phase in the superconductor is set to 
zero) 

gN = gnh + /jvif 2 , (2-10) 

where Tj, i = 1,2,3 are Pauli matrices and [A, B] = AB — BA is the commutator for any matrices A and B. 

We consider a flat S/N interface normal to the x-axis. The normal metal occupies the region < x < d/v- We 
assume that in the normal metal N there is no electron-electron interaction (Ajv = 0, see Ea. l|2.4|l ) so that in this 
region the superconducting order parameter vanishes, A^v = 0. In the superconductor the matrix As has the structure 
A s = Aif 2 . 

At large distances from the S/N interface the Green's functions gs of the superconductor do not depend on 
coordinates and the first term in Eq. H2.9J) can be neglected. Then we obtain a simpler equation 

e[T$,gs] + &[in,gs] = 0- (2.11) 

The solution for this equation satisfying the normalization condition (|2.7|) is 

9BCS = e/U f B cs = A/& , (2.12) 

where £ e = y/e 2 — A 2 . Eq. Ij2.12|l is just the BCS solution for a bulk superconductor. 

In order to find the matrix g(x) both in the S and N regions, Ea. H2.9fl should be complemented by boundary 
conditions and this is a non-trivial problem. Starting from the initial Hamiltonian H to t, Eq. (|2.22() . one does not 
need boundary conditions at the interface between the superconductor and the ferromagnet because the interface can 
be described by introducing a proper potential in the Hamiltonian. In this case the self-consistent Gor'kov equations 
can be derived. 

However, deriving the Usadel equation, Eq. i|A18|) . we have simplified the initial Gor'kov equations using the 
quasiclassical approximation. Possible spatial variation of the interface potential on a very small scale, due to the 
roughness of the interface cannot be included in the quasiclassical equations. Nevertheless, this problem is avoided 
deriving the quasiclassical equations at distances from the interface exceeding the wavelength. In the diffusive case 
one should go away from the interface to distances larger than the mean free path I. In order to match the solutions 
in the superconducting and non-superconducting regions one should solve exact the equations near the interface and 
compare the asymptotic behavior of this solution at large distances with the solutions of the Usadel equation. This 
procedure i s equiva l ent to solvi ng the quasiclassical equat i ons wi th some boundary conditions. These conditions were 
derived by IZaitsevI l|l984l) and iKunrianov and Lukichevl |l988h (see also Appendix EI where these conditions are 
discussed in more details). For the present case they can be written as 

2js,N(gdg/dx) s ,N = [<?s, <?iv]|x=o (2-13) 

where js,N = Rb&s,N, Rb, measured in units flcm 2 , is the S/N interface resistance per unit area in the normal state, 
and a$ n are the conductivities of the S and N metals in the normal state. 

We assume that the thickness of the normal metal cLm is smaller than the characteristic penetration length £at(c) = 
y5jv/e for a given energy e, that is 2 e << Dn /d 2 N = Et- Then the functions g and / remain almost constant over 



2 The quantity Ej<h = £*iv / d% is the so called Thouless energy 
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the thickness of the metal, and for finding them, one can average the Usadel equation over the thickness. In other 
words, we assume that the thickness dpf of the N layer satisfies the inequality 



d N « y/D N /e, e ~ e bN (2.14) 

(ebN is a characteristic energy in the DOS of the N layer) and average Eq. i|2.9[) over the thickness djv considering 
gpf as a constant in the second term of this equation. Using the boundary condition, Eq. l|2.13|) . the first term in Eq. 
I|2.9(l can be replaced after the integration by the commutator [gs, gN]\x=o- At x — (In the product (gdg / 8x)n is 
zero because t he barrier resistance R b {d,N) is infinite (the current cannot flow into the vacuum). Finally, we obtain 
<|Zaitsevl fl990^ 

(e + ie b g s {0))[f 3 ,g N } + e bN ifs(0)[if 2 ,g N ] = . (2.15) 

where e b N — (f jv /27jvrfjv) is a new characteristic energy that is determined by the S/N interface resistance R b . This 
equation looks similar to Ea. (|2.11|l after making the replacement gs — > c/n- The solution is similar to the solution 

firm 

9n = e/le! /jv = hN/L , (2-16) 

where e = e + iebjv3s(0), £ c = \Jt 2 — e' bN , e b N = ewv*/s(0). Therefore the Green's functions in the N layer gN and 
Jn are determined by the Green's functions on the S side of the S/N interface <?s(0) and /s(0). In order to find the 
values of gs(0) and /s(0), one has to solve Eq. (|2.9I) on the superconducting side (x < 0). However, provided the 
inequality 

In /is = <J N /a s « 1 (2-17) 

is fulfilled one can easily show that, in the main approximation, the solution in the S region coincides with the solution 
for bulk superconductors Q2.12p . If the transparency of the S/N interface is not high, e b N << A, the characteristic 
energies e ~ e b N are much smaller than A and the functions gs(0) and /s(Q) are equal to: gs(0) w ffscs(O) ~ e/»A, 
fs(0) w /bcs(0) ~ 1/i- For these energies the functions gN and Jn have the same form as the BCS functions gscs 
and Jbcs H2.12|) with the replacement A — > e b N 

V 6 _ e bN V 6 _ £ bAf 

where ebAr = Dn / (2R b <JNdN)- The energy ehAr can be represented in another form 

i* 2 ,Rq , k v F v F T b 

ebN = Y { R^ F )h d^ = ^ ( T } • (2 - 19) 

J?q = fi/e 2 is the resistance quantum, v F and kp are the Fermi velocity and wave vector. When obtaining the latter 
expression, we used a relation between the barrier resistance R b an d an effective coefficient of transmission T b through 
the S/N interface ijKuprianov and Lukichevl . Il988l IZaitsevl Il984h : R b a n = (2/3)(l/T b ), where I = vft is the mean 
free path, T b =< T{&) cos 8/ (1 — T(0) >, 9 is the angle between the momentum of an incoming electron and the vector 
normal to the S/N interface, and T{&) is the angle dependent transmission coefficient. The angle brackets mean an 
averaging over 6. 

An imp ortant result follo ws from Ea. (|2.18|l : the DOS is zero at |e| < e b N, ie., e b N is a minigap in the excitation 
spectrum (jMcMillanl . fl968). Remarkably, in the considered limit e b N « A the value of e b N does not depend on A, 
but is determined by the interface transparen cy or, in other w ords, by the interface resistance R b . The appearance of 
the minigap is related to Andreev reflections ijAndreevl Il964j) . 

Eq. (|2.19|) for the minigap is valid if the inequalities (|2.14|) and e b N < A are fulfilled. Both inequalities can be 
written as 

(D N /A)/d b <d N <d b (2.20) 

where d b = 2R b aN is a characteristic length. In the case of a small interface resistance R b or a large thickness of the 
N layer, that is, if the condition y/ Dn / A , d b < dN is fulfilled, the value of the minigap in the N layer is given by 
l|Golubov and Kuprivanovt ^996) 

e bN = cx^f (2.21) 
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where c\ is a factor of the order 1. This result has been obtained from a numerical solution o f the Usadel equation . 
The DOS for the case of arbitrary thickness d/v and interface transparency was calculated bv lPilgram et al\ (|2000). 

The situation changes in the clean limit. Let us consider, for example, a normal slab of a thickness in contact 
with an infini te supe r condu ctor. If the Thouless energy Exh = vp/dx i s less than A, th en discrete energy levels e n 
appear l)Saint-Jamesl [1964') in the N region due to Andreev reflect ions llAndreevl |l9 64). As a result, the DOS has 
sharp peaks at e = e n (for a recent review see iDeutscherl ((2005)). If Etk is much larger than A, the DOS v{e) 
is zero at e = and increases with increasing the energy e (no gap). However, this is true only for such a simple 
geometry. For samples of more complicated shape s the behavior of the DOS v(e) depends on whether the electron 
dynamics in the N region is chaotic or integrable ijB eenakkeil Il997t lLodder and Yu.V.Nazarovi Il998t iMelsen et a/1 
119961 iPikram et all l2000l iTaras-Semchuk and Altlandl 1200 ll) 

Finally, it has been shown bv lAltland^^IT|2000|) and lOstrovskv et al\ l)200l|) that mesoscopic fluctuations smear 
out the singularity in the DOS at |e| = ebN and the DOS in the diffusive limit is finite, although small, for |e| < e^. 
The minigap discuss e d above has been obser ved on a Nb/Si bilayer system and on a Pb/Ag granular system by 
iHeslinga etall (|l994|k iKouh and Vallesl l|2003h . respectively. 

From this analysis we see that the proximity effect changes the DOS of the normal metal which acquires supercon- 
ducting properties. In the next section we will focus our attention on the case that the normal metal is a ferromagnet. 
We will see that new interesting physics will arise from the mutual interaction of superconductivity and magnetism. 



B. Superconductor-ferromagnet structures with an uniform magnetization 

In this section we consider the proximity effect between a superconductor S and a ferromagnet F. We assume that 
the ferromagnet is a metal and has a conduction band. In addition, there is an exchange field due to spins of electrons 
of other bands. 

As has been already mentioned, the effective exchange field acts on spins of the conduction electrons in the fer- 
romagnet, and an additional term H ex describing this action appears in the total Hamiltonian (for more details see 
Appendix El 

Htot = H + H ex (2.22) 



H ex = - J d 3 r^+ (r) (h (r) a af} ) $ (r) dv (2.23) 

where ip + (ip) are creation and destruction operators, h is the exchange field, a a p are Pauli matrices, and a, f3 are 
spin indices. The Hamiltonian H stands for a non-magnetic part of the Hamiltonian. It includes the kinetic energy, 
impurities, external potentials, etc. and is sufficient to describe all properties of the system in the absence of the 
exchange field. 

The energy of the spin-up electrons differs from the energy of the spin-down electrons by the Zeeman energy 2h. Due 
to the presence of the term H ex describing the exchange interaction all functions, including the condensate Green's 
function /, are generally speaking non-trivial matrices in the spin space with non-zero diagonal and off-diagonal 
elements. 

The situation is simpler if the direction of the exchange field does not depend on coordinates. In this case, choosing 
the z-axis along the direction of h one can consider electrons with spin "up" and "down" separately. In this Section 
we concentrate on this case. This can help the reader to understand several interesting effects and get an intuition 
about what one can expect from the presence of the exchange field. The results of this section will also help in 
understanding which effects in the superconductor-ferromagnet structures can be considered as rather usual and what 
kind of behavior is "exotic" . We will see that the exotic phenomena occur in cases when the exchange field is not 
homogeneous and therefore postpone their discussion until the next chapters. 

If the exchange field h is homogeneous the matrix / describing the condensate / is diagonal and can be represented 
in the form 

/ = h&3 + (2.24) 

where f% is the amplitude of the singlet component and /o is the amplitude of the triplet component with zero 
projection of the magnetic moment of Cooper pairs on the z axis (S z = 0). Note that in the case of a S/N structure 
the condensate function has a singlet structure only, i.e. it is proportional to 03. The presence of the exchange field 
leads to the appearance of the triplet term proportional to ctq 
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The amplitudes of the singlet and tripl et components are related to the correlation functions (ipaipp) as follows 
llLeggetl Il975t IVollhardt and WolfleL ll99CD 

h{t) ~ (V» T (*)Vi(0))-(^(t)^T(0)) , 

fo(t) ~ (V» T (t)Vl(0)) + (^(*)^ T (0)) . (2-25) 

One can see that a permutation of spins does not change the function /3(C)), whereas such a permutation leads to a 
change of the sign of /q (0). This means that the amplitude of the triplet component taken at equal times is zero in 
agreement with the Pauli exclusion principle. Later we will see that in the case of a non-homogeneous magnetization 
all triplet components including (ip^(t)ip^(0)) and (V^ WVuW) differ from zero. 

Once one determines the condensate function, Eq. I|2.24|l . one is able to determine physical quantities ,as DOS, the 
critical temperature Tc, or the Joscphson critical current through a S/F/S junction. 

Next paragraphs are devoted to a discussion of these physical properties in F j S systems with homogeneous mag- 
netization. 



1. Density of states (DOS) 

In this section we discuss the difference between the DOS in S/N and S/F structures. General equations for the 
quasiclassical Green's functions describing the system can be written but they are rather complicated (see Appendix 
|A"|) . In order to simplify the problem and, at the same time, give the basic idea about the effects it is sufficient to 
consider some limiting cases. This will be done in the present section leaving the general equations for the Appendix 

m 

In the case of a weak proximity effect, the condensate function / is small outside the S region. We consider again 
the diffusive limit. Then, the general Eq. (|A18(1 can be linearized and one obtains an equation for the matrix / 
similar to Ea. (|2.5|) but containing an extra term due to the exchange field h 

D F d 2 f F /dx 2 + 2i(ea Q + ha 3 )f F = . (2.26) 

The subscript F stands for the F region. 

In the absence of the exchange field h, Eq. (|2.26l) reduces to Eq. I|2.5|l . It is important to emphasize that Eq. (|2.26() 
is valid for a homogeneous h only. Any variation of h in space makes the equation much more complicated. 

Eq. H2.26J) should be complemented by boundary conditions which take the form (see Appendix 0) 

lF df F /dx = -f s (2.27) 

where jp — Rb(J F , Rb is the boundary resistance per unit area, a F is the conductivity of the F region, /f,s are 
the condensate matrix functions in the F and S regions. Since we assume a weak proximity effect, a deviation of 
the fs from its BCS value fscs = frsfscs is small. Therefore on the right-hand side of Eq. l|2.27|) one can write 
fg Rj o"3/bcS7 where Jbcs is defined in Ea. l|2.9|l . At the ferromagnet /vacuum interface the boundary condition is 
given by the usual expression d x f F — 0, which follows from the condition Rf, — > 00. 

Using Eq. I|2.27|l . one can easily solve Eq. I|2.26|l . We assume, as in the previous section, that the normal metal 
(ferromagnet) is in a contact with the superconductor at x = (x is the coordinate perpendicular to the interface). 
The other boundary of the ferromagnet is located at x — d F and the space at x > d F is empty. 

The proper solution for the diagonal matrix elements f± = A 1(22) can be written as 

f _i_ fees cosh(K £± (x-d F )) n , 
f±{x) = { K =±7F sinh(« E± d F ) U<-X<d F 

x > d F 



Here K, e ± = y/— 2i(e ± h)/Dp is a characteristic wave vector that determines the inverse penetration depth of the 
condensate functions /o,3 into the ferromagnet. 

Usually, the exchange energy h is much larger than the energy e (e cx max{A, T}). This means that the condensate 
penetration depth £p = y/Dp/h is much shorter than the penetration depth into a normal (non- magnetic) metal ^n- 
The strong suppression of the condensate in the ferromagnet is caused by the exchange interaction that tries to align 
the spins of electrons parallel to the magnetization. This effect destroys the Cooper pairs with zero total magnetic 
moment. 
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It is worth mentioning that the condensate function f± experiences oscillations in space. Indeed, for a thick F layer 
{dp » £p) we obtain from Eq. 



/± = ± 



E t K F ±lF 



exp(— x/£,f)[cos(x/1;p) ± ism(x/£p)]- 



(2.29) 



where E e = \/ e 2 — A 2 , Kp± = K e ±(e) at e = 0. The damped oscillations of f± lead to many interesting effects and, 
in particular, to a non-monotonic dependence of the critical temperature on the thickness dp of a F/S bilayer which 
will be discussed in the next section. 

In order to calculate the DOS we have to use the normalization condition, Eq. (|2.7|l . which is also valid for the matrix 

1 + f±, which can be written for small f± as g± w 1 + /±/2. 



(2.30) 



elements f± and g±. Thus, for g± we obtain g± = 
Then the correction to the normalized DOS in the F region bvp =vp — \ takes the form 



Substituting Eq. 



into Eq. 

8vp(dp) 



6v F (x)=Re(fl + fl)/4. 
, we obtain finally the DOS variation at the edge of the F film 

(l/ 4 )Re((^£) 2 [( K£+ sinh( Ke+ d i ,))- 2 + 
+(k £ _ sinh(K £ _d F ))~ 2 ] } . 



(2.31) 



In FIG. |21 we plot the function 8vp(e) for different thicknesses dp and h/A = 20. It is seen that at zero energy e = 
the corr ection to DOS 5vp is positive for F films with d F = 0.8£o while it is negative for films with dp = 0.5£o where 
£o = y/Dp/ A. Such a behavior of the DOS, which is typical for S/F systems, has been observed experimentally 




FIG. 3 Calculated change of the local density of states for a S/F bilayer at the outer F interface. The solid line corresponds to 
a F thickness d F = 0.5£o, where £o = U Df/A, while the dashed one corresponds to d F = 0.8£o- The latter curve is multiplied 
by a factor of 10. 



bv iKontos et all {2001) in a bilayer consisting of a thin PdNi film (5nm < dp < 7.5nm) on the top of a thick 
superconductor. The DOS was determined by tunnelling spectroscopy. This type of dependence of 8vn on d^ can 
also be obtained in N/S contacts but for finite energies e. In the F/S contacts the energy e is shifted, e — > e ± h 
(time-reversal symmetry breaking) and this leads to a non-monotonic dependence of vp on the thickness dF even at 
zero energy. On th e other hand, a non osc illatory behavior of the DOS v(€) has been found recently in experiments on 
Nb/CoFe bilayers l|Revmond et ad feOQOh The discrepancy between the existing theory and the experimental data 
may be due to the small thicknesses of the ferromagnetic layer (0.5nm < dp < 2.5nm) which is comparable with the 
Fermi wave length A_f w 0.3nm. Strictly speaking, in this case the Usadel equation cannot be applied. 

The DOS in F/S structures was studied theoretically in many papers. lHalterman and Vallsl l)2002bf) studied the 
D OS variation numeric a lly for ballistic F/S st r uctures . Th e DOS in quasiballist ic F/S structures was investigated 
bv iBaladie and Buzdir] J200ll) . iBergeret et~al\ l)2002b|) and IZarevan et all l)200ll) and for dirty F/S structures by 
iFazio aruiLucr ^roni| ij 19 991) and|Bi^dirTl|200of) . The subgap in a dirty S/F/N structure was investigated in a recent 
publication bv lGoTubov^ all l|2005j) . 

It is interesting to note that in the ballistic case (rh » 1, r is the momentum relaxation time) the DOS in the 
F layer is constant in the main approximation in the parameter l/(rh) while in the diffusive case (rh << 1) it 
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experiences the damped oscillations. The reason for the constant DOS in the ballistic case is that both the parts of / 
, the symmetric and antisymmetric in the momentum space, contribute to the DOS. Each of them oscillates in space. 
However, while in the diffusive case the antisymmetric part is small, in the ballistic case the contributions of both 
parts to the DOS are equal to each other, but opposite in sign, thus compensating each other. 

Finally, we would like to emphasize that both, the singlet and triplet components, contribute to the DOS. As it is 
seen from Eq.JESJ, the changes of the DOS can be represented in the form 8vf = Rc(/q +/| )/4, which demonstrates 
explicitly this fact. 

2. Transition temperature 

As we have seen previously, the exchange field affects greatly the singlet pairing in conventional superconductors. 
Therefore the critical temperature of the superconducting transition T c is considerably reduced in S/F structures 
with a high interface transparency. 

The critical temperature for S/F bilayer and multilayered structures was calculated in many works 
(iBagrets et a,l\. l2003t iBaladie and BuzdinL I20fl3t i Baladie et all IpOOli iBuzdin and Ku privanovi Il99ll iDemler et all 
I l997t iFominov et all 12001 l2003t iKhusainov and Proshiii Il997t IProshin et all l200li iProshin and Khusainovt 11991 
Il999t iRadovic et adll99ltlTagirovill998tlTollis et adl2005trYou et aZll2004ft. Experimental studies of the T r were also 
repor ted in many publications l)Aartir"erarnT99"7llo^ et all l2006T : iMuhge et oil 

1998j). A good agreement between theory and experiment has been achieved in some cases (see FIG.0J. 
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FIG. 4 Dependence of the superconducting transition temperature on the thickness of the Fe layer as determined by resistivity 
measurements. The das hed line is a fit assuming a perfect interface transparency while the solid line corresp onds to a non-perfect 
interface. Adapted from lLazar. Westerholt. Zabel. Taeirov. Gorvunov. Garif'vanov. and GarifullirJ (§2000). 

One has to mention that, despite of many papers published o n this subject, the prob lem of the transition temperature 
T c in the S/F structures is not completely clear. For example. iJiang et al\ l)l995ft and lQgrin et al\ l)2000ft claimed that 
the non-monotonic dependence of T c on the thickness of the fe rromagnet observed on Gd/Nb samples was due to the 
oscillatory behavior of the condensate function in F. However. I Aarts et al. in an other experiment on V/FeV 

have shown that the interface transparency plays a crucial role in the interpretation of the e xperimental data that 
show ed both non-monotonic and monotonic dependence of T c on ((If)- hi other experiments ({Bourgeois and Dvnesl 
2002) the critical temperature of the bilayer Pb/Ni decreases with increasing dp in a monotonic way. 

From the theoretical point of view the T c problem in a general case cannot be solved exactly. In most papers it 
is assumed that the transition to the superconducting state is of second order, i.e. the order parameter A varies 
continuously from zero to a finite value with decreasing the temperature T. However, generally this is not so. 

Let us consider, for example, a thin S/F bilayer with thicknesses obeying the condition: dp < £f, ds < £s, where 
dF,s are the thicknesses of the F (S) layer. In this case the Usadel equation can be averaged over the thickness (see for 
instance, iBergeret et ai\ (|2001bD ) and reduced to an equation describing an uniform magnetic superconductor with 
an effective exchange energy h and order parameter A. 
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This problem can easily be solved. The Green's functions g± and f± are given by 

9 ± = ^, / ± = |, (2-32) 

where _E e = (e ± h) 2 — A 2 , ft, = rp/i, A = rsA, r> = 1 — r$ = vpd.F / {vpdp + f^ds). In this case the Green's 
functions are uniform in space and have the same form as in a magnetic superconductors or in a superconducting film 
in a parallel magnetic field acting on the spins of electrons. 

The difference between the S/ F bilayer system and a magnetic superconductors is that the effective exchange energy 
h depends on the thickness of the F layer and may be significantly reduced in comparison with its value in a bulk 
ferromagnet. A thin superconducting film in a strong magnetic field H = h/fj,s is an effective Bohr magneton) is 
described by the same Green's functions. The behavior of th ese systems and, in parti c ular, the critical temperature o f 
the su perco n ducting transition T c , was studied long ago bv lFulde and Ferrelil l)l965ft : ILarkin and Ovchinikovl l)l964ft : 
iMakil lll968MSarmal lTl963). It was established that both first and second order phase transitions may occur in these 
systems if h is less or of the order of A. If the effec tive exchange field h exceed s the value A/y/2 « 0.707A, the system 
remains in the normal s tate ( the IClogston ] dl962[) and jCh andrasekharl l|l962ft limit). Independently from each other 
lLarkin and Ovchinikovl l)l964ft and iFuldeanoTl^errelll l)l965ft found that in a clean system and in a narrow interval 
of h the homogeneous state is unstable and an inhomogeneous state with the order parameter varying in space is 
established in the system. This state, denoted as the Fulde-Ferrel-Larkin-Ovchinnikov (LOFF) state, has not been 
observed yet in bulk superconductors. In bilayered S/F systems such a state cannot be realized because of a short 
mean free path. 

In the case of a first order phase transition from the superconducting to the normal state the order parameter A 
drops fro m a finite val ue to zero. The study of this transition requires the use of nonlinear equations for A. It was 
shown bv lTollisI l)2004ft that under some assumptions both the first and second order phase transitions may occur in 
a S/F/S structure. 

In the case of a second order phase transition one can linearize the corresponding equations (the Eilenberger or 
Usadel equation) for the order parameter and use the Ginzburg-Landau expression for the free energy assuming that 
the temperature T is close to the critical temperature T c . Just this case was considered in most papers on this 
topic. The critical temperature of an S/F structure can be found from an equation which is obtained from the 
self-consistency condition Eq. I|2.4|) . In the Matsubara representation it has the form 

ln^ = ( 7 rT*)^( r ^-i/ w /A) ) (2.33) 

where T c is the critical temperature in the absence of the proximity effect and T* is the critical temperature with 
taking into account the proximity effect. 

The function f u is the condensate (Gor'kov) function in the superconductor; it is related to the function /53(e) 
as fs3(iu> n ) = fu, where oj n = ir(2n + 1) is the Matsubara frequency. Strictly speaking, Ea. (|2.33l) is valid for a 
superconducting film with a thickness smaller than the coherence length £5 because in this case / w is almost constant 
in space. 

The quasiclassical Green's function f u obeys the Usadel equation (in the diffusive case) or the more general Eilen- 
berger equation. One of these equations has to be solved by using the boundary conditions at the S/ F interface (or 
S/F interfaces in case of multilayered structures). This problem was solved in different situations in many works 
where an oscillation of T c as a function of the F thickness was predicted (see FIG. |H . In most of thes e papers it was 
assumed that magnetization vectors M in different F layers are collinear. Onlv tFominov et al\ ()2003ft considered the 
case of an arbitrary angle a between the M vectors in two F layers separated by a superconducting layer. 

As mentioned previously, in this case the triplet components with all projections of the spin S of the Coopers pair 
arise in the F/ S/F structure. It was shown that T c depends on a decreasing from a maximum value T cmax at a = 
to a minimum value T cm ; n at a = tt. We will not discuss the problem of T r .for S/ F structures in detail because this 
problem is discussed in other review articles l|Buzdinl . l2005al llzvumov et oH . 12002ft . 



3. The Josephson effect in SFS junctions 

The oscillations of the condensate function in the ferromagnet (see Ea. (|2.29|l ') lead to interesting peculiarities not 
only in the dependence T c (dp) but also in the Josephson effect in the S/F/S junctions. Although, as it has been 
mentioned in the previous section, the experimental results concerning the dependence T c (dp) are still controversial, 
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there is a more clear evidence for these oscillations in experiments on the Josephson current measurements that we 
will discuss here. 

It turns out that under certain conditions the Josephson critical current I c changes its sign and becomes negative. 
In this case the energy of the Josephson coupling Ej = (HI c /e)[l — cosf] has a minimum in the ground state when 
the phase difference ip is equal not to 0, as in ordinary Josephson junc tions, but to it (the so called 7r— junction). 

This effect was predicted for the first time bv lBulaevskii et "all l)l977t) . The authors considered a Josephson junction 
consisting of two superconductors separated by a region c ontaining magnetic impurities. The Josephson current 
through a S/FJ S junction was calculated for the first time bv lBuzdin et al] 11 1982). Di fferen t aspects of the Josephson 
effect in S/F/S structures has been studied in many subsequent papers (Barash e^^j2002; Buzdin and Kuprivanovl 
' 19911 IChtchelkatchev et aP . l200lt iFogelstroml Eo<M iGolubov et ad l2002allHeikkila et adl2000HRadovic et ad 12003 
Zjuzin^^i 20031 e.g). Recen t experiments confirmed the 0-tt transition of the critical current in S/F/S junctions 
l)Bauer et aEl booi iBlum et all. l2002HKontos et all . 120021: IRvazanov efall . l200lHSellier et all\200^ . 

In the experiments of IRvazanov et all (1200 if) a nd iBlum et all ^200^\ . Nb is used as a superconductor and a 
Cu x Ni\- x alloy as a ferromagnet. iKontos et a,l\ (j2002i^ used a more complicated Si/F/I/S structure, where Si 
is a Nb/Al bilayer, S is Nb, I is the insulating AI2O3 layer and F is a thin (40A< g?f < 150A) magnetic layer of a 
PdNi alloy. All these structures exhibit oscillations of the critical current I c . In FIG. the temperature dependence 
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FIG. 5 Measurement of the critical current I c as a functi on of the temperature for a Nb/Cun.48Nip.52/Nb junction. The 
thickn ess of the CuNi layer is cLf = 22nm. Adapted from IRvazanov. Oboznov. Rusanov. Veretennikov. Golubov. and Aartsl 
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of I c measured bv IRvazanov et all l)200l|) is shown. It is seen that the critical current in the junction with aV = 27nm 
turns to zero at T « 2 K, rises again with increasing temperature and reaches a maximum at T 5.5 K. If tem- 
pe rature increases fu rther, I c decreases. In FIG.Elwe also show the dependence of I c on the thickness cLf measured 
bvl Blum et all d2Q02h . The measured oscillatory de pendence is well fitted with the theoretical dependence calculated 
bv lBuzdin et all l |l982|) and lBergeret et The 7r-state in a Josephson junction leads to some observable 
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FIG. 6 Critical current of a Nb/Cu/Ni/Cu/Nb jun ction as a func t ion of the Ni layer thickn ess d. The squa r es are t he measured 
points. The t heoretical fits are done according t olB uzdin et~all il982fl (dashed line) and iBergeret et al ] teOOld) (solid line). 
Adapted from lBlum. A. Tsukernik. and Palevskil i2002f) 

phenomena. As was shown bv iBulaevskii et HH l|l977f) . a spontaneous supe rcurrent may arise in a superconducting 
loop with a ferromagnetic 7r-junction. This current has been measure d bvlBau er et al] l|2004l) . Note also that the 
fractional Shapiro steps in a ferromagnetic 7r-junction were observed bv lSellier et all l)2004|) at temperatures at which 
the critical current I c turns to zero. 

Oscillations of the Josephson critical current I c are related to the oscillatory behavior of the condensate function 
/ in space (see Eq. (|2.29|l ). The critical current I c in a S/F/S junction can easily be obtained once the condensate 
function in the F region is known. We use the following formula for the superconducting current Is in the diffusive 
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limit which follows in the equilibrium case from a general expression (see Appendix |A"|) 

Is = L y L z a F (mT/Ae)J2 T r(Mf+d f+/dx+ f-d f-/dx]), (2.34) 

UJ 

where L y L z is the area of the interface and erp is the conductivity of the F layer. 

In the considered case of a non-zero phase difference the condensate functions f± are matrices in the particle- 
hole space. If in Ea. H2.34J) instead of f± we write a 4x4 matrix for /, then A is given by A = A(iT2 cos(<£>/2) =p 
ifi sin(<yp/2))<73. We set the phase of the right (left) superconductor equal to ±<p/2. For simplicity we assume that 
the overlap between the condensate functions f± induced in the F region by each superconductor is small. This 
assumption is correct in the case (If » £f- Under this assumption the condensate function may be written in the 
form of two independently induced / functions 

f±(x) = (l/CeKF±7F)«T 2 [Aiexp(-K e± (x + d_F/2) 

+ A r exp(-K e ±(-x + d F /2)]- (2.35) 

Here A r ^ is the order parameter in the right (left) superconductor. Substituting Ea. H2.35|) into Ea. (|2.34() . we get 

I s = I c sin(ip) = 4TrT(L y L z )(TF/(KFjF) e:x -P(-dF/t,F) 

A 2 

cos(rf F /£f)^ A2 2 siny>. (2.36) 

When deriving Ea. H2.36|l . it was assumed that the exchange energy h is much larger than both T and A. 
Calculating the sum in Eq. (|2.36|) . we come to the final formula for the critical current 

I c = Atanh(A/2T)cr jF /(K F 7|) exp(— dp/^F) cos(d F /£ F ). (2.37) 

As expected, according to Ea. H2.37[) the critical current oscillates with varying the thickness of the ferromagnet dp- 
The period of these oscillations gives the value of £p an d therefore the valu e of the exchange energy h. For example, 
according to the experiments on Nb/CuNi performed bv lBlum et all l)2002ft h fts HOmeV, which is a quite reasonable 
value for CuNi. 

The non-monotonic dependence of the critical current on temperature observed bv lRvazanov et~al can be 

obtained only in the case of an exchange energy h comparable with A (at least, the ratio h/A should not be too 
large). If the exchange energy were not too l arge, the effectiv e penetration length £,F.eff would be temperature 
dependent. According to estimates presented bv lRvazanov et al\ h ~ 30 K } which means that th e exchange energy in 
this experiment was much smaller than in the one performed bv lBlum et all and bv lKontos et al\ (in the last reference 
h « 35meV). 

The conditions under which the ir— state is realized in S/F/S Josephson junctions of different types were studied 
theoretically in many papers f Buzdin and Baladic, 2003: Buzdin and Kuprivanov, 199l| iBuzdin and Vuiicid 1992; 
IChtchelkatchev et all . l200lt iKrivoruchko and KoshinaL 12001 allLi et all . 12002). In all these papers it was assumed that 
the ferromagnet consisted of a single domain with a mag netization M fixed in s pace. The case of a S/F/S Josephson 
junction with a two domain ferromagnet was analyzed bv lBlanter and Hekkind {2004). The Josephson critical current 
I c was calculated for parallel and anti-parallel magnetization orientations in both ballistic and diffusive limits. It 
turns out that in such a junction the current I c is larger for the anti-parallel orientation. 

A sim ilar effect arises in a S/F/S junction with a rotating in space magnetization, as it was shown bv lBergeret et all 
l|2001c|) . In this case not only the singlet and triplet component with projection S z — 0, but also the triplet component 
with S z = ±1 arises in the ferromagnet. The last component penetrates the ferromagnet over a large length of the 
order of £/y and contributes to the Josephson current. In FIG.[7|the temperature dependence of the critical current is 
presented for different values of Ql, where Q = 2ir/L m , L m is the period of the spatial rotation of the magnetization 
and I is the mean free path. It is seen that at Q = (homogeneous ferromagnet) and low temperatures T the critical 
current I c is negative (n— state), whereas with increasing temperature, I c becomes positive (0— state). If Q increases, 
the interval of negative I c gets narrower and disappears completely at Ql w 0.04, that is, the S/F/S structure with 
a non-homogeneous M is an ordinary Josephson junction with a positive critical current. 

It is interesting to note that the 7r-type Josephson coupling may also be realized in S/N/S junctions provided 
the distribution function of quasiparticles n(e) in the N region deviates significantly from the equilibrium. This 
deviation may be achieved with the aid of a non-equilibrium quasiparticles injection through an additional electrode 
in a multiterminal S/N/S junction. The Josephson current in such a junction is again determined by Ea. H2.34|) in 
which one has to put h = 0, /+ = /_ and replace tanh(e/3) = (1 — 2n(e)) by (l/2)[tanh((e + eV)0) + tanh((e — eV)0)\, 
where V is a voltage difference between N and S electrodes. 
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FIG. 7 Dependence of the critical current on T for hr = 0.06,Aor 
momentum relaxation time. 



0.03, d/l — -K and different values of Ql. Here r is the 



At a certain value of V the critical current changes sign. Thus, there is some analogy between the sign reversal 
effect in a S/F/S junction and the one in a multiterminal S/N/S junction under non-equilibrium conditions. 

Indeed, when calculating Is in a multiterminal S/N/S junction one can shift the energy by eV or —eV. Then 
the function (l/2)[tanh((e + eV)f3) + tanh((e — eV)/3)} is transformed into tanh(e/3) while in the other functions one 
performs the substitution e — > e ± eV. So, we see that eV is analogous to the exchange energy h that appears in the 
case of a S/F/S junction. 

The sign reversal effect in a multiterminal S/N/S jun ction under non-equilibrium condition s has been observed 
bv iBaselmans et all l)l999ft and studied theoretically bv IVolkovl ljl995ft : IWilhelm etaR l)l998|k lYipl (|1998|) . Later 
Hcikkila et al studied theoretically a combined effect of a non-equilibrium quasiparticle distribution on the 

current I c in a S/F/S Josephson junction. 

Conclud i ng thi s Sec tion we note t hat t he experimental results bv lRvazanov et all l|200lf) . iKontos etal\ J2002), 
iBhim et all (|2002ft and IStrunk eTal\ |^94) seem to confirm the theoretical prediction of an oscillating condensate 
function in the ferromagnet and the possibility of switching between the and the 7r-state. 



III. ODD TRIPLET SUPERCONDUCTIVITY IN S/F STRUCTURES 
A. Conventional and unconventional superconductivity 

Since the development of the BCS theory of superconductivity bv iBardeen. Cooper, and Schriefferl l)l957|) . and 
over many years only one type of superconductivity was observed in experiments. This type is characterized by the 
s-wave pairing between the electrons with opposite spin orientations due to the electron-phonon interaction. It can be 
called conventional since it is observed in most superconductors with critical temperature below 20 K (the so-called 
lo w temperature superconduc tors) . 

iBednorz and Miillerl l|l986j) discovered that a La2- x Sr x CuOi compound is a superconductor with a critical tem- 
perature of 30K. This was the first known high-T c copper-oxide (cuprate) superconductor. Nowadays many cuprates 
have been discovered with critical temperatures above the temperature of liquid nitrogen. These superconductors (the 
so called high T c superconductors) show in general a <i-wave symmetry and, like the conventional superconductors, 
are in a singlet state. That is, the order parameter A a p is represented in the form: A a ^ = A • (ia3) a p, where 03 is 
the Pauli matrix in the spin space. The difference between the s and d pairing is due to a different dependence of the 
order parameter A on the Fermi momentum pp = Kkp. In isotropic conventional superconductors A is a k- (almost) 
independent quantity. In anisotropic conventional superconductors A depends on the kp direction but it does not 
change sign as a function of the momentum orientation in space. In high T c superconductors where the d-wave 
pairing occurs, the order parameter A(kp) changes sign at certain points at the Fermi surface. 

On the other hand, the Pauli principle requires the function A(kp) to be an even function of kp, which imposes 
certain restrictions for the dependence of the order parameter on the Fermi-momentum. For example, for <i-pairing 
the order parameter is given by A(kp) = A(0)(fc^ — ky), where k x , y are the components of the kp vector in the Cu—O 
plane. This means that the order parameter may have either positive or negative sign depending on the direction. 

The change of the sign of the order parameter leads to different physical effects. For example, if a Josephson 
junction consists of two high T c superconductors with properly chosen crystallographic orientations, the ground state 
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of the system may correspond to the phase difference <p = tt (it— ju nction). In some high T a s uperconductors the order 
parameter may consist of a mixture of s- and <i-wave components ijTsuei and Kirtlevl 1200 . 

Another type of pairing, the spin-triplet superconductivity, has been discovered in materials with strong electronic 
correlations. The triplet supercond uctivity has been found in h eavy fermion intermetallic compounds and also in 
organic materials (for a review see iMineev and Samokhinl l)l999t) ) . Recently a lot of work has been carried out to 
study the superconducting properties of strontium ruthenate Sr2RuO±. Convincing experimental data have been 
ob tained in favor of th e tri plet, p-wave supercon ductivity. For more details we refer the reader to the review articles 
bv lMaeno et all l)l994|) and lEremin et all (|2004^ . 

Due to the fact that the condensate function < ip a (r,t)ip/3(r' ,t') > must be an odd function with respect to the 
permutations a <-» /3, r <-> r' (for equal times, t = £'), the wave function of a triplet Cooper pair has to be an odd 
function of the orbital momentum, that is, the orbital angular momentum L is an odd number: L — 1 (p-wave), 3 
etc. Thus, the superconducting condensate is sensitive to the presence of impurities. Only the s-wave (L — 0) singlet 
condensate is not sensitive to the scattering by nonmagnetic impurities (Anderson theorem). In contrast, the p-wave 
condensate in an impure material is suppresse d by impuritie s and therefore the order parameter A Q ^ = J2k ^-afsO^F) ~ 
J2k < ''Pair, t)ij)p(r' , t) >k is also suppressed l|Larki That is why the superconductivity in impure S^RuOi 

samples has not been observed. In order to observe the triplet p-wave superconductivity (or another orbital order 
parameter with higher odd L), one needs to use clean samples of appropriate materials. 

At first glance one cannot avoid this fact and there is no hope to see a non-conventional superconductivity in 
impure materials. However, another nontrivial possibility for the triplet pairing exists. The Pauli principle imposes 
restrictions on the correlation function < ip a (r,t)ip(}(r' ,t) >k for equal times. In the Matsubara representation this 
means that the sum ^ w < ip a (r, T)ipg(r', r') >k,u must change sign under the permutation r <-> r' (for the triplet 
pairing the diagonal matrix elements (a — j3) of these correlation functions are not zero). This implies that the sum 
XL < ''Pair, T)ipp(r', t') >k,uj has to be either an odd function of k or just turn to zero. The latter possibility does 
not mean that the pairing must vanish. It can remain finite if the average < ip a (r, T)if>p(r', t') >k,u is an odd function 
of the Matsubara frequency co (in this case it must be an even function of k). Then the sum over all frequencies is 
zero and therefore the Pauli principle for the equal-time correla tion functions is not violated. 

This type of pairing was first suggested bv iBerezinskiil l)l975l) as a possible mechanism of superfluidity in 3 He. He 
assumed that the order parameter A{uj) oc ^ w k < ip a (r, T)ipff(r', r') >k.^ is an odd function of u) : A(uj) = — A(— u>). 
However, experiments on superfluid 3 He have shown that the Berezinskii state is only a hypothetical state and the 
p-pairing in 3 He has different symmetries. As it is known nowadays, the condensate in 3 He is antisymmetric in the 
momentum space and symmetric (triplet) in the spin space. Thus, the Berezinskii hypothetical pairing mechanism 
remained unrealized for few decades. 

However, in recent theoretical works it was found that a superconducting state si milar to the one suggested b y 
Berezinskii might be induced in conventional S/F systems due to the proximity effect l)Bereeret et aZJ . l2001aL 120(33^ . 
In the next sections we will analyze this new type of the superconductivity with the triplet pairing that is odd in the 
frequency and even in the momentum. This pairing is possible not only in the clean limit but also in samples with a 
high impurity concentration. 

It is important to note that, in spite of the similarity, there is a difference between this new superconducting state 
in the S/F structures and that proposed by Berezinskii. In the S/F structures both the singlet and triplet types of 
the condensate / coexist. However, the order parameter A is not equal to zero only in the S region (we assume that 
the superconducting coupling in the F region is zero) and is determined there by the singlet part of the condensate 
only. This contrasts the Berezinskii state where the order parameter A should contains a triplet component. 

Note that attempts to find conditions for the existence of the odd superconductivity were undertaken 
in se v eral papers in connecti o ns wi t h the pairing mechan i sm in high Ij, superconductor s llAhrahams et all 
1993t IBalatskv and Abrahams! Il992t iBalatskv et all. Il995l IColeman etaUL Il993alb! Il995l iHashimotol 1200(1 
Krrck^aTrfok^md^e4rtX1l992t iKirkpatrick^rndBeTlta Il99l|) . In these papers a singlet pairing odd in frequency 
and in the momentum was considered. 

We would like to emphasize that, while theories of unconventional superconductivity are often based on the presence 
of strong correlations where one has to use a phenomenology, the triplet state induced in the S/F structures can be 
studied within the framework of the BCS theory, which is valid in the weak-coupling limit. This fact drastically 
simplifies the problem not only from the theoretical, but also from experimental point of view since well known 
superconductors grown in a controlled way may be used in order to detect the triplet component. 

We summarize the properties of this new type of superconductivity which we speak of as triplet odd superconduc- 
tivity: 

• It contains the triplet component. In particular the components with projection S z — ±1 on the direction of 
the field are insensitive to the presence of an exchange field and therefore long-range proximity effects arise in 
S/F structures. 
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• In the dirty limit it has a s-wave symmetry. The condensate function is even in the momentum p and therefore, 
contrary to other unconventional superconductors, is not destroyed by the presence of non-magnetic impurities. 

• The triplet condensate function is odd in frequency. 

Before we turn to a quantitative analysis let us make the last remark: we assume that in the ferromagnetic regions no 
attractive electron-electron interaction exists, and therefore A = in the _F-regions. The superconducting condensate 
arises in the ferromagnet only due to the proximity effect. This will become more clear later. 

Another type of triplet s uperc onductivity in the S/ F structures that differs from the one considered in this review 
was analyzed bv lEdelsteinl 1)20011) . The author assumed that spin-orbit interaction takes place at the S/F inte rface due 
to a s trong electric field which exists over interatomic distances (the so-called Rashba term in the Hamiltonian (Rashba, 
I960)). It was also assumed that electron-electron interaction is not zero not only in the s-wave singlet channel but 
also in the p-wave triplet channel. The spin-orbit interaction mixes both the triplet and singlet components. Then, 
the triplet component can penetrate into the F region over a large distance. 

However, in contrast to odd superconductivity, the triplet component analyzed bv lEdelsteinl is odd in the momentum 
and therefore must be destroyed by scattering on ordinary nonmagnetic impurities. This type of triplet component was 
also studied in 2-dimensional systems and in S/N s tructures in the presence of the Rashba-type spin-orbit interaction 
l|Edelsteiiill989ll200ltlGor'kov and Rashbal 12001^ . 



B. Odd triplet component (homogeneous magnetization) 

As we have mentioned in section TlLBI even in the case of a homogeneous magnetization the triplet component with 
the zero projection S z = of the total spin on the direction of the magnetic field appears in the S/F structure. Unlike 
the singlet component it is an odd function of the Matsubara frequency u>. In order to see this, we look for a solution 
of the Usadel equation in the Matsubara representation. In this representation the linearized Usadel equation for the 
ferromagnet takes the form 

D F d 2 f F /dx 2 - 2(M«7 - iK& 3 )f F = , (3.1) 

where u> = irT(2n + 1) is the Matsubara frequency and h u — sgn(uj)h. 
The solution of Eq. I|3.1|l corresponding to Eq. 1)2.29(1 can be written as 

f±(uj) = ±(A/i^K ± (uj)j F ) exp(-n±(uj)x) . (3.2) 

where 



«±(w) = y/2(\u\T iK)/D F (3.3) 

and ^ = Vu 2 + A 2 . 

For the amplitudes of the triplet (/o = (/+ + /_)/2) and singlet (fa = (/+ — /_)/2) components we get in the 
ferromagnet 

/ % , . ,„ j. ..exp(—K + (oj)x) , exp(— K-(lo)x), ,„ , s 

/ 3 ,o(w,x) = A/2i£, 7F yy +\> ' ± PV \ J ' ] . (3.4) 

Eqs. (|3.2I) and i|3.4[) show that both the singlet and the triplet component with S z — of the condensate functions decay 
in the ferromagnet on the scale of Rck + having oscillations with Im/? + . Taking into account that k + (uj) = K-(—u>), 
we see that f^(oj) is an even function of u, whereas the amplitude of the triplet component, fo(u)), is an odd function 
of ui. The mixing between the triplet and singlet components is due to the term proportional to h^&s in Ea. H3.1|l . 
This term breaks the time-reversal symmetry. 

Due to the proximity effect the triplet component /o penetrates also into the superconductor and the characteristic 
length of the decay is the coherence length £g. The spatial dependence of this component inside the superconductor 
can be found provided the Usadel equation is linearized with respect to a deviation of the fs matrix from its bulk 
BCS form fscs- In the presence of an exchange field the Green's functions g are 4x4 matrices in the particle-hole 
and spin space. In the case of a homogeneous magnetization they can be represented as a sum of two terms (the f 
matrices operate in the particle-hole space) 



9 = 9h + fih , 



(3.5) 



20 



where g and / are matrices in the spin space. 

In a bulk superconductor these matrices are equal to 

9bcs = 9BCs{u)a \ f BCS = fBCs{u)a z , (3.6) 

where 

9bcs{u) = /bcs(w) = A/i^. (3.7) 

and = V^ 2 + A 2 . 

We linearize now the Usadel equation with respect to a small deviation Sgs = Sgsf^ + SfstT2 = gs ~ 9bcs an d 
obtain for the condensate function Sfs in the superconductor the following equation 

(d 2 /dx 2 )Sf s - 4Sf s = 2i(SA/D s )g% cs a 3 , (3.8) 



where n 2 s = 2^J (w 2 + A 2 ) / Dg and 5 A (x) is a deviation of the superconducting order parameter from its BCS value 
in the bulk. 

A solution for Ea. (|3.8|) determines the triplet component Sfso and a correction J/53 to the singlet component. To 
find the component J/53 is a much more difficult task than to find Sfso because SA(x) is a function of a: and, in its turn, 
is determined by the amplitude J/53. Therefore, the singlet component 5/53 obeys a non-line ar integro-differential 
equation. That is why the critical temperature T r can be calculated only approximat ely |BaCTetf^^z^j ! _|200 i 3i 
i Baladie and Buzdinl l2Q03l iBuzdin and Kuprivanovl 1199(1 iDemler et oil Il997t lizvumov et oil l2002t iRadovic et all 



Baladic and Buzdin, 2003; Buzdin and Km 
19911 iTaeirovl Il998[) . iFominov et al\ <|200 



proposed an analytical trick that reduces the T c problem to a form 
allowing a simple numerical solution. 

On the contrary, the triplet component Sfso proportional to ao can be found exactly (in the linear approximation). 
The solution for Sfso (0) takes the form 

Sfso(x) = <5/5o(0) exp(-Ks(w)x). (3.9) 

The constant Sfso(0) can be found from the boundary condition (see Appendix lAl 

dSf S0 (x)/dx\ x=0 = f FO (0)hs ■ (3.10) 

As follows from this equation, the triplet component in the superconductor Sfso has the same symmetry as the 
component fFo, that is, it is odd in frequency. So, the triplet component of the condensate is inevitably generated 
by the exchange field both in the ferromagnet and superconductor. Both the singlet component and the triplet 
component with 5 = decay fast in the ferromagnet because the exchange field h is usually very large (see Eq. (|3.3|) ). 
At the same time, the triplet component decays much slower in the superconductor because the inverse characteristic 
length of the decay ks is much smaller. 

To illustrate some consequences of the presence of the triplet component in the superconductor, we use the fact 
that the normalization condition g 2 — 1 results in the relation 

9093 = /3/0 (3.11) 
The function go entering Eq. I|3.11|l determines the change of the local DOS 

is(e)=Reg (e) (3.12) 
while the function 53 determines the magnetic moment M z of the itinerant electrons (see Appendix |A"|) 

M z = (i B i>iirT g 3 (u) (3.13) 

We see that the appearance of the triplet component in the superconductor leads to a finite magnetic moment in the 
S'-region, which can be spoken of as an inverse proximity effect. This problem will be discussed in more detail in 
section IY.BI 

Thus, even in the case of a homogeneous magnetization, the triplet component with S z = arises in the S/F 
structure. This fact was overlooked in many papers and has been noticed for the first time bv lBereeret et all (J2003) . 
This component, as well as the singlet one, penetrates the ferromagnet over a short length £f because it consists of 
averages of two operators with opposite spins < V'tV'l > an d is strongly suppressed by the exchange field. The triplet 
component with projections S z = ±1 on the direction of the field results in more interesting properties of the system 
since it is not suppressed by the exchange interaction. It can be generated by a non-homogeneous magnetization as 
we will discuss in the next section. 



21 



C. Triplet odd superconductivity ( inhomogeneous magnetization) 

According to the results of the last section the presence of an exchange field leads to the formation of the triplet 
component of the condensate function. In a homogeneous exchange field, only the component with the projection 
S z = is induced. 

A natural question arises: Can the other components with S z = ±1 be induced? If they could, this would lead to a 
long range penetration of the superconducting correlations into the ferromagnet because these components correspond 
to the correlations of the type < ipl^t > with parallel spins and are not as sensitive to the exchange field as the other 
ones. 

In what follows we analyze some examples of S/ F structures in which all projections of the triplet component are 
induced. The common feature of these structures is that the magnetization is nonhomogeneous. 

In order to determine the structure of the condensate we will use as before the method of quasiclassical Green's 
functions. This allows us to investigate all interesting phenomena except those that are related to quantum interference 
effects. 

The method of the quasiclassical Green's functions can be used at spatial scales much longer than the Fermi wave 
length 3 . As we have mentioned already, in order to describe the S/F structures the Green's functions have to be 4 x 4 
matrices in the particle-hole and spi n space. Such 4x4 matrix Gr een's functions (not necessarily in the quasiclassical 
form) have been used long ago by ijMakil ll969HVaks et a/.l . ll962|) . Equations for the quasiclassical Green's functions 
in the presence of the exchange field similar to the Eilenberger and Usadel equations can be derived in the same way 
as the one used i n the non-magnetic case (see Append ix |A" |) . For example, a generalization of the Eilenberger equation 
was presented bv lBergeret. Efetov. and LarkirJ |2000) and applied to the study of cryptoferromagnetism. 



1. F/S/F trilayer structure 



We start the analysis of the non-homogeneous case by considering the F/S/F system shown in FIG. [HI The structure 
consists of one S layer and two F layers with magnetizations inclined at the angle ±a with respect to the z-axis (in 
the yz plane). 



-a 



a 



-(d s+ d F ) -dg 



d s+ d F x 



FIG. 8 Trilayer geometry. The magnetization of the left (right) F layer makes an angle a (—a) with the z-axis. 



We want to demonstrate now that the triplet component with S z = ±1 inevitably arises due to the overlap of the 
triplet components generated by the ferromagnetic layers in the S layer. It is not difficult to understand why it should 
be so. 

As we have seen in the previous section, each of the layers generates the triplet component with the zero total 
projection of the spin, S z = 0, on the direction of the exchange field. If the magnetic moments of the layers are 
collinear to each other (parallel or antiparallel) , the total projection remains zero. However, if the moments of the 
ferromagnetic layers are not collinear the superposition of the triplet components coming from the different layers 
should have all possible projections of the total spin. 



3 Note that, as was shown by iGalaktionov and ZaikirJ <2002ft ; IShelankov and Ozanal 1200QD . in the ballistic case and in the presence of 
several potential barriers some effects similar to the quantum interference effects may be important. We do not consider purely ballistic 
systems assuming that the impurity scattering is important. In this case the quasicl assical approach is applicable. T he applicability of 
the quasiclassical approximation was discussed long ago by Larkin and Ovchinnikov iLarkin and O vchinnikov. 1963)- 
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From this argument we can expect the generation of the triplet component with all projections of the total spin 
provided the thickness of the S layer is not too large. Since the only relevant length in the superconductors is 
£s w y/D s /nT c , we assume that the thickness of the superconducting layer S does not exceed this length. 

Now we perform explicit calculations that support the qualitative conclusion about the generation of the triplet 
component with all projections of the total spin. We consider the diffusive case when the Usadel equation is applicable. 
This means that the condition 

hr « 1 (3.14) 

is assumed to be fulfilled (t is the elastic scattering time). 

The linearized Usadel equation in the F region takes the form (see Appendix 

dl J - 4/ + ^ {fofo, /]+ cosa ± r 3 [a 2 , /] sin a} = , (3.15) 

where / is a 4 x 4 matrix (condensate function) which is assumed to be small and [173, /]+ = 63./ + /.03. The wave 
vectors k u and entering Eq. l|3.15[l have the form 

kI = 2\u\/D f (3.16) 

and 

k\ = 2hsgmj/D F (3.17) 

The magnetization vector M lies in the (y, z)-plane and has the components: M = A/{0, ± sina, cos a}. The sign 
"+" ("— ") corresponds to the right (left) F film. We consider here the simplest case of a highly transparent S/F 
interface and temperatures close to the critical temperature of the superconducting transition T c . In this case the 
function /, being small, obeys a linear equation similar to Ea. (|3.8|l 

(d 2 f/dx 2 ) - 4/ = 2i(6A/D s )g 2 BCS , (3.18) 

where n 2 s = 2\w\/Dg. 

The boundary conditions at the S/F interfaces are 

fx=d s +o — fx=d s -o (3.19) 
-y(df/dx)\ F = (df/dx)\ s . (3.20) 

where 7 = of/cts and ap (as) is the conductivity in the ferromagnet (superconductor). 

The first condition, Eq. I|3.19|l . corresponds to the continuity of the condensate function at the S/F interface with 
a hig h transparency, whereas Eq H3.2()(l ensures the continuity of the current across the S/F interface l(Volkov et all 
12003) . 

A solution for Eas. (|3. 1513. 181) with the boundary conditions H3.19l3.20f) can easily be found. The matrix / can be 
represented as 

/ = «r 2 (g) / a + in (8 A , (3.21) 

where 

/1 = bi(x)ai, f 2 = b 3 (x)a 3 + b (x)a , (3.22) 

In the left F layer the functions bk (x) are to be replaced by bt (%) ■ For simplicity we assu me that the thickness of the 
F films dp exceeds £p (the case of an arbitrary dp was analyzed bv lBergeret et al\ ((2003) ) . Using the representation, 
Eqs. H3.21I3.22|) . we find the functions bi(x) and bi(x). They are decaying exponential functions and can be written 
as 

6fc(x) = b k exp(-/«(x - d s )), b k (x) = b k exp(«(a; + ds)) (3.23) 

Substituting Ea. H3.23|l into Eas. (|3.15|l - H3.18|l . we obtain a set of linear equations for the coefficients b k . The condition 
for the existence of non-trivial solutions yields an equation for the eigenvalues k. This equation reads 

(k 2_ k 2 )[(k 2_ k 2 )2+k 4 ]=0 (3 ^ 24) 
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Eq. H3.24fl is of the sixth order and therefore has six solutions. Three of these solutions should be discarded because 
the corresponding to bk (x) grow when going away from the interface. The remaining three solutions of Eq. (|3.24|) 
give three different physical values of k. 

If the exchange energy h is sufficiently large (h >> {T, A}), the eigenvalues are 

k = (3.25) 
k ± (l±i)K h (3.26) 

We see that these solutions are completely different. The roots k± proportional to Kh (c/. Eq. I|3.17(l ). are very 
large and therefore the corresponding solutions bk (s) decay very fast (similar to the singlet component). This is the 
solution that exists for a homogeneous magnetization (collinear magnetization vectors). 

In contrast, the value for k given by Eq. (|3.25|) . is much smaller (see Eq. i|3.16[l % ) and corresponds to a slow decay 
of the superconducting correlations. The solutions corresponding to the root given by Ea. H3.25|) describe a long-range 
penetration of the triplet component into the ferromagnetic region. For each root one can easily obtain relations 
between the coefficients bk(x). As a result, we obtain 



b x {x) = b w e- K ^ x - d ^ 

-sin a \b 3+ e- R +^- ds "> - b 3 ^e- K - {x - ds) 



(3.27) 



b (x) = -tana-fc^^"^ 5 ' 

- cosa \b 3+ e~ K+( - x - ds) - b 3 -e- K - {x - ds) ] (3.28) 

and 

b 3 (x) = b 3+ exp(-K+(x - d s )) + b 3+ exp(-K_(x - ds)) (3.29) 

The function b\ (x) is the amplitude of the triplet component penetrating into the F region over a long distance of the 
order of kJ 1 ~ £at. Its value as well as the values of the other functions bk(x) is found from the boundary conditions 
(|3.19I3.20|I at the S/F interfaces. 

What remains to be done is to match the solutions for the superconductor and the ferromagnets at the interfaces 
between them. The solution for the superconductor satisfies Ea. (|3.18|l and can be written as 

/3(a) = A/iE u +a 3 cosh(K S x) (3.30) 
/o(x) = a cosh(K S x) (3.31) 
fi(x) — di sinh(K5x) , (3.32) 

where E w = \J up- + A 2 . 

Matching these solutions with Eqs. (|3.27I3.29|I at the S/F interfaces we obtain the coefficients bk and bk as well 
as ak- Note that b 3 ± = b 3 ± and 6 W = — b w . Although the solution can be found for arbitrary parameters entering the 
equations, we present here for brevity the expressions for 63^ and b u in some limiting cases only. 

Let us consider first the case when the parameter jKh/fts is small, that is, we assume the condition 



JK h /Ks~^J-^«l (3.33) 
Vs V F) s ixL c 

to be fulfilled. 

Here vp g is the density of states in the ferromagnet and superconductor, respectively (in the quasiclassical approx- 
imation the DOS for electrons with spin up and spin down is nearly the same: h << ej?). The condition, Eq. I|3.33|l . 
can be fulfilled in the limit Dp « Ds- Taking, for example, vp sa 1/5, lp w 30 A and Is ~ 300 A, we find that h 
should be smaller than 30T C . 

In this limit the coefficients &i,3± and ai can be written in a rather simple form 

2A jK h sinacos 2 a 

bu^—Eri ) ■ r , on s » 3.34) 

E u k s sinh(29 s ) 

b 3+ w 63- w (3.35) 
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«:i = 



A JK h 1 

'iE u k s sinh(26 s )' 



(3.36) 



where O — Ksds- 

As follows from the first of these equations, Eq. I|3.35|l . the correction to the bulk BCS solution for the singlet 
component is small in this approximation and this justifies our approach. 

At the S/ F interface the amplitude of the triplet component 6 W is small in comparison with the magnitude of the 
singlet one 63+ . However the triplet component decays over a long distance £jv while the singlet one vanishes at 
distances exceeding the short length £p. The amplitudes b u and b^± become comparable if the parameter 7«/ 1 /ks is 
of the order of unity. 

It follows also from Eq. I|3.34|l that the amplitude of the triplet component b^ is zero in the case of collinear vectors 
of magnetization, i.e. at a = or a = tt/2. It reaches the maximum at the angle a m for which sina m = l/y/3. 
Therefore the maximum angle-dependent factor in Eq. I|3.34[) is sin a m cos 2 a m = 2/3V3 « 0.385. 

One can see from Ea. H3.34|) that b w becomes exponentially small if the thickness ds of the S films significantly 
exceeds the coherence length £5 w \J Ds/tzT c . This means that in order to have a considerable penetration of the 
superconducting condensate into the ferromagnet one should not make the superconducting layer too thick. 

On the other hand, if the thickness ds is too small the critical temperature T c is suppressed. In order to avoid this 
suppression o ne has to use, for i nstanc e, an F/S/F structure with a small width of the F films. Similar systems were 
considered bv iBeckmann et all l)2004j) . where non-local effects of Andreev reflections in a S/F nanostructure were 
studied . 

Another limiting case that allows a comparatively simple solution is the limit of small angles a l|Volkov et all . l2003fl 
but an arbitrary parameter jKh/^s, Eq. I|3.33|) . At small angles a the amplitudes of the triplet and singlet component 
are given by the following formulae 



sin a(jKh/ K s) tanhOs 



E u cosh 2 Q s \ tanhOs + (jK h / k s )\ 2 (1 + [iKh/Ks) tanhOg) ' 



(3.37) 



A 



2iEu 1 + (7«±/«s) tanh <d s ' 



(3.38) 



One can see from Eqs. (|3. 3713. 38(1 that, provided the parameter given by Eq. (|3.33() is not small and a, \&s\ ~ 1> the 
amplitudes b^ and 6 3 ± are again comparable with each other. 

The amplitudes of the triplet and singlet components were calculated by iBergeret et all l)2003t ) in a more general 
case of an arbitrary S/F interface transparency and a finite thickness of the F films. 



-(ds+dF) 



-ds ds 



ds+dF 



FIG. 9 The spatial dependence of Im(SC) (dashed line) and the long-range part of Re(TC) (solid line) 
o F I a s = 0.2, h/T c = 50, a F R b /^ F = 0.05, d F ^/T c /Ds = 2, d s ^/Tc/D s = 0.4 and a 



We have chosen 



7r/4. The discontinuity of the TC a t 

the S/F interface is because the short-range part is not shown in this figure. Taken from Bergcrct, Volko v. and Efetovl lH)03) 



In FIG-IHlwe plot the spatial dependence of the triplet (TC) and singlet (SC) components in a F/S/F structure. 
It is seen from this figure that, as expected, the triplet component decays slowly, whereas the singlet component 
decays fast over the short length For this reason, in a multilayered S/F structure with a varying direction of the 
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magnetization vector M and thick F layers (£/, << <1f) a Josephson-like coupling between neighboring layers can 
be realized via the odd triplet component. In this case the in-plane superconductivity is caused by both triplet and 
singlet components. Properties of such S/F multilayered structures will be discussed in the next chapter. 

Let us mention an important fact. The quasiclassical Green's function g(-d) in the diffusive case can be expanded 
in spherical harmonics. In the present approach, only the first two terms of this expansion are taken into account 
such that 

9 = 9sym + 9as cos I? (3.39) 

where $ is the angle between the momentum p and the x-axis, and g as — —lg S ymdg sym /dx is the antisymmetric part 
of g($) and g sy m is the isotropic part of <?($), which does not depend on $. The antisymmetric part of g determines 
the electric current in the system. 

Higher order terms in the expansion of g are small in the diffusive limit and can be neglected. In the case of a weak 
proximity effect the antisymmetric part of the condensate function in the F region can be written as 

f as cos $ -lf 3 a sgmvdf sym /dx cos d . (3.40) 

This expression follows from the fact that go ss — T3 <g) crnsgnu; (corrections to go are proportional to /g). Eq. (|3.40() 
holds for both the singlet and triplet components. 

As we have clarified previously, the symmetric part /o is an odd function of u. Thus, according to Eg. (|3. 3811 the 
antisymmetric part is an even function of ui so that the total condensate function / = /n + fx cos 1? is neither odd nor 
even function of lu. However, in the diffusive limit it is still legitimate to speak about the odd superconductivity since 
the symmetric part is much larger than antisymmetric part of /. 

If the parameter hr is not small, i.e. the system is not diffusive, the symmetric and antisymmetric parts are 
comparable, and one cannot speak of the odd supercond uctivity. All this di stinguishes the superconductivity in S/F 
structures from the odd superconductivity suggested by iBerezinskiil l)l975|) who assumed that the order parameter 
A(w) was an odd function of u>. In our discussion it is assumed that the order parameter A is an uj— independent 
quantity and it is determined by the singlet component of the condensate function / . 

2. Domain wall at the S/F interface 

In the previous section we have seen how the generation of the triplet component takes place. The appearance of 
this component leads to long range effects in a structure where the angle between the directions of magnetization 
in the different layers can be changed experimentally. This is an example of a situation when the long range triplet 
component of the superconducting condensate can be produced under artificial experimental conditions. 

In this section we show that the conditions under which the triplet long range superconducting correlations occur are 
considerably more general. It is well known that the magnetization of any ferromagnet can be quite inhomogeneous 
due to the presence of domain walls. They are especially probable near interfaces between the ferromagnets and 
other materials. Therefore, making an interface between the ferromagnets and superconductors one produces almost 
inevitably domain walls, and one should take special care to get rid off them. In this section we consider a domain 
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FIG. 10 S/F structure with a domain wall in the region < x < w. In this region a = Qx, where Q is the wave vector 
which describes the spiral structure of the domain wall. For x > w it is assumed that the magnetization is homogeneous, i.e. 
a w = a(x > w) — Qw. 

wall like structure and show that it will also lead to the triplet long range correlations. This structure is shown 
schematically in FIG. E3 It consists of a S/F bilayer with a non-homogeneous magnetization in the F layer. We 
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assume for simplicity that the magnetization vector M = M(0, sin a(x), cos a(x)) rotates in the F film starting from 
the S/F interface (x = 0) and the rotation angle has a simple piece-wise x-dependcnce 

, « / Qx, < x < w /q ai\ 

a(x) = < „ (3-41) 

v ' \ a w = Qw, w < x v ' 

This form means that the M vector is aligned parallel to the z-axis at the S/F interface and rotates by the angle 
a w {— a(w)) over the length w (w may be the width of a domain wall). At x > w the orientation of the vector M is 
fixed. 

We calculate the condensate function in the F region and show that it contains the long range triplet component 
(LRTC). As in the preceding section, we assume that the condensate function in the F region is small. The smallness 
of / in this case is either due to a mismatch of the Fermi velocities in the superconductor and ferromagnet or due to 
a possible potential barrier at the S/F interface. In such cases the transparency of the interface is small and only a 
small portion of the superconducting electrons penetrates the ferromagnet. 

Due to the smallness of the transparency of the interface the function / can experience a jump at the S/ F interface, 
which contrasts the preceding case. The boundary condition for the 4x4 matrix / has the same form as in Eq. (|2.27() 

iFdJ = -f s (3.42) 

The function fs on the right-had side is the condensate matrix Green's function in the superconductor that, in the 
limit considered here, should be close to the bulk solution 

fs = fBcsih ® <5- 3 (3.43) 

We have to solve again Ea. H3.15|l with the boundary conditions (I3.42|) . Therefore we assume that the domain wall 
thickness w is larger than the mea n free path I and the condition, Eq. I|3.14|l is fulfilled (dirty lim it). This case was 
analyz ed bv lBereeret et al\ l)2001bj) . Another case of a thin domain wall (w < I) was considered bv lKadigrobov et al\ 
(2001. 

The problem of finding the condensate functions in the case of the magnetization varying continuously in space is 
more difficult than the previous one because the angle a depends now on x. However, one can use a trick that helps to 
solve the problem, namely, we exclude the dependence a{x) introducing a new matrix f n related to / via an unitary 
transformation (a rotation in the particle-hole and spin-space) 

/ = U.f n .U+ (3.44) 

where U = exp(if 3 ® a\a{x)/2). 

Performing this transformation we obtain instead of Eq. i|3.15[) a new equation 

(d 2 x -Q 2 /2)f n - n 2 J n +iK 2 h [a 3 J n } + 
Q 2 

— ^-(^l/no-i) + iQh[vi,d x fn]+ = (3.45) 

Correspondingly, the boundary condition, Eq. 1)3. 42[) . takes the form 

7H(Q/2)if 3 [a 1) /„] + + df n /dx} = -f s (3.46) 

Ea. H3.45fl complemented by this boundary condition has to be solved in the region < x < w. In the region w < x 
one needs to solve Ea. (|3.15|) with Q = 0. Both the solutions should be matched at x = w under the assumption that 
there is no barrier at this point. Therefore, the matrix /„ and its "generalized" derivative should be continuous at 

x = w 

fn \x=w — = fn \x=w+0 (3.47) 
~2^3[^l,fn]+ + d x f n \x=w-0— dxfn U=m>+0 (3.48) 

In this case the solution has the same structure as Eq. l|3.21|) but small changes should be done. The eigenvalues k 
obey the equation 

l(K 2 ~Q 2 - kI) 2 + 4Q 2 K 2]( K 2 _ K 2 j + k 4 [k 2 _ g2 _ K 2j = Q (3^9) 
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where k 1 h are determined in Eas. (|3. 1613. 17(1 . The eigenvalue given by Eq. I|3.25[l changes. Now it is equal to 

4 - Q 2 + 4 . (3-50) 
while the eigenvalues K±, Eq. I|3.26(l . remain unchanged provided the condition 

Q, k u << K h (3.51) 

is fulfilled. 

In the opposite limit of large Q > > Kh, the eigenvalues k± take the form 

k ± = ±iQ[l t iK 2 h /V2Q 2 }, (3.52) 

Thus, in this limit k± is imaginary in the main approximation, which means that the function f n {x) oscillates fast in 
space with the period 2ir/Q. In this case the eigenvalues H3.5()|l change also and have the form 

= < + % ( 3 - 53 ) 

Therefore the limit of a very fast rotating magnetization (k^/Q — > 0) is analogous to the case of a normal metal, 
i.e. when the condensate penetrates the ferromagnet over the length kJ 1 ~ \J Df/2ttT which is the characteristic 
penetration length of the condensate in a S/N system. 

More interesting and realistic is the opposite limit when the condition l|3.51|l is fulfilled and the long-range pene- 
tration of the triplet component into the ferromagnet becomes possible. 

In the limit of large Kh, (Eq. (|3.51|l '). the singlet component penetrates the ferromagnet over a short length of the 
order £/, = while the LRTC penetrates over the length ~ 1/kq. As follows from Eq. (|3.5U|) . this penetration 

length is about 1/Q (provided w/a w is smaller than the length £jv). 

Now let us find the amplitude of the LRTC. The solution for Ea. lj3.45p in the interval < x < w is determined by 
Eas. l|3~2"Tll3~2"2")) with the functions h (x), i = 0, 1, 3 given by the following formulae 

bi(x) = b Q exp(nQx) + b Q exp(-n Q x) (3.54) 

bo(x) = — 63+ exp(— k + x) + &3_ cxp(— K-x) (3.55) 

and 

bs(x) — 63+ exp(— k+x) + &3_ cxp(— k_x) (3.56) 

In the region w < x the solution for the condensate function /„ takes the form 

fn = ih ® aic UJ exp(-K, UJ (x - wj) (3.57) 

where c u is a coefficient that has to be found by matching the solutions at x = w. 
Terms of the order of Q/nh are small and they are omitted now. 

Then we find from the matching conditions at the S/F interface, Eq. (|3.48|l . the following relations for the 
coefficients 

b 3 ± = (3.58) 



and 



-b Q = (Q/ KQ )(b 3+ -b 3 _) (3.59) 



(the parameter 7^ given by Eq. I|3.42|l l 

One can see from the above equations that the condensate function |/| is small provided parameter |7j?/i±| is large. 
It follows from Eq. (|3.59|) that the amplitude of the LRTC, 6q, is not zero only if the magnetization is nonhomogeneous, 
i.e., Q ^0. 

Matching the solutions (|3.54I3.57|I at x = w, we find for the amplitude of the LRTC 

c " = — o ; r 1 [^5 ) 3.60) 

27f k q sinn a w + n u cosh a w \n + |^ReK + 
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where a w = Qw is the total angle of the magnetization rotation. As it has been mentioned, the amplitude of the 
LRTC is an odd function of u>. 

As one can see from the last expression the amplitude c w increases from zero when increasing Q, reaches a maximum 
at Qmax corresponding a certain angle ot max and then exponentially decreases at a w » a max . 

The maximum of c w is achieved at 




At a w — a max the ratio in the square brackets in Ea. (|3.60[) is equal to w 0.68 . This means that the amplitude of the 
LRTC is of the order of the singlet component at the S/F interface. The width w should not be too small because in 
deriving the expression for cq we assumed the condition w » £h- 

In FIG. ^2 we represented the dependence of \c^\ on a w for a fixed w. The spatial dependence of the LRTC and 
the singlet component is shown in FIG. El It is seen that for the parameters chosen the LRTC is larger than the 
singlet component and decays much slower with increasing the distance a;. 




FIG. 11 Dependence of the amplitude of the triplet component on a w = Qw. We have chosen wk„ = 0.01. 




FIG. 12 Spatial dependence of the amplitudes of the singlet (dashed line) and triplet (solid line) components of the condensate 
functi on in the F wire for different values o f a w . Here w = L/5, e = Et and Ii/Et = 400. Et = D/L 2 is the Thouless energy 
(From iBergeret. Volkov. and Efetovl l|2001aM . 



If the magnetization vector M rotates by the angle ir (a domain wall) over a small length w so that Q ~ ir/w » k w , 
then the ratio in brackets in Ea. (|3.6()|l is equal to 

( )w Q/(Qsinh tt) ps 0.087 (3.62) 

kq sinn a w + n u cosh a w 

which shows that the amplitude of the LRTC in this case is smaller than the amplitude of the singlet component. 

We can conclude from this analysis that in order to get a large LRTC, a small total angle of the rotation of the 
magnetization vector is more preferable. 
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The amplitude of the condensate function calculated here enters different physical quantities. In section UlI.DI we 
discuss how the long-range penetration of the triplet component into the ferromagnet affects transport properties of 
F/S structures. 

It is interesting to note that the type of magnetic structure discussed in this section differs drastically from the one i n 
the case of an in-plane rotating magnetization. The latter was considered recently bv lChampel and Esch rig (200 5albf) . 
It was assumed that the magnetization vector Mp was parallel to the S/F interface and rotates; that is, it has the form 
Mf_=Mq{0, sin(Qy), cos(Qy)} (the x-axis is normal to the S/F interface plane). As shown bv lChampel and Eschrid 
l)2005bj) . the odd triplet component arises also in this case but it penetrates into the ferromagnetic region over a short 
distance of the order of £h- 

3. Spin-active Interfaces 

In almost all papers containing discussions of the S/F structures it is assumed that the S/F interface is spin-inactive, 
i.e. the spin of an electron does not change when the electron goes through the interface. 

Although in many cases it is really so , one can imagine another situation when the spin of an electron passing 
through the interface changes. One can consider a region with a domain wall at the interface also as a " spin-active 
interface" provided the width w of the domain wall is very small but the product Qw is of the order unity. As we 
have seen in section UlLC. 21 at such type of interfaces the triplet condensate arises. 

Boundary co nditions at spin-active S/F interfaces for the quasiclassical Green's functions were d erived in a number 
of pub lications ()Kopu et q?l 12004 iMillis et al.l .ll988) and were used in studying different problems. iKulic and EndresI 
(2000 ) employed these b oundary conditions in the study of a system similar to the one shown in FIG. |S| Contrary 
to iBereeret et all l|2003|) . they assumed that the ferromagnets F are insulators so that the condensate does not 
penetrate them. Nevertheless, the calculated critical temperature T c of the supercondu cting transition depends o n 
the mutual orientatio n of the m a gnetiz ation Mp in the ferromagnets. In accordance with lBaladie and Buzdinl J2003); 
Fom inov et al .] l)2002(l : iTagirovl (1998) where metallic ferromagnets were considered in a F/S/F structure, Kulic 
and Endres found that the critical temperature T c was maximal for the antiparallel magnetization orientation. If 
the directions of magnetization vector Mp are perpendicular to each other, a triplet component also arises in the 
superconductor. The authors considered a clean case only, so that the influence of impurity scattering on the triplet 
component rem ained unclear. 

According to iHuertas- Hernando et all l)2002|) a spin-active N/F interface plays an important role in the absolute 
spin-valve effect which can take place in a S/N/F mesoscopic structure. The authors considered a structure with a 
thin normal metal layer (N) and a ferromagnetic insulator F. The DOS variation i n a co nventional superconductor 
w hich is in contact wit h a ferromagnetic insulator was analyzed bv lTokuvasu et~cd\ l)l988f) . 

lEschrig et al\ l)2003|) considered a clean S/F/ S Josephson junction in which the ferromagnet F was a half metal 
so that the electrons with only one spin orientation (say the spin- up | electrons) existed in the ferromagnet. In this 
case only the triplet component corresponding to the condensate function < ijj^^ > may penetrate the ferromagnet. 
Assuming the p-wave triplet condensate function, the authors have calculated the critical Josephson current I c . They 
showed that the tt— state (negative critical current I c ) is possible in this junction. The dc Josephs on effect in 
a junc tion consisting of two superconductors and a spin-active interface between them was analyzed by IFogelstroml 
l|2000D . 

It would be of interest to analyze the influence of impurities on the critical current in such type of Josephson 
junctions because, as we noted, in a clean case the singlet component can penetrate the ferromagnet (not a half 
metal) over a large distance. 

D. Long-range proximity effect 

In the last decade transport properties of mesos copic superconductor/normal metal S/N structu res were inten- 
sively studied (see for example the review articles bv lBeenakkerl l)l997Tj : [Lambert and Raimondil l)l998|) and references 
therein). In the course of these studies many interesting phenomena have been discovered. Among them is a non- 
monotonic voltage and temperature dependence of the conductance in S/N mesoscopic structures, i.e. structures 
whose dimensions are less than the phase coherence length L v and the inelastic scattering length l £ . This means that 
the resistance R of a S/N structure changes non-monotonically when the temperature decreases below the critical 
temperature T c . 

This complicated behavior is due to the fact that there are two contributions to the resistance in such systems: the 
one coming from the S/N interface resistance and the resistance of the normal wir e itself. The experimentally observed 
changes of the resistance can be both positive (SR > 0) and negative (8R < 0) ijQuirion et all . 120021: IShapira et all . 
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2000). The increase or decrease of the resistance R depends, in particular, on the interface resistance Rs/n- If the 
latter is very small, the resistance of the S/N structure is determined mainly by the resistance of the N wire Rjy. 
This resistance decreases with decreasing the temperature T, reaches a minimum at a temperature of the order of the 
Thouless energy Djq / L 2 N , and increases again returning to the value in the normal state Rn(T c ) at low T, where Dpj is 
the diffusion coefficient and L at is the length of the iV wire. This is the so called re-entrant behavior observed in many 
experiments ilCharlat etai\ Il99fit IChien and Chandrasekharl Il999t Ipimoulas et all . Il995t iGubankov and Margolin! 
119791 iPetrashov et a/.lll995HPothier et adll994HShapira et all . \200tfi 

Theoretical explanations for the non-monotonic behavior of the resistance variation as a f unction of the t emper - 
ature T or voltage V in S/N structure s have been presented by Artemenko et all l)l979(k iGolubov et all l)l997f) : 
iNazarov and Stood l|l996|) : IShapfra et all l|2000|k IVolkov et all j|l996l1 l993). Such a variation of the resistance of the 
normal metal wire can be explained in terms of the proximity effect that leads to the penetration of the condensate 
into the N wire. Due to this penet ration there are two types of contributions to the conductance Gn ijGolubov et all . 
119971 IVolkov and Pavlovskiil Il99fi). One of them reduces the DOS in the N wire and therefore reduces the conduc- 
tance Gn- The other term, similar to the Maki-Thompson term l|Golubov et oi.iri997UVolkov and Pavlovskiii ri996). 
leads to an increase of the conductance of the N wire. 

In principle, the magnitude of the conductance variation 6Gn may be comparable with the conductance Gn- So, 
there are no doubts that the proximity effect plays a very important role in many experiments on S/N structures. 

Recently, similar investigations have been carried out also on mesoscopic F/S structures in which ferromagnets (F) 
were used instead of normal (nonmagnetic) metals. According to our previous discussion, the depth of the condensate 
penetration into an impure ferromagnet equals £f = yj hD Jh. This length is extremely short (5 — 50A) for strong 
ferromagnets like Fe or Ni. Therefore one might expect that the influence of the proximity effect on the transport 
properties of such structures should be negligibly small. 

It was a great surprise that experiments carried out recently on F/S structures showed that the resistance 
variation SR were quite visible (varying from about 1 to 10%) when decreasing the temperature be l ow T a 
l)Aumentado and Ohandrasekhall200lUGiroud eit"aHll998tlLawrence and GiordanoLll996albl : lPetrashov et a/l ll999). 
For example in the experiments bv lLawrence and Giordano! I|1996afbf) . where an Sn/Ni structure was studied, the 
effective condensate penetration length estimated from the measured resistance was about 400A. This quantity ex- 
ceeds £f by ord e r of m agnitude. Similar results have been obtained by iGiroud et all 1 199th on Co/Al structures, by 
IPetrashov et all (^999) on a Ni/Al structures and bv lAumentado and Chandrasekhai pOOl) on Ni/Al structures. 

It is worth mentioning that the change of the resistance was bot h positive and negative. In some experi ments the 
variation SRp was rela t ed to a change of the interface resistance llAum gntaxic^aiKiCh andrasekharl 12001^1 , whereas 
in others ijGiroud et all Il998l lLawrence and Giordano! Il996albl : lPetrashovT^rril999j) to the resistance variation of 
the ferromagnetic wire 5Rp. 
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FIG. 13 Reduction of the resistance of a Ni wire attached to a superconductor(Al). Adapted from IPetrashov et all l|l99flt) . 



In F IG.IT3l we show the temperature dependence of the resistance of a Ni wire attached to an Al bank measured 
bv IPetrashov et all i|l999ft . According to estimates of £_f performed in this experiments, the observed SRf is by two 
orders of magnitude larger than it might be expected from the conventional theory of S/F the contacts. Therefore 
these results cannot be explained in terms of the penetration of the singlet component. 

In FIG. El we show similar data from the experiment on Co/Al structures performed bv lGiroud et all (^.998). 
In this experiments a reentrance behavior of SR was observed. In the limit of very low temperatures T — > the 
resistance was even larger than in the normal state. 

The final explanation of this effect remains until now unclear. However, the long range proximity effects considered 
in the previous sections may definitely contribute to the conductance variation. In order to support this point of 
view we analyze qualitatively the changes of the conductance due to the LRTC penetration into the ferromagnet and 
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FIG. 14 Temperature dependence of the resistance of a Co wire attached to a superconductor (Al) measured bv lGiroud et all 
(1998). Note that at low temperatures the authors observed a reentrance behavior. 



demonstrate that the LRTC may lead to the conductance variation comparable with that observed in the experiments. 

However, before presenting these calculations it is reasonable to understand if one can explain the experiments in 
a more s imple way. Actually , the resistance of the SI F structures has been analyzed in many theoretical works. For 
example. iBelzig et all (|200oj) : lde Jong and Beenakkerl l)l99^) : lGolubovl l)l999l) analyzed a ballistic S/F contact. It was 
shown that at zero exchange field (h = 0), the contact conductance Gp/s is twice as large as its conductance Gp/N 
in the normal state (above T c ), as it should be. This agrees with a conductance in a N/S ballistic contact according 
to theoretical predictions. At the same time, it drops to zero at h — Ep, wh ere Ep is the Fe rmi energy. 

The conductance of a diffusive point contact Gp/s has been calculated bv iGolubovl 1^.999) who showed that Gp/s 
was always smaller than the conductance Gp/N in the normal state. In the case of a mixed conductiv ity mechanism 
(partly diffusive and partly ballistic) the conductance Gp/s has been calculated bv lBelzig et all (J2000). According to 
their calculations it can be both larger or smaller than the conductance in the normal state Gp/M- 

The resis t ance flp of a ferromagnetic wire attached to a superconductor was calculated bv lBergeret et all (12002a - ): 
Falko et al. (1999); J e dema et all (jl999|) and let us shortly describe what happens in such a system. 

The proximity effect was neglected in these works but a difference in the conductivities <7-|-j for spin-up and down 
electrons was taken into account. The change of the conductance (or resistance) 6Gp is caused by a different form of 
the distribution functions below and above T c because of Andreev reflections. 

The conductance Gp(T c ) of the F wire in the normal state (T > T c ) is given by the simple expression 

G F (T C ) = G T + G x , (3.63) 

where Gjj = a-^LpA, Lp and A are the length and cross-section area of the F wire. 

This means that the total conductance is the sum of the conductances of the spin- up and down channels. In this 
case not only the electric current but also the spin current is not zero. It turns out that below T c (T < T c ) the 
conductance decreases and at zero temperature it is equal to 

Gf(0) = 4G|Gj/(G| + Gi) (3.64) 

Eq. I|3.64f> shows that the zero-temperature conductance Gp(0) for the system considered is smaller than the normal 
state conductance Gp(T c ). 

It is possible to obtain the explicit formulae not only in the limiting cases, Eq. I|3.63l l3*U4"|l . but also to describe 
the system at arbitrary temperatures. The general formula for the conductance of the F wire can be written as 

G F {T) = G F (0) tanh(A/2T) + G F (T C )(1 - tanh(A/2T)) (3.65) 

Eqs. (|3.63l) and (|3.64|) can be obtained from Eq. H3.65fl by putting A or T to zero. Eas. (|3.63l3.65[) are valid provided 
the length Lp satisfies the condition 

111 < L F < L S o, L in , (3.66) 

where is the mean free path of spin-up and spin-down electrons, while Lso and L>i n are the spin-orbit and inelastic 
relaxation length, respectively. 
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The resistance of mul titerminal S/F structures was calculated by iMelirJ (12001ft : iMelin and Feinberg (2004); 
IMelin and Pevssonl l)2003ft on the basis of the tunnel Hamiltonian metho d. The influence of s u perconducting c o ntacts 
on giant magnetoresistance in multilayered structures was studied bv iTaddei et al\ l(200lft . iTkachov et al\ (,2002) 
studied an enhancement of Andreev reflection at the S/F interface due to inelastic magnon-assisted scattering. 

One can conclude from the works listed above that neglecting the penetration of the LRTC into the F wire an 
increase in the conductance Gf cannot be explained. Therefore, let us discuss the consequences of the LRTC 
penetration into the ferromagnetic wire. In order to avoid the consideration of the S/F interface contribution to the 
total resistance, we consider a cross geometry (see FIG. I15|) and assume that the resistance of the inter face between the 
F wir e and F reservoirs is negligible. Such a geometry was used, for example, in the experiments bv lPetrashov et al\ 
( 1995). The structure under consideration consists of two F wires attached to the F and S reservoirs. We assume that 
there is a significant mismatch between parameters of the superconductor and ferromagnet so that the condensate 
amplitude induced in F is small and is determined by Eas. (|3.37fl or i|3.60[) . 



FIG. 15 The cross geometry used to measured the changes of the resistance of a F wire due to the proximity effect. 

According to our results obtained previously the long range proximity effect is possible provided there is a domain 
wall near the interface between the superconductor and ferromagnet and we assume that this is the case for the setup 
shown in FIG. El Another possibility to generate the triplet condensate would be to attach to the superconductor 
an additional ferromagnet with a non-collinear magnetization. 

The conductance can be found on the basis of a general formula for the current (see for example the book bv lKoprrml 
(|200lft and Appendix |XJ 

I = {l/We)(L y L z )a F Tra ® f 3 o / de[g R d x g K + g K d x g A ] (3.67) 

where a> is the conductivity of the F wire in the normal state. 

The matrix Green's function g K — g R F — F g A is the Keldysh function related to a matrix distribution function 
F. The distribution function consists of two parts, namely, one of them is symmetric with respect to the energy e, 
the other one is antisymmetric in e and determines the dissipative current. 

In the limit of a weak proximity effect the retarded (advanced) Green's function has the form 

3 fl W«±7a®* + f R{A) , (3-68) 

where f R ^ is given by Eas. (|3~57|> or l|3~oT)jl . 

We have to find the conductance of the vertical F wire in FIG. El I n the main approximation the distribution 
function in this F wire is equal to 

F = F Q • t ® CT + F 3 • f 3 ® ct , (3.69) 

where F , 3 = [tanh((e + V)/2T) ± tanh((e - V)/2T)]. 

The distribution function F% symmetric in e determines the current /. The differential conductance Gd = dl/dV 
can be represented as 

G d = Go + SG , (3.70) 

where Go = apLpA is the conductance in the normal state (here we neglect for simplicity the difference between o-\ 
and ai). 

The normalize d correction to the co nductance due to the proximity effect SS(T) = SG/Gq can be found using a 
general formula l)Bereeret et all l2001aft 



SS(T) = (32Ty 1 Tra J de < (f R - f A f > F^e) 



(3.71) 
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where 

F{,(e) = [ costr 2 ((e + eV)/2T) + cosh" 2 ((e - eV)/2T)}/2. 

The angle brackets < ... > denote the average over the length of the ferromagnetic wire between the F (or N) 
reservoirs. The functions f R(A) are given by expressions similar to Ea. (|3.60|l . This formula shows that if T < Dp/L 2 , 
on the order of magnitude SS(T) ~ |/ tr | 2 , where L is the length of the ferromagnetic wire and |/ tr | is the amplitude 
of the triplet component at the S/F interface at a characteristic energy e c h ~ min{T, Dp/L}. According to Eg. 1)3.60(1 
the amplitude of the triplet component is of the order of c\{p^h/Rb) where p is the resistivity of the ferromagnet and 
ci is determined by the factor in the square brackets, that is, by the characteristics of the domain wall. In principle 
the amplitude \ftr\ may be of the order of 1. 

Strictly speaking, both the singlet and triplet components contribute to the conductance. However if the length 
Lp much exceeds the short length £p only the contribution of the LRTC is essential. 
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FIG. 16 The SG(T) dependence. Here 7 = p£ h /Rb- A/E T > 1 and w/L = 0.05 rfrom lBergeret. Volkov. and Efetovl l|2001al^ . 



In FIG. El w e present the temperature dependence of the correction to the conductance SG(T). It is seen that 
with increasing the temperature SGp(T) decreases in a monotonous way. This dependence differs from the re-entrant 
behavior discussed above that occurs in the S/N structures. The reason for this difference is that the time-reversal 
symmetry in S/F structures is broken and this leads to a difference in transport properties. In a S/N system, 
a relation/ fl (e) = f A (e)\ e= o holds and this equality is a consequence the time-reversal symmetry. That is why 
5G(0) = 6G(T C ) = in S/N structures, whereas in a S/F structure f R (e) ^ f A (e)U=o and that is why 6S(T) T =o + 0. 

Although the LRTC may be the reason fo r the enhancement of the conductivity in the S/F structures (this 
possibility was also pointed out in the work bv lOiroud et al\ (|2003)), our understa nding is based on the assumption 
that the magnetic moment is fixed and does not change with the temperature. iDubonos et al\ (|2002^) suggested 
another mechanism based on an assumption about a domain redistribution when the temperature drops below T c . 
The ferromagnetic wires (or strips) used in different experiments may consist of many domains. Their form and 
number depend on the sample geometry and parameters of the system. When the temperature decreases below 
T c , stray magnetic fields excite the Meissner currents in the superconductor attached to the F wire. Therefore the 
demagnetizing factors change, which may lead to a new domain structure. At the same time, the total conductance 
(or resistance) Gp depends on the form and the number of domains. So, one might expect that the conductance 
Gp(T) below T r w ould differ from Gp in the normal state. This idea was supported by measurements carried out by 
IDubonos et all l)2002|) . In this work a structure consisting of a two-dimensional electron gas and five Hall probes was 
used. An F/ S system (Ni+Al disks) was placed on top of this structure. Measuring the Hall voltage, the authors 
were able to probe local magnetic fields around the ferromagnetic disks. They found that these fields really changed 
when the temperature dropped below T c . 

On the other hand, the Meissner currents and, hence, the effect of the redistribution of the domain walls may be 
considerably reduced in wires, as discussed previously. Changing the thickness of the superconducting wires in a 
controlled way and measuring the conductance could help to distinguish experimentally between the contribution to 
the conductivity of the triplet condensate and the effects of the redistribution of the domain walls. 

An experiment in which the domain redistribution was excluded has been performed bv lNueent et all l)2004|) . The 
authors measured the resistance variation of a ferromagnetic wire (Ni 1 _ 2 ,Cu a; ) lowering the temperature T below the 
critical temperature T c of the superconductor (Al or Pb), which was attached to the middle part of the ferromagnetic 
wire. A magnetic field, strong enough to align all domains in the ferromagnet in one direction but not too strong to 
suppress the superconductivity, was applied to the system. Under these conditions a small increase in the resistance 
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(5R/R w 3 • 10 -3 ) was observed when the temperature T drops below T c . The analysis presented above shows that 
the triplet component leads to an increase of the conductance but not in the resistance of the ferromagnetic wire. 
Therefore this particular experiment can hardly be explained in ter ms of the long-range p roximity effect. Perhaps the 
small increase in the resistance of the ferromagnetic wire observed in lNugent et all ((2004) was related to the "kinetic" 
mechanism discussed above (see Eq. 113. 65H) or to weak l ocalization corrections caused by the triplet Cooperons 
l)MnCa.nn et all . 12(10(1) . According to iMcCann et all l|2000l) the change of the resistance of the ferromagnetic wire 
is positive (contrary to the contribution of the LRTC) and its order of magnitude is (e 2 /K)Rf, where Rf is the 
resistance of the F wire in the normal state. In order to clarify the role of the LRTC in the transport properties of 
S/F structures, further theoretical and experimental investigations are needed. Note that strong ferromagnets like 
Fe are not suitable materials for observing the contribution of the LRTC into the conductance variation because of 
the strong exchange field h. In this case, according to Ea. H3.34() and Eq. (|3.(jQ[l . the amplitude of the LRTC is small 
because it contains h in the denominator. 



IV. JOSEPHSON EFFECT IN S/F SYSTEMS (INHOMOGENEOUS MAGNETIZATION) 

As we have mentioned above, one of the most interesting issues in the S/F structures is the possibility of switching 
between the so called 0- and 7r-states in Josephson S/F/S junctions. The 7r-state denotes the state for which the 
Josephson critical current I c becomes negative. This occurs for a certain thickness <1f and temperature T. In this 
state the minimum of the Josephson coupling energy Ej — (HI c /e)(l — cos</>) corresponds to a phase difference of 
<j> = 7r but not to 4> = as in conventional Josephson junctions. 

The reason for the sign reversal of I c is the oscillatory dependence of the condensate functions / on the thickness 
d,F (see Ea. (|2.37[) ). Since the critical current I c is sensitive to the phase of the condensate function at the boundary, 
the 7r-state is a rather natural consequence of the os cillations. 

The possibility of the ir— sta te was predicted by jBulaevskii et all | |l977l) an d iBuzdin et all 1^)82), and studied 
later in many other works (e.g. IBuzdin and Vuiicial [l992TTR^dovice^ all . Il99l|) . Experimentally, this phenomenon 
manifests in a non-monotonic dependence of the critical temperature on the thickness of the junction observed in many 
experiments and discussed in Section lil.B.2l Another manifestation of the transition from the 0- state to the ir- state 
is the sign reversal of the critical current observed in the experiment bv lRvazanov et all l|200ll ) on Nb/Cu x Nii- x /Nb 
Josephson junctions (see FIG.J5J). The proper choice of an alloy with a weak ferromagnetic coupling was crucial for 
the observation of the effect . 

Subsequent experiments, iKontos et all II2002D . iBlum et~aH l)2002|) and iGuichard et all l)2003|) . corroborated the 
observed change of the sign of the Josephson coupling varying the thickness of the intermediate F layer. Qualitatively, 
the experimental data on the Josephson effect in the S/F/S structures are in agreement with the theoretical works 
above mentioned. However, a more accurate control and understanding of the 0-tt transition demands knowledge of 
the magnetic structure of the ferromagnetic materials. 

Almost in all theoretical papers very simplified models of the S/F/S junction were analyzed. For example, 
iBlanter and Hekkind l)2004|) assumed that the F layer consisted either of one domain or two domains with the 
collinear orientation of the magnetization. In this case and according to the discussion of section 1111. CI the LRTC is 
absent in the system. 

If the F layer is a single domain layer, the critical current I c is maximal at a non-zero external magnetic field FL ex ± 
equal to —AitMf, where Mf is the magnet i zation of the F layer. At the same time , in ex periments ijBlum et alJ . l2002t 
IKontos et all . 120021 l200lt iRvazanov etkll . l200lt ISellier et all 12004 IStrunk et all . Il994|) a decrease of the current I c 
with increasing field H ext was observed and it was maximal at H ext = 0. This means, as it was assumed in these 
experimental works, that the F layer in real junctions contains many magnetic domains. In this case the Josephson 
critical current I c may change sign in the S/F/S junctions with a multidomain magnetic structure even if the local 
Josephson current density j c is always positive. The reason for the sign reversal of I c in this c ase is a spatial modulation 
of th e phase difference <fi(x) due to an alternating magnetization M (x) in the domains l|Volkov and Anishchankal 
120041) . In order to determine the mechanism that leads to the sign reversal of the critical current further experiments 
are needed. 

In this chapter we discuss a new phenomenon, namely, how the Josephson coupling between the F layers in the 
S/F structures is affected by the LRTC. 

First, we consider a planar S/F/S Josephson junction with a ferromagnet magnetization Mp rotating in the 
direction normal to the junction plane. This model is an idealization of a real multidomain structure with different 
magnetization orientations. In this case, as we discussed in Section UlI.DI the LRTC arising in the structure affects 
strongly the critical current I c . 

Next, we will analyze a multilayered S/F/S/ ... structure in which the vector Mp has a different direction in 
different F layers. Again, in this case the LRTC arises in the system. Interestingly, if the thickness of the F layers 
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dp is much larger than the penetration length £p of the singlet component but less or of the order of £jy, then the 
Josephson coupling between the F layers is realized due to the LRTC (odd triplet superconductivity in the transverse 
direction). At the same time, the in-plane superconductivity is due to the conventional singlet superconducting 
pairing. 

Finally we will discuss the dc Josephson effect in a SF/I/FS junction (here SF is a superconductor-ferromagnet 
bilayer and / is a thin insulating layer). In this structure, the exchange field may lead not only to a suppression of 
the Josephson coupling as one could naively expect but, under a certain condition, to its enhancement. 

Let us consider first a planar S/F/S Josephson junction. We assume the following spacial dependence of the 
magnetization vector in the F layer: Mj? = Mp(0, sin(Qx), cos(Qx)), where the x-axis is normal to the junction 
plane. 

In this case, as we have seen in Section IIII.CJ.2I the LRTC arises. Due to the long range penetration into the 
ferromagnet the triplet component can give a very important contribution to the Josephson current. A general 
expression for the Josephson current can be written in the form 

I j = (L y L z /4e)a F (TrT)Tr(& Q ® f 3 . ^ f u dj u ) . (4.1) 

We assume that the impurity concentration is sufficiently high and therefore the condensate function f u should be 
found from the Usadel equation. In the limit of a weak proximity effect (the S/F interface transparency is not too 
high) this equation can be linearized and solved exactly. The solution for the f u matrix in the F region can be found 
in a similar way as it was done in Section IIII.CJ.2I Due to the rotation of the ma gnetization the condensate function 
contains the LRTC. We obtain for the Josephson current ijBergeret et a/.t feOQlcj) the following expression 

Ij = I c sin <j> (4.2) 

where the critical current I c is equal to 

I c = (L y L z a F /l)j 2 F ReY,f! 

and nf = 2\uj n \/ A + Q 2 . The parameter jp = (3/4) < fiT(fi) > is an effective, averaged over angles, transmittance 
coefficient which characterizes the S/F interface transparency and k+ is defined in Eq. (|3.52|) 

The first term in the brackets containing the parameter K + corresponds to Eq. (|2.36() . It decays by increasing 
the thickness dp over the short characteristic length £p = y/Dp/h and can change the sign. The second term in Eq. 
(|4.3() originates from the rotation of h along the x-axis. It decays with the thickness dp over another characteristic 
length that can be much larger than £p. Therefore this term results in a drastic change of the critical current. 

The presence of the second term in Eq. (|4.3|l is especially interesting in the case when the thickness dp of the 
ferromagnetic spacer between the superconductors obeys £p < dp < n^ 1 . Then the main contribution to the 
Josephson coupling comes from the long-range triplet component of the condensate. Another important feature 
of this limit is that for sufficiently large values of Ql the critical current is always positive (no possibility for the 
7r-contact). This can be seen from FIG.[7| 

The fact that the superconductivity looses its "exotic properties" at large Q is quite understandable. The super- 
conductivity is sensitive not to the local values of the exchange field but to its average on the scales of the order of 
the superconducting coherence length. If the exchange field oscillates very fast such that the period of the oscillations 
is much smaller than the superconducting coherence length, its average on this scale vanishes and therefore all new 
properties of the superconductivity originating from the presence of the exchange field become negligible. 

To conclude this introduction we summarize the results known for S/F/S Josephson junctions. 

When the magnetization in the ferromagnetic F is homogeneous, we have to distinguish between two different 
cases. 

In the dirty limit (hr <C I) the change of the sign of the critical current occurs if the thickness of the F layer d F is 
of the order of £p. The condensate function in the F layer decays exponentially over this ^ and oscillates with the 
same period. 

In the opposite clean limit, hr 3> 1, the condensate function oscillates in space with the period vp/h and decays 
exponentially over the mean free path /. 

Finally, if the ferromagnetic region contains a domain wall described by a vector Q, a long-range component of 
the condensate appears. It decays in the F film over a considerably larger length of the order £jy = \JD/2-kT that 
can greater exceed the characteristic length (~ D/h) in a homogeneous F layer [Q = 0). In this case the coupling 
between the superconductors survives even if the thickness of F is larger than £p. 

It is clear that the presence of a domain wall between the superconductors is something that cannot be controlled 
very well experimentally. Therefore in the next section we discuss a possible experiment on S/ F multilayered structures 
that may help in detecting the LRTC by measuring the Josephson critical current. 
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A. Josephson coupling between S layers via the triplet component 

In this subsection we analyze another type of multilaycred S/F structure in which the LRTC arises. This is a 
multilayered periodic ...S/ F n _i/ S/ F n /S/ F n+ i/ S... structure with alternating magnetization vector M^„ in different 
F layers. We assume that the vector Mj? in is turned with respect to the vector M^„_i by an angle 2a, such that 
the angle increases monotonously with increasing n. We call this arrangement of the magnetization the one with a 
positive chirality. 

In an infinite system the magnetization vector Mi? averaged over n is equal to zero (it rotates when one moves 
from the n — th to the (n + 1) — th, layer etc.) . Another type of chirality (negative chirality) is the arrangement when 
the angle between vectors Mir n and Mf,„_i is equal to 2a(— l) n . In this case the averaged vector Mj? is not zero. 

In Section Till. C. II we have seen that in a F/S/F structure with a non-collinear orientation of the magnetization 
vectors in the F layers the LRTC arises. If one assumes that the thickness of the F layers dp is larger than the coherence 
length in the normal metal £/v, the overlap of the condensate functions created in a F layer by neighboring S layers 
is weak, and the solutions given by Eos. (|3. 2113. 29fl remain valid for the multilayered S/F structure. 

Using these solutions one can calculate the Josephson current between neighboring S layers. As the thickness dp 
is assumed to be much larger than £f (as usual, we assume that £f << £,n), the Josephson coupling between the S 
layers is solely due to the LRTC. So, in such systems we come to a new type of the superconductivity: an odd triplet 
out-of-plane superconductivity and the conventional singlet in-plane superconductivity (the triplet component gives 
only a small contribution to the in-plane superconductivity). 

Using the general Eqs. I|3.2H3.29|) one can perform explicit calculation s for this case w i thout considerable difficulties. 
As a result, the Josephson critical current I c can be written as follows l|Bergeret et q/J.l2003|) 

eR F Ic = ±2-KT^K UJ d F b\{a) (l + tan 2 a) e - rf ^ , (4.4) 



where 



. k 2 s (k + - K_)sgnw 

bi(a) = -/Bcssma- 



cosh" 9 5 (A/+T_ + M_T+) (g BC s + If^u tanh8 F ) 

6 S = K s d s , Of = Kud F , k± = K±/(g B cs + Jfk±), k = K u /(g B cs + 1f^ tanh 8f), k s = ks/{9bcsi) and 

M± = T±(ks coth Qs + k tanh Op) + tan 2 a C±(ks tanh 8s + k tanh Op) 
T± — Kg tanh 0£ + k± 
C± = ks coth 85 + k± . 

Rp is defined as Rp — 2dp j {L v L z ap\ Eq. (|4.4|l describes the layered systems with both the positive ("+" sign) and 
negative ("-" sign) chirality. 

One can see from Eq. (|4.4(l that in the case of positive chirality the critical current is positive, while if the 
chirality is negative the system is in the 7r-state (negative current). This means that changing the configuration of 
the magnetization, one can switch between the and 7r state. 

It is important to e mphasize that th e nature of the 7r-cont act obtained here differs from that predicted by 
iBulaevskii et al\ l)l977i) and observed bv iRvazanov et al\ l)200l[) . In the latter cases the transition is due to the 
change of the values of either the exchange field, the temperature or the thickness of the F film. In the case consid- 
ered in this section, the negative Josephson coupling originates from the presence of the triplet component and can 
be realized in S/F structures with negative chirality. Since for the positive chirality the Josephson current is positive, 
the result obtained gives an unique opportunity to switch experimentally between the and 7r-states by changing the 
angles of the mutual magnetization of the layers. 

A similar dependence of the Josephson current I c on the chirality was predicted i n a Josephson junction S m IS m (I is 
an insulator) between two magnetic superconductors S m bv lKulic and Kulid For the magnetic superconductors 

considered in that work, the magnetization vector M rotated with the angle of rotation equal to a — xQ ■ n^, where 
Q is the wave vector of the x- dependence of the a ngle a, n T is the unit vector normal to the insulating layer I. 
Therefore the chirality (or spiral helicity, in terms of Ku lic and KulicT) in this case is determined by the sign of the 
product Q/f-Qi, where Ql,r is the wave vector in the left (right) magnetic superconductor. 

However, there is an essential difference between the multilayered S/F structure discussed here and the magnetic 
superconductors. In the magnetic superconductors with a spiral magnetization the triplet component also exists but, 
in contrast to the S/F structures, the singlet and triplet components cannot be separated. In particular, in the case 
of a collinear alignment of M, the Josephson coupling in the S/F structures with thick F layers disappears, whereas 
it remains finite in the S m IS m system. 
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FIG. 17 FIG. 16: Dependence of th e critical current on the angle a. The value of the current is given in arbitrary units (from 
iBergeret. Volkov. and Efetovl J2003D ). 

FIG. 1171 shows the dependence of the Josephson current I c on the angle a given by Eq. I|4.4|l . If the mutual 
orientation of M is parallel (a — 0) or antiparallel (2a = n) the amplitude of the triplet component is zero and 
therefore there is no coupling between the neighboring S layers, i.e. I c = 0. For any other angles between the 
magnetizations the amplitude of the triplet component is finite and this leads to a non-zero critical current. At 
2a = 7r/2 ( perpendicular orientation of M) the Josephson current I c reaches its maximum value. 

Another possible experiment for detecting the long ra nge triplet compone nt is the measurement of the density of 
states in the F/S/F system as it is shown in FIG. ^] Kontos et al\ l|200lJ) determined the spatial changes of the 
DOS in a PdNi/Al structure with the help of planar tunnelling spectroscopy. This method could also be used in 
order to detect the LRTC. It is clear that if the thickness of the F layer in FIG. ^] is larger than the penetration of 
the short-range components, then any change of the DOS at the outer boundary of the F layer may occur only due 
to the long range penetration of the triplet component. If the magnetizations of both F layer are collinear no effect 
is expected to be observed, while for a non-collinear magnetization a change of the DOS should be seen. 



FIG. 18 Schematic : Measurement of the change of the density of states at the outer F interface by tunnelling spectroscopy. 
iKontos et al\ i200lf) performed such experiments on S/F structures. 



B. Enhancement of the critical Josephson current 

Another interesting effect in the S/ F structures that we would like to discuss is the enhancement of the Josephson 
critical current by the exchange field. The common wisdom is that any exchange field should reduce or destroy the 
singlet superconductivity. In the previous sections we have seen that this is not always so and the superconductivity 
can survive in the presence of a strong exchange field. But still, it is not so simple to imagine that the superconducting 
properties can be enhanced by the exchange field. 

Surprisingly, this possibility exists and we will demonstrate now how this unusual phenomenon occurs. Although the 
LRTC is not essential to get the critical current enhancement, the short-range triplet component arises in this case and 
it plays a certain role in this effect. The enhancement of the Jos ephson current i n the S /F/I/F/S tunnel structures 
(J stands for an ins ulating layer, se e FIG. |H))| was predicted bv lBergeret et al\ l)2001b|) and further considered in a 
subsequent work by iGolubov et al\ (|2002bTl . As we will see below, if the temperature is low enough and the S/ F 
interface transparency is good, one can expect an enhancement of the critical current with increasing the exchange 
field provided the magnetizations of the F layers are antiparallel to each other. This surprising result can be obtained 
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FIG. 19 The S/F/I/F/S system. I is an insulating thin layer. The relative magnetization of the F layers can be switched. 



in a quite simple way in the limit when the thicknesses ds and d F of the S and F layers are smaller than the 
superconducting coherence length £g ~ y/ D /2ttT c and the penetration length of the condensate into the ferromagnet 

£f ~ y/D/h, respectively. In this case one can assume that the quasiclassical Green's functions does not depend on 
the space coordinates and, in particular, the superconducting order parameter A is a constant in space. Moreover, 
instead of considering the dependence of the exchange field h on the coordinates one can replace it by a homogeneous 
effective exchange field h e ff with a reduced value. Therefore we use in our calculations effective fields A e ff and h e /f 
denned as 

A e ff/A = v s ds{vsds + v F d F )~ 1 , (4.5) 
h eff/ h = v f&f {"sds + vpdp)' 1 , (4.6) 

where i>s and vp are the densities of states in the superconductor and ferromagnet, respectively. 

With this simplification, we can write the Gor'kov equations for the normal and anomalous Green's functions in 
the spin-space as 

(iw + i - ah) G u + AF+ = 1 (4.7) 
(-iuj + £ - ah) F u + AG^ = 0, (4.8) 

where a = [a±, a 2 , 1J3) are the Pauli matrices and £ = e (p) — ep, ep is the Fermi energy, e (p) is the spectrum, and 
oj = (2n + 1) ttT are Matsubara frequencies. (We omit the subscript eff in Eqs. H4.7I4.8(I and below). 

In order to calculate the Josephson current Ij through the tunnel junction represented in FIG. IT^l we use the well 
known standard formula 

/./ = (2nT/eR)TrY,Mh 1 )f u (h 2 )sm<p , (4.9) 

n 

where 




(4.10) 



is the quasiclassical anomalous Green's function, ip is the phase difference between both the superconductors, R is 
the barrier resistance and hyp are the exchange fields of the left and the right f-layers. 

The only difference between Eqs. Q4.9I I4.10[) and the corresponding equations in the absence of the exchange field 
is the dependence of the condensate function / w on h. This dependence can be found immediately from Eqs. I|4.7I 
Oil . 

U = A ((w + iahf + A 2 ^) (4.11) 

What remains to be done is to insert the condensate function / into Eq. (|4.9|l for certain exchange fields hi and h 2 
and to calculate the sum over the Matsubara frequencies to. Although it is possible to carry out these calculations for 
arbitrary vectors hi and /12, we restrict our consideration by the cases when the absolute values the magnetizations 
hi and h 2 are equal but the magnetizations are either parallel or antiparallel to each other. This simplifies the 
computation of the Josephson current but, at the same time, captures the essential physics of the phenomenon. 
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Using Eqs. i|4.9l4.11 || and assuming; fir st that hi and /12 are parallel to each other we write the expression for the 
critical current as (Bcre eret etaR l2001b|) 



I j = I c sin if 



(4.12) 



A 2 (T) 4ttT 



uj 2 + A 2 (T, h) - h 2 



eR ^ (to 2 + A 2 (T,h)~h 2 ) 2 +Au 2 h 21 
The corresponding equation for the antiparallcl configuration is different from Eq. (|4.13|) and can be written as 



(4.13) 



(o)= A 2 (T)4nT 
eR 



u y] (uj 2 + A 2 (T, ft) - /i 2 ) 2 + 4w 2 /i 2 



(4.14) 



One can easily check that the critical current Ic P ^ for the parallel configuration, Eq. 14. 13L is always smaller than 
the current Ie for the antiparallel case. These two expressions are equal to each other only in the absence of any 
magnetization. 
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FIG. 20 Dependence of the normalized critical current on h for different temperatures in the case of a parallel orientation. 
Here eV c = eRI c , /if is the effective exchange field, t — T/Ao and Ao is the superconducting order parameter at T = and 
h = 
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FIG. 21 The same dependence as in FIG. I2()| in the case of an antiparallel orientation. 

In FIGS. 1201 and [21 we represent the dependence of the critical current on the strength of the exchange field. We 
sec from FIG. |201 that for the parallel configuration the exchange field reduces the value of the Josephson current and 
this is exactly what one could expect. At the same time, the critical current grows with the exchange field for the 
antiparallel configuration at low temperatures, which is a new intriguing effect (see FIG. I21|) . 

This unexpected result can be understood from Eq. I|4.14|l rather easily without making numerics. In the limit 
T — > the sum over the Matsubara frequencies can be replaced by an integral and one can take for the superconducting 
order parameter A the values A = An if h < An, and A = i f h > A n, where Ao is the BCS order parameter in the 
absence of an exchange field (see, e.g. lLarkin and Ovchinikovl l)l964j) '). 



40 



Inserting this solution into Eq. (|4.14() one can see that the Josephson critical current I c a ' grows with increasing 
exchange field. Moreover, formally it diverges logarithmical when h — ► A 

/W(/ l ^A )^^l n (A /wo) , (4.15) 

TT 

where I c (0) is the critical current in the absence of the magnetic moment at T — 0, and u>q is a parameter needed to 
cut the logarithm at low energies. 

When deriving Eqs. l|4.13ll4~T4*|l the conventional singlet superconducting pairing was assumed. The electrons of a 
Cooper pair have the opposite spins. This picture of a superconducting pairs with the opposite spins of the electrons 
helps in the understanding of the effect. 

If the magnetic moments in both the magnetic layers are parallel to each other, they serve as an obstacle for the 
Cooper pair because the pairs located in the region of the ferromagnet demand more energy. This leads to a reduction 
of the Josephson current. However, if the exchange fields of the different layers are antiparallel, they may favor the 
location of the Cooper pairs in the vicinity of the Josephson junction. A certain probability exists that one of the 
electrons of the pair is located in one layer, whereas the other is in the second layer. Such a possibility is energetically 
favorable because the spins of the electrons of the pair can now have the same direction as the exchange fields of the 
layers. Then it is more probable for the pairs to be near the junction even in comparison with a junction without 
exchange fields and, as a result, the critical current may increase. 

The results presented above have been obtained for the SF/I/FS structure bv lBergeret et all l)2001bl) . Earlier a 
for mula for the Josephson critical current in the S m IS m (S m is the magnetic superconductor) junction was presented 
bv lKulic and Kuiicl l)200lh . From that formulae one could, in principle, derive an enhancement of the critical current 
for the antiparallel M orientation in magnetic superconductors S m . Unfortunately, the authors seem to have missed 
this interesting effect. 

Some remarks should be made at this point: 

1) The results presented above are valid in the tunnelling regime, i.e. when the transparency of the insulating 
barrier I is low enough. IColubov et al\ (|2002bll have shown that a smearing of the singularity of Ic^ is provided by 
a finite temperature or a not very low barrier transparency. The maximum of the critical current for the antiparallel 
configuration Ic^ decreases with decreasing resistance of the I layer. The effect becomes weaker as the thickness of 
the F layer grows. 

2) We assumed that the S/F interface was perfect. In a structure with a large S/F interface resistance Rs/f the 
bulk properties of the S film are not considerably influenced by the proximity of the F film (to be more precise, the 
condition Rs/f > (^f<^f /^sds)pFS,F must be satisfied, where pf is the specific resistance of the F film). Then, as one 
can readily show (see section llLB|) . a minigap e^F = (Dp) F / (2Rs/f<1f) arises in the F layer. The Green's functions 
in the F layers have the same form as before with A replaced by e&F- The singularity in I c (h) first occurs at h equal 
to e bF . 

A physical explanation for the singular behavior of the critical current 1^ was given bv lColubov et al\ |2002b) 
These authors noticed that the density of states in the F layer has a singularity when h = e^p- At this value of h the 
maximum of Ic is achieved due to an overlap of two eT 1 ! 2 singularities. This leads to the logarithmic divergency of 
the critical current in the limit T — > in analogy with the well known Riedel peak in SIS tunnel junctions for the 
v oltage difference 2A. In the latter case the shift of the energy is due to the electric potential. 

iGolubov et d\ (|2002b) have also shown that for the parallel configuration, at h = e^F t he critical current changes 
its signs, i.e. there is a transition from to a tt junction. Similar results were obtained bv iKrivoruchko and Koshinal 
(2001a b). The c ase of an arbitrary S/F tr ansparency was also studied bv lBarash et al\ l)2002|) : Iclitehelklrichev^^/] 
l|2002ft : iLi etaR (2002). In the paper bv iBarash et al\ l)2002|) the authors calculated the Josephson current as a 
function of the angle between the magnetizations in the F film. 

V. REDUCTION OF THE MAGNETIZATION DUE TO SUPERCONDUCTIVITY: INVERSE PROXIMITY EFFECT 

Until now we have been studying the superconducting properties of different S/F structures for a fixed magnetiza- 
tion. This means that we assumed a certain value for this quantity and its dependence on coordinates. The implied 
justification of this assumption was that the ferromagnetism is a stronger phenomenon than the superconductivity 
and the magnetic moment of conventional fcrromagnets can hardly be affected by the superconductivity. 

This assumption is certainly correct in many cases but not always. Often the presence of the superconductivity can 
drastically change magnetic properties of the ferr oma gnets even if they are s trong. 

Experiments performed bv lMuhge et all l|l998j) and iGarifullin et a^J l|2002|) showed that the total magnetization of 
certain S/ F bilayers with strong ferromagnets decreased with lowering the temperature below the critical supercon- 
ducting transition temperature T c . As an explanation, it was suggested that due to the proximity effect domains with 
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different magnetization appeared in the magnetic materials a nd this could reduce the tota l magnetization. At the 
same time, quantitative estimates based on an existing theory ijBuzdin and Bulaevskiiri988|) led to a conclusion that 
this mechanism was not very probable. 

In this Chapter we address the problem of the reduction of the magnetic moment by the presence of a superconductor 
assuming again that, in the absence of the ferromagnet, we would have the conventional singlet superconducting 
pairing. It turns out that two different and independent mechanisms that lead to a decrease of the magnetization in 
S/F heterostructures due to the proximity effect exist and we give a detailed account of them. 

In order to study the magnetic properties we have to choose a model. One can distinguish two different types of the 
ferromagnetism: a) itinerant ferromagnetism due to the spin ordering of free electrons and b) ferromagnetism caused 
by localized spins. Most of ferromagnetic metals show both of the types of ferromagnetism simultaneously, i.e. their 
magnetization consists of both contributions. 

We consider a model in which the conducting electrons interact with the localized moments via an effective exchange 
interaction. The corresponding term in the Hamiltonian is taken in the form (see Appendix lA"|l : 



This term is suitable to describe s — d or s — / interaction between the s and localized d and / electrons. We also 
consider the ferromagnetic interaction between the localized moments. This interaction can be very complicated and 
to determine it, one should know the detailed band structure of the metal as well as different parameters. However, 
all these details are not important for us and we write the interaction between the localized spins phenomenologically 
as 



It is assumed that J is positive. This interaction, Eq. 1)5. 2fl , is responsible for the ferromagnetic alignment of the 
localized moments and is known as the Heisenberg Hamiltonian. 

So , we consider a metallic ferromagnet in which the conduction electrons interact with localized magnetic moments. 
The ferromagnetic interaction l|5.2[) assures a finite magnetic moment of the background. The total magnetization is 
the sum of the background magnetization (localized moments) and the magnetization of the polarized free electrons. 

In the next two sections we discuss the two different mechanisms that lead to a decrease of the magnetization at 
low temperatures. In Section IV.AI we consider a possibility of changing the magnetic order of the localized magnetic 
spins in a F film deposited on top of a bulk superconductor. The contribution from free electrons to the magnetization 
is first assumed to be small. We will see that for not too strong ferromagnetic coupling J the proximity effect may 
lead to an inhomogeneous magnetic state. Contrary to this case, we consider in Section lV.BI an itinerant ferromagnet 
in which the main contribution to the magnetization is due to free electrons. We will show that the magnetization of 
free electrons may decrease at low temperatures due to a so me kind of spin s c reening. Thus, both effec ts may lead 
to the decrease in the magnetization observed in experiments (jGarifullin et all l2002t iMiihee et a0 . ll998|) . 

A. Cryptoferromagnetic state 

In ll959irAnderson and Suhl| suggested that su perconductivity could coexist with a nonhomogeneous magnetic order 
in some type of materials. lAnderson and Suhll called this state cryptoferromagnetic state. 

The reason for this coexistence is that, if the magnetization direction varies over a scale smaller than the super- 
conducting coherence length, the superconductivity may survive despite the ferromagnetic background. This is due 
to the fact that the superconductivity is sensitive to the ferromagnetic moment averaged on the scale of the size of 
Cooper p airs rather than to its lo cal values. 

In 1988 lBuzd in and Bulaevskiil discussed properties of a bilayer system consisting of a conventional superconductor 
in contact with a ferromagnet. They have shown that the magnetic ordering in the magnet might take the form of 
a structure consisting of small size domains, such that the superconductivity is not destroyed. Of course, as follows 
from Eq. (|5.2() . the formation of a domain- like structure costs a magnetic energy but this is compensated by the 
energy of the superconductor that would have been lost if the magnetic order remained ferromagnetic. 

This is only possible if the stiffness of the magnetic order parameter {J) is not too large. For instanc e this 
nonhomogeneous magnetization occurs in magnetic superconductors as those studied bv lBulaevskii et al\ (^985). But 
can one see it in the heterostructures containing strong ferromagnets like Fe or Ni in contact with conventional 
superconductors? 

At first glance, it seems impossible, since the Curie temperature of, for example, iron is hundred times or more 
larger than the critical temperature of a conventional superconductor. Therefore any change of the ferromagnetic 




(5.1) 




(5.2) 
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order look much less favorable energetically than the destruction of the superconductivity in the vicinity of the S/F 
interface. 

This simpl e argument was howev er questioned in the experiments performed bv lMuhge et al\ (|1998^I on Fe/Nb bi- 
layers and bv lGarifullin et all i|2002|) on V / Pd\- x Fe x structures. Direct measurements of the ferromagnetic resonance 
has shown that in several samples with thin ferromagnetic layers the average magnetic moment started to decrease 
below the superconducting transition temperature T c . 

Of course, one can reduce the influence of the f erromagnet on the supe r condu c tor by dimini s hing t he thickness of 
the ferromagnet. Using the formulae obtained bv lBuzdin and Bulaevskiil (|1988f ) . iMiihge et all 1)19981) estimated the 
thickness of the ferromagnet for which the superconductivity was still possible and got a value of the order of lA, 
which created a doubt on the explanation of the experim ent in this way. 

At the same time, the use of the formulae derived bv lBuzdin and Bulaevskiil was not really justified because the 
calculations were done for thick but weak ferromagnets assuming a strong anisotropy of the ferromagnet that was 
n ecessary for a formation of the doma in walls with the magnetization vector changing its sign but not its axis. 

iBergeret. Efetov. and Larkml 1)200(1) investigated theoretically the poss ibility of a cr ypt oferromagnetic-li ke (CF) 
state in S/F bilayers with parameters corresponding to the experiments bv lMuhge et al\ and lOarifullin et all In that 
work a CF state with a magnetic moment that rotates in space was considered. This co rresponds to a weak anisotropy 
of the ferromagnet, which was the case in the samples studied in lMuhge et al\ ()l998j) . In particular. IBergeret et al\ 
(2000) studied a phase transition between the CF and the ferromagnetic (FM) phases. The calculations were carried 
out in the limit dp <C 6i = v o/h, T c <C h <C e , v and e ar e the Fermi v e locity and Fermi energy, resp ectively. This 
limit is consistent with the parameters of the experiment of lMiihge et al\ II1998D . iGar ifullin et al. (2002(). We present 
here the main ideas of this work. 
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FIG. 22 A S/F bilayer consisting of a thin ferromagnet attached to a bulk superconductor. The ferromagnet may be either in 
the (a) ferromagnetic or the (b) cryptoferromagnetic phase. 



The Hamiltonian describing the bilayer structure in FIG. l2*2*l can be written as 

H ( 7 ) = Ho + Hbcs - 7 / dr*+(r) [h(r)<r] Q/3 ^(r) + H M , (5.3) 

where the integration must be taken in the region — d < x < 0. Here Hq is the one-particle electron energy (including 
an interaction with impurities), Hbcs is the usual term describing the conventional BCS superconductivity in the 
superconductor S and the third term describes the interaction between localized moments and conduction electrons, 
where 7 is a constant that will be put to 1 at the end (see Appendix |A"|) . 

The term Hm describes the interaction between the localized moments in the ferromagnet (c/. Ea. H5.2fl ). We assume 
that the magnetization of the localized spins is described by classical vectors and take into account the interaction 
between neighboring spins only. In the limit of slow variations of the magnetic moment in space with account of Eq. 
(15.21) . the Hamiltonian Hm can be written in the form 



H 



M 



J 



(vs B ) 2 + {vs y y + (ys z y 



dV, 



(5.4) 



where the magnetic stiffness J characterizes the strength of the coupling between the localized moments in the F 
layer and the Si are the components of a unit vector that are parallel to the local direction of the magnetization. 

We assume that the magnetic moments are directed parallel to the S/F interface and write the spin vector S as 
S = (0, — sin (9, cos 9). A perpendicular component of the magnetization would induce strong Meissner currents in the 
superconductor, which would require a greater additional energy. 

The condition for an extremum of the energy Hm, Eq. I|5.4|) can be written as 



A9 = 



(5.5) 



Solutions of Eq. I|5.5|l can be written in the form 9 = Qy, where Q is the wave vector characterizing the rotation in 
space (see FIG. I22H . The value Q = corresponds to the ferromagnetic state. 
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What we want to do now is to compare the energies of the ferromagnetic and cryptoferromagnetic states. The 
latter will be considered for the case with a rotating in space magnetic moment 9 — Qy. This should be energetically 
more favorable than the domain-like structure one provided the ma gnetic anisotropy of F is low. Such a CF state 
corresponds to a so called Neel wall (see for example I AharonJ l)l996|) ). 

Strictly speaking, one has to take into account also a magnetostatic energy due to a purely magnetic interaction of 
the magnetic moments. However, if the condition 

ik >>d ' (5 ' 6) 

where M 8 is the mag netic moment p er volume, is fulfilled one can neglect its contribution with respect to the one of 
the exchange energy ljAharonl[l99q) . 

Taking typical values of the parameters for Fe: M s — 800emu/cm 3 and J = 2.10 _6 erg/cm one can see that 
the condition (|5.6|l requires that the thickness d of the ferromagnet is smaller than lOnm, which corresponds to 
comparatively thick layers. Throughout this section this condition is assumed to be fulfilled. 
In this case the magnetic energy Qm (per unit surface area) is given by the simple expression 

9, M = JdQ 2 . (5.7) 

In order to calculate the superconducting energy £1$ one has to take into account the fact that the order parameter 
should be destroyed, at least partially, near the contact with the ferromagnet. This means that the order parameter 
A is a function of the coordinate x perpendicular to the interface. As we want to minimize the energy we should look 
for a non-homogeneous solution for A (x) of non-linear equations describing the superconductivity. Near the critical 
temperature T c one can use Ginzburg-Landau equations. The proper solution of these equations can be written in 
the form 



iW = A " tant l7ifcrm +c J (5 - 8) 

where Ao the value of the order parameter in the bulk, and £g l is the correlation length of the superconductor defined 
in Eq. (|2.2|) . Near T c this length can be much larger than the length £g. The parameter C in Eq. 1)5. 8fl is a number 
that has to be found from boundary conditions. 

The solution for A (x) , Eq. (|5.8I) is applicable at distances exceeding the length and therefore we cannot use it 
near the interface. 

Having fixed the constant C one can compute the decrease of the superconducting energ y due to the suppr ession 
of superconductivity in the S layer using the Ginzburg-Landau free energy functional Ce.g. lde GennesH T966). The 
decrease of the superconducting energy f2g per unit area at the F/S interface is a function of C and can be written 
as 

n s ^^=\r\^(2 + K)(l-K)\ (5.9) 

where K = tanhC, and r = (T - T c )/T c . 

It remains only to determine the contribution from the third term of the Hamiltonian l|5.3|) . The corresponding 
free energy flM/s is given by the expression: 



Vm/S = -MrTVoy^ J q d ^ J d3r ( hCT )(5>o . 



(5.10) 



where vq is the density of states and (g) is the quasiclassical Green's function averaged over all directions of the 
Fermi velocity. 

Since the exchange field h in a strong ferromagnet may be much higher than the value of r _1 (here t is the 
momentum relaxation time), one has to solve the Eilenberger equation in the F region and the Usadel equation in the 
S regi on. Solutions for these equations in both the superconductor and ferromagnet were obtained bv lBergeret et ail 
(2003). _ 

Thus, the total energy is given by f2 = Qm + + ^m/s> Eqs. (|5.71 15. 9| I5.10fl . As a result, one can express the 
free energy as a function of two unknown parameters, K and Q. One can find these parameters from the condition 
that the free energy must be minimal, which leads to the equations 



d tt/dK = dn/dQ = 



(5.11) 
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One can show that the CF-F transition is of second order, which means that near the transition the parameter 
Q is small. At the transition it vanishes and this gives an equation binding the parameters. Solving the equation 
numerically we come to the phase diagram of FIG. 1231 determining the boundary between the ferromagnetic and 
cryptoferromagnetic states. The parameters a and A used in FIG. E31 are defined as 

2^ _ Jd F 7C(3) 
fl =^V' X = VF y/ZT c D* W ' (5 - 12) 

where r/ is the ratio between the Fermi velocities Vq /v§ . It is clear from Eqs. 1)5. 12[) that the parameter a is related 
to the exchange energy h, while A is the related to the magnetic stiffness J . 

The conclusion that the phase transition between F and CF states should be of the second order was drawn 
negl ecting the magn etostatic interaction. The direct magnetic interaction can change this transition to a first order 
one (iBuzdid - EfflSH) . However, in the limit of Eq. i|5.6|) . this first order transition will be inevitably close to the 
second order one. Such a modification of the type of the phase transition is out of the focus of this review. 

Let us make estim ates for the materials used in the experiments. Performing ferromagnetic resonance measurements, 
Miih ge et al\ l|l998f ) have observed a decrease of the effective magnetization of a Nb/Fe bilayer. The stiffness J for 
materials like Fe and Ni is w 6QK/A. The parameters characterizing Nb can be estimated as follows: T c = 10 K , 
v F = 10 8 cm/s, and Z = 100A. The thickness of the magnetic layer is of order rf=10A, and the exchange field h = 10 A K 
which is proper for iron. 

Assuming that the Fermi velocities and energies of the ferromagnet and superconductor are close to each other, we 
obtain a w 25 and A ~ 6.10 -3 . It is cl ear from F IG. |2*3lthat the cryptoferromagnetic state is hardly possible in the 
Fe/Nb samples used in the experiment IMiihge et al\ l)l998l) . 

However, one can in principle explain the observed, decrease of th e magnetization taking a closer look at the 
structure of the S/F interface. In the samples analyzed bv IMiihge et al\ the interface between the Nb and Fe layers is 
rather rough. So, one can expect that in the magnetic layers there were "islands" with smaller values of J and/or h. A 
reduction of these parameters in the Fe/Nb bilayers is not unrealistic because of the formation of non-magnetic "dead" 
layers that can also affect the parameters of the ferromagnetic layers. If the cryptoferromagnetic state were realized 
only on the islands, the average magnetic moment would be reduced but would remain finite. Such a conclusion 
correlates with what one observes experimentally. One can also imagine islands very weakly connected to the rest of 
the layer, which would lead to smaller en ergies of a non-homoge neous state. 

Let us now consider the experiment bv lGarifullin et all J2002) on Pdo.97Feo.03fV. Due to the low concentration of 
iron, the magnetic stiffness and the exchange field of t he F-layers is much lowe r than the one in the case of a pure iron. 
For this system, one estimates the parameters as (see IGarifullin et al\ (|2002T)^ J ~ QQK/nm , h ~ 100K. Assuming 
again that the Fermi velocities of V and Pd\- x Fe x are close to each other. IGarifullin et al\ i|2002|) obtained for the 
sample with d F = 1.2nm the following values of the parameters a ~ 1.2 and A ~ 1.3.1CP 3 . 

Using these values for a and A one can see from the phase diagram in FIG. I23l that there can be a transition from 
the F to the CF state at |r| ~ 0.2, which corresponds to T ~ 2AK. The decrease of the effective magnetization 
M e ff with decreasing temperature was not observed in samples with larger F thickness d F ■ M e ff was a temperature- 
independent constant for the sample with d F — AAnm and ds = 37.2nm. In the sample with d F — l.lnm and 
ds = iOnm the effective magnetization M e ff decreased by ~ 50% w ith cooling from T w 4K to T s» 1.5-ftf. This fact 
is again in accordance with the predictions of lBergeret et all J2000) . 

The results of this section demonstrate that not only ferromagnets change superconducting properties but also 
superconductivity can affect ferromagnetism. This result is valid, in particular, for strong ferromagnets, although the 
thickness of the ferromagnetic layers must be small in this case. 
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The exchange interaction between the superconducting condensate and the magnetic order parameter reduces the 
energy of the system if the direction of the magnetization vector Mp is not constant in space but oscillates. Provided 
the energy of the anisotropy is small, this interaction leads to the formation of a spiral magnetic structure in the F 
him. 

As we will see in the next section the appearance of the CF-state is not the only effect that leads to a reduction 
of the effective magnetization in S/F structures. We will show that the proximity effect may also lead to a change 
of the absolute value of the magnetic moment Mp in the ferromagnet and to an induced magnetization Ms in the 
superconductor. 



B. Ferromagnetism induced in a superconductor 

In the previous section we have seen that the superconductivity can affect the magnetic ordering changing the 
orientation of magnetic moments in the ferromagnetic him. In this section we want to demonstrate that another 
mechanism for a change of the total magnetization of a S/F system exists. In contrast to the phenomenon discussed 
in the previous Section, the orientation of the magnetic moments in the F film does not change but the magnitude of 
the magnetization both in the F and S films does. 

This change is related to the contribution of free electrons both in the ferromagnet (8Mp) and in the superconductor 
(Ms) to the total magnetization. On one hand, the DOS in the F film is reduced due to the proximity effect and 
therefore SMp is reduced. On the other hand, the Cooper pairs in the S him are polarized in the direction opposite 
to Mp, where Mp is the magnetization of free electrons in the ferromagnet. 

Let us consider first a bulk ferromagnet and derive a relation between the exchange field and the magnetization 
of the free electrons. The exchange held h = JS in the ferromagnet can be due to the localized moments (see Eq. 
I|5.1|) l or due to the free electrons in the case of an itinerant ferromagnet 4 In some ferromagnets both the localized 
and itinerant moments contribute to the magnetization. 

The magnetization of the free electrons is given by 



M = \ m J £ J Wf Tif ^ ^ R - ^ np ' (5 - 13) 

where fj,p is an effective Bohr magneton and np is the Fermi distribution function of the free electrons. The expression 
in front of np in Eq. I|5.13|l determines the DOS that depends on the exchange held h. We assume that the 
magnetization is oriented along the z-axis. 

Using Eq. (|5.13(l one can easily compute the contribution of the free electrons to the magnetization in a bulk 
ferromagnet. In the simplest case of a normal metal with a quadratic energy spectrum we have 



(/ ~ f d 3 p ., // ,, / , 



M F = ^ I P 2 dp [n(C p -h)~ n(C p + h)} , (5.14) 



where £ p = p /2m — ep. At T = the magnetization is given by: 



MF0 =2^k)( P +- p ^ (5 - 15) 



where p± — \j2m[ep ± h) are the Fermi momenta for spin up and spin down electrons. In the quasiclassical limit it 
is assumed that h <C ep, and therefore 

MpQ = HBvh, (5.16) 

where v — ppom/w 2 is the density of states at the Fermi level, and ppo — \J2mep is the Fermi momentum in the 
absence of the exchange field 5 . For the temperature range T -C h we are interested in, one can assume that the 
magnetization of the ferromagnet does not depend on T and is given by Eq. I|5.16[) . 



4 In many papers the exchange "field" h is defined in another way (h = —JS) so that the energy minimum corresponds to orientation of 
the vector < a > antiparallel to the vector h. In this case the magnetic moment m = —fj,g < a > is parallel to h. Both definitions lead 
to the same results. 

5 Actually Eq. l5.lb1 is valid not only in the case of a quadratic spectrum but also in a more general case. 
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FIG. 24 S/F structure and schematic representation of the inverse proximity effect. The dashed curves show the local magne- 
tization. 

Now let us consider a S/F system with a thin F layer (see FIG. I24fl and ask a question: Is the magnetization of 
the itinerant electrons modified by the proximity effect? We assume that the exchange field of the ferromagnet F is 
homogeneous and aligned in the z- direction, which is the simplest situation. 

At first glance, it is difficult to expect anything interesting in this situation and, to the best of our knowledge, such 
a system has not been discussed until recently. 

However, physics of this heterostructure is actually very interesting and is general for any shape of the S and F 
regions. It turns out that the proximity effect reduces the total magnetization of the system and this effect can be 
seen as a certain kind of "spin screening" . 

Before doing explicit calculations we would like to explain the phenomenon in simple words. If the temperature is 
above T c , the total magnetization of the system M to t equals Mopdp, where cLf is the thickness of the F-layer. When 
the temperature is lowered below T c , the S layer becomes superconducting and the Cooper pairs with the size of the 
order of £s = \J Dg /2irT c arise in the superconductor. Due to the proximity effect the Cooper pairs penetrate the 
ferromagnet. In the case of a homogeneous magnetization the Cooper pairs consist, as usual, of electrons with the 
opposite spins, such that the total magnetic moment of a pair is equal to zero. The exchange field is assumed to be 
not too strong, otherwise the pairs would break down. 

It is clear from this simple picture that pairs located entirely in the superconductor cannot contribute to the 
magnetic moment of the superconductor because their magnetic moment is simply zero, which is what one could 
expect. Nevertheless, some pairs are located in space in a more complicated manner: one of the electrons of the pair 
is in the superconductor, while the other moves in the ferromagnet. These are the pairs that create the magnetic 
moment in the superconductor. This follows from the simple fact that the direction along the magnetic moment M in 
the ferromagnet is preferable for the electron located in the ferromagnet (we assume a ferromagnetic type of exchange 
field) and this makes the spin of the other electron of the pair be antiparallcl to M. So, all such pairs with one 
electron in the ferromagnet and one in the superconductor equally contribute to the magnetic moment in the bulk of 
the superconductor. As a result, a ferromagnetic order is created in the superconductor, the direction of the magnetic 
moment in this region being opposite to the direction of the magnetic moment M in the ferromagnet. Moreover, the 
induced magnetic moment penetrates the superconductor over the size of the Cooper pairs £s that can be much larger 
than d.F- 

This means that although the magnetization Ms induced in the superconductor is less than the magnetization in the 
ferromagnet Mpo, the total magnetic moment in the superconductor Ms = J s d 3 rMs(r) may be comparable with the 
magnetic moment of the ferromagnet in the normal state Mpo = MfoVf, where Vf = dp in the case of a flat geometry 
(Mfo is the magnetic moment per unit square) and Vf — Anap/i is the volume of the spherical ferromagnetic grain. 
It turns out that the total magnetic moment of the ferromagnetic region ( film or grain) Mp = \xb v f^Vp due to 
free electrons is compensated at zero temperature by the total magnetic moment M$ induced in the superconductor. 
This statement is valid if the condition 



is fulfilled. If the thickness of the F film (or radius of the F grain) is not small in comparison with the correlation 
length £s, the situation changes: the induced magnetic moment M$ is much smaller than Mpo but a variation of 
the magnetic moment of the ferromagnetic film (or grain) SMp becomes comparable with Mfq. The latter is caused 
by a change in the density of states of the ferromagnet due to the proximity effect. However, the case of a large 
ferromagnet size is less interesting because the exchange field h should be smaller than A (the full screening of Mp 
occurs only if the second condition in Ea. H5.17[) is fulfilled). 

Using similar arguments we can come to a related effect: the magnetic moment in the ferromagnet should be 
reduced in the presence of the superconductivity because some of the electrons located in the ferromagnet condensate 
into Cooper pairs and do not contribute to the magnetization. 



A << h « E Th = Dp /dp 



(5.17) 
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^,From this qualitative and somewhat oversimplified picture one can expect that the total magnetization of the S/F 
system will be reduced for temperatures below T c . Both the mechanism studied here and that of the last section 
lead to a negative change of the total magnetization. Thus, independently of th e origin of ferromagnet ism, they can 
explain, at least qualitatively, the experimental data of lMiihge etah §998) and IGarifullin et all (|2002f> . 

The ideas presented above can be confirmed by calculations based on the Usadel equation. In order to determine 
the change of the magnetization it is enough to compute the quasiclassical Green's functions g R ( A ^ = (i/ir) J d£G R ( A ^ 
and, in particular, the component proportional to 73173. 

The matrix Green's function has the form (we write g in Matsubara representation: g(uj) — g R (iui) for positive u>) 

9 = T39 + if 2 f ■ (5.18) 
In the ferromagnet we represent, for convenience, the matrix / in the spin-space as 

f S I ) < 5 - 19 > 

The diagonal form of the matrix is a consequence of the uniformity of the exchange field h. The matrix g has the 
same form. 

In order to find the function 33 that determines the magnetization, we have to solve the Usadel equation i|A18(l in 
the F and S region and to match the corresponding solutions with the help of the boundary conditions (|A21(1 . 

The simplest case when the Usadel equation can be solved analytically is the case of a thin F layer. We suppose 
that the thickness c?f of the F layer is small compared with the charact eristic length £f o f the condensate penetration 
into the ferromagnet (this condition is fulfilled in the experiments bv IGari fullin et al. (2002)). In this case we can 
average the exact Usadel equation <|A18|) over x in the F layer assuming that the Green's functions are almost constant 
in space. In addition, provided the ratio ap/as is small enough, the Green's functions in the superconductor are 
close to the bulk values Jbcs and gscs- This allows us to linearize the Usadel equation in the superconductor. The 
component of the Green's function in S that enters the expression for the magnetization can be obtained from the 
boundary condition (|A21|) and is given by 

9S3(x) = — (-gBcsfFo + fBCS9F3) e KsX , (5.20) 

is us 

where k| = 2y/ui 2 + A 2 /D$, /fo = (/+ + /-)/2, gF3 = {g+ — <7-)/2 and g± and f± are the components of the 
matrices g and /. They are defined as 

g F ± = Q±/Cu±, !f± = ±£bFfBcs/(uj± , (5.21) 



where u>± = uj + f.bF9BCS T it 1 , Cj± = y^± — ( e bF J bc 's) 2 £bF — -Df/(27f^f)- The magnetization variation is 
determined by the expression 

00 

SM = — iirvT Tr (9 ■ o-s) , (5.22) 

UJ — — 00 

Using Eqs. I|5. 2015.22(1 for Tr(<? • a 3 )/2 = g 3 = (g + - g~)/2, one can easily calculate SM. In FIG. |2S1 we show the 
change of the magnetization SM induced in the superconductor as a function of the temperature. We see that for 
low enough temperatures the decrease of t he magnetization can be very large. At the same time, the change of the 
magnetization in the ferromagnet is small ijBergeret et all . 120044 - 

It is interesting to calculate the total magnetic moment SM$ induced in the superconducting film and compare it 
with the total magnetization of the ferromagnet Mp dp (as we have mentioned, the magnetization variation SMp in 
the ferromagnet is small and can be neglected). 

The total magnetization of the superconductor is given by 

SM S = f dxSM s (x) . 

J-d s 

Assuming that h -C e bF — D F /(2^ F dF) or h <C [D F / (2d 2 F )](pFdF / Rb) , we can easily compute the ratio 

SM S D s v s ^ 2 T 



-7T 



Mpodp dFlsVFtbF ^ {oj 2 + A 2 ) 3 / 



S („,2 I A2W2 ~ 1 ' (5 ' 23) 
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FIG. 25 Change of the magnetization of a F/S bilayer as a function of the temperature. 

where pf is the resistivity of the F region. 

We see that in the case of a thin ferromagnet at low temperatures and a not too strong exchange field the mag- 
netization induced in the superconductor compensates completely the magnetization in the ferromagnet. This result 
follows from the fact that the magnetization induced in the superconductor (it is proportional to gss) spreads over 
distances of the order of £g. In view of this result one can expect that the magnetic moment of a small ferromagnetic 
particle embedded in a superconductor should be completely screened by the Cooper pairs. We discuss the screening 
of a ferromagnet particle by the Cooper pairs in the next subsection. 

It is worth mentioning that the problem of finding the magnetization in a S/F structure consisting of thin S 
{ds < £s) and F (dp < £f) layers is equivalent to the problem of magnetic superconductors where ferromagnetic 
(exchange) interaction and superconducting correlations coexist. If we assume a strong coupling between the thin S 
and F layers, we can again average the equations over the thickness of the structure and arrive at the Usadel equation 
for the averaged Green's function with an effective exchange field h — hdp/d and an effective order parameter 

A = Adg/d, where d — ds + dp- In this case the magnetization is given by M = gfiB^Vh 2 — A 2 <d(h — A), where 
(x) is the step function. This means that the total magnetization M is zero for h < A. This result agrees with those 
obtained bv lKarchev et al\ l)200l|) : lshen et all l)2003|) who studied the problem of the coexistence of superconductivity 
and itinerant ferromagnetism in magnetic superconductors. 

One of the assumptions made for obtaining the previous results is the quasiclassical condi tion h/eF <C 1. For 
some m aterials the latter is not fulfilled and one has to go beyond the quasiclassical approach. lHalterman and Vallsl 
(2002a) studied the imbalance of spin up and spin down electrons in pure S/F structures (i.e. without impurities) 
in the case of strong exchange fields (h/ep < 1). In that case superconductivity is strongly suppressed at the S/F 
interface. Solving the Bogoliubov-de Gennes equations numerically the authors showed that there was a magnetic 
"leakage" from the ferromagnet into the superconductor, which lead to a polarization of the electrons in S over the 
short length scale Xp. The direction of the induced magnetic moment in the superconductor was parallel to that in 
the ferromagnet, which contrasts our finding. 

At the same time, the limit of a very strong exchange field considered by lHalterman and Vails! l)2002aj) differs 
completely from ours. It is clear that due to the strong suppression of the superconductivity at the S/F interface, 
the magnetic moment cannot be influenced by the superconductivity and therefore thick ferromagnetic layers with 
exchange energies of the order of the Fermi energy are not suitable for observing the reduction of the magnetization 
described above. 

The DOS for sta tes with spin-up and spin-do wn electrons in a S/F structure has been calculated on the basis of the 
Usadel equation bv lFazio and Lucheronil l)l999|) . The authors have found that the DOS of these states was different in 
the superconductor over the length of the order £g. However, the change of the magnetization has not been calculated 
in this work. 

I has been done later by iKrivoruchko and Koshinal l)2002(l for a S/F structure. Using the Usadel equation, 
the authors numerically calculated the magnetization induced in the superconductor. They found that the magnetic 
moment leaked from the F layer into the S layer and changed the sign at some distance of the order of £g, thus 
becoming negative at sufficiently large distances only. In our opinion, the "leakage" of the magnetic moment Ms 
obtained in that paper is a consequence of the use by the authors of a wrong expression for the magnetic moment. 
They did not add to the formula obtained in the quasiclassical approximation a contribution from the energies levels 
located far from the Fermi energy. The latter contribution is not captured by the quasiclassical approach and should 
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be written additionally. 

We have seen that under certain conditions a finite magnetic moment is induced inside the superconductor. Does this 
magnetic moment affect the superconductivity? The magnetic field B$ in the superconductor equals the magnetization 
AttM s . The induced magnetization in the superconductor Ms is smaller than the magnetization in the ferromagnet: 
Ms — Mp man(dp /{£s, ds})- The critical field for superconducting thin films is given by the expression H c ~ 
(Al/ 1 ds)Hbuik 1 where Xl is the London penetration depth, and Hbuik is the critical field of the bulk material. The 
superconductivity is not affected by the induced field Bs if the field Bs ~ AitMp(dp / Xs) (we set ds ~ £s) is 
smaller than H c . Therefore the condition AnMp < dp)Hb u ik should be satisfied. If we take Xl « 1/j.m and 
dp « 50A, we arrive at the condition AirMp < 200Hi, u ik- This condition is fulfilled easily for the case of not too 
strong ferromagnets. Due to the presence of the magnetization in the ferromagnet and superconductor spontaneous 
currents ar ise in the system. The s pontaneous Me i ssner currents induced by the magnetization in S/F structures were 
studied bv lBergeret et all l)2001c|) : iKrawiec et all ll2004|) . 

The phenomenon discussed in this section can be c onsid ered as an alternative mechanism of the decrease of the 
total magnetic moment observed by iGarifullin et all (|2002f) . In order to clarify which of these two effects is more 
important for the experimental observations one needs more information. 

The most direct check for the cryp toferromagnetic phase would be measurements with polarized neutrons. In 
a recent work by (St ahn et all 12005^) . in which a multilayered S/F/S/F... structure was studied. This structure 
consists of the high T c superconductor YBa2Cu30r (S layer) and of the ferromagnet La^/^Cax/^MnOz (F layer). 
Two samples with the S and F layers of the same thickness were used. Layers of sample 1(2) are 98A(160A) 
thick. The Curie temperature of the ferromagnet and the temperature of the superconducting transition are equal to 
165!^ and 75K respectively. By using neutron reflectometry the authors obtained an information about the spatial 
distribution of the magnetic moment in the structure. Analyzing the temperature dependence of the Bragg peaks 
intensity they came to the conclusion that the most probable scenario to explain important features of this dependence 
observed was the assumption that an induced magnetization arises in the S layers. If this explanation was correct, 
the sign of the induced magnetization had to be opposite to the sign of the magnetization in the F layers. It is quite 
reasonable to think that the mechanism discussed above for conventional superconductors should be present also in 
high T c superconductors and then the theoretic scenario analyzed in this section can serve as an explanation of the 
experiment. 

C. Spin screening of the magnetic moment of a ferromagnetic particle in a superconductor 

Let us consider now a ferromagnetic particle (grain) embedded into a superconductor (see FIG. I26fl . As in the 
previous subsection, we analyze the magnetic moment induced in the superconductor around the particle and compare 
it with the magnetic moment of the F particle (Ana 3 /3)Mp o (we assume that the particle has a spherical form and 
radius a). 




FIG. 26 Ferromagnetic grain embedded into a superconductor. Due to the inverse proximity effect the magnetic moment of 
the grain is screened by the electrons of the superconductor. 

It is well known that the superconducting currents (Meissner currents) in a superconductor screen a magnetic field 
that decays from the surface over the London penetration length Al and vanishes in the bulk of the superconductor. 
The same length characterizes the decay of the magnetic field created by a ferromagnetic (F) grain embedded in a 
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superconductor if the radius of the grain a is larger than A^. However, if the radius a is small, the Meissner effect 
can be neglected and a stray magnetic field around the grain should decay, as in a normal metal, over a length of the 
order a. We consider now just this case. 

Above the critical temperature T c the stray magnetic field polarizes the spins of free electrons and induces a magnetic 
moment. This magnetic moment is very small because the Pauli paramagnetism is weak (fi%u ~ 10 -6 ). In addition, 
the total magnetic moment induced by the stray magnetic field is zero. The penetration depth Al can be of the order 
of hundreds of interatomic distances or larger, so that if a is smaller or of the order of lOnm, the Meissner effect can 
be neglected. 

The screening of the magnetic moment is a phenomenon specific for superconductors. It is usually believed that in 
a situation, when the screening due to the orbital electron motion can be neglected (small grains and thin films), the 
total magnetic moment is just the magnetic moment of the ferromagnetic particle and no additional magnetization is 
induced by the electrons of the superconductor. 

This common wisdom is quite natural because in conventional superconductors the total spin of a Cooper pair is 
equal to zero and the polarization of the conduction electrons is even smaller than in the normal metal. Spin-orbit 
interactions may lead to a finite magnetic susceptibility of the superconductor bu t it is positive and smaller anyway 
than the one in the normal state l)Abrikosovlll988|: lAbrikosov and Gor'kovt Il962j) . 

Let us now take a closer look at the results of the last subsection. We have seen that the proximity effect induces in 
the superconductor a magnetic moment with the sign opposite to the one in the ferromagnet. In view of this result it 
is quite natural to expect that the magnetic moment of a small ferromagnetic particle embedded in a superconductor 
may be screened by the Cooper pairs as it is sketched in FIG. [213 So, let us consider this situation in more detail. 

We consider a ferromagnetic grain of radius a embedded in a bulk superconductor. If the size of the particle 
is smaller than the length £p we can again assume that the quasiclassical Green's functions in the F region are 
almost constant and given by Eq. (|5.21|) . where now CbF = %Dp / (2~fpa). In the superconductor we have to solve the 
linearized Usadel equation for the component gs3 determining the magnetization 

V 2 5 S3 - 4ffS3 - , (5.24) 

where V 2 = d rr + (2/r)d r is the Laplace operator in spherical coordinates. 

Using the boundary conditions Eq. i|A21|) we write the solution of this equation as 

IbCS g — k s0 — a ) 

9S3 = (gBCsfFO - fBCS9F3) -t-. , (5.25) 

7S 1 + nsa r 

where f F0 = (f F+ + /f-)/2 and g F 3 = {qf+ ~ 9f-)/2. 

We assume again that the transmission coefficient through the S/ F interface is not small and the condition A < < 
h < (D F /a 2 ) is fulfilled. In this case the expression for gs3 drastically simplifies. Indeed, in this limit 17^3 = 
fFofBCs/gBCS and f F o = ihfBesgBes/tbF- Therefore Eq. 15.25(1 acquires the form 

IS r e bF 

This solution can be obtained from Eg. Q5. 24(1 if one writes down the term AnA5(r) on the right-hand side of this 
equation with A = fg ( yga 2 ih/('js e bF)- This means that the ferromagnetic grain acts on Cooper pairs as a magnetic 
impurity embedded into a dirty superconductor. It induces a ferromagnetic cloud of the size of the order £g with a 
magnetic moment ~ —fj,B^hV F . 

In order to justify the assumptions made above we estimate the energy Dp /a 2 assuming that the mean free path 
is of the order of a. For a — 30A and vf — 10 8 cm/ sec we get Df/o 2 ~ lOOOif;. This condition is fulfilled for 
ferromagnets with the exchange energy of the order of several hundreds K. 

In the limit of low temperatures the calculation of the magnetic moment becomes very easy and we obtain for the 
magnetic moment Ms induced in the superconductor the following expression 

Ms = 1 (5.27) 



M F0 {^a 3 /3) 



This is a remarkable result which shows that the induced magnetic moment is opposite in sign to the moment of 
the ferromagnetic particle and their absolute values are equal to each other. In other w ords, the magnetic moment of 
the ferromagnet is completely screened by the superconductor l)Bergeret et aLlf2004bf) . The characteristic radius of 
the screening is the coherence length £g, which contrasts the orbital screening due to the Meissner effect characterized 
by the London penetration depth Al- 
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To avoid misunderstanding we emphasize once again that the full screening occurs only if the magnetization (per 
unit volume) of the ferromagnetic grain Mp is given by Eq. l|5.16|) . which means that the ferromagnetic grain is 
an itinerant ferromagnet. If the magnetization of the ferromagnet is caused by both localized moments (Mi oc ) and 
itinerant electrons (Main), the full screening is not achieved. Moreover, the magnetization M; oc may be larger than 
Main and have opposite direction. In this case we would have an anti-screening l|Bergeret and Garciall2004l) . 

Actually, we have discussed the diffusive case only. However, it turns out that the spin screening occurs also in the 
clean case provided the exchange field is not too high : h << vp/dp, where vp and dp are the Ferm i velocity and the 
thickness (radius) of the ferromagnetic film or grain ijBergeret et all l2005t iKharitonov et alj l2005|) . 

The energy spect r um of a superc ondu ctor with a poi nt-like classical magnet ic moment was stud ied many years 
ago bv lShibal l)l968|) . ISakural l)l970|) and lRusinovl l)l969l) . and more recently bv lSalkola et all l)l997j) . The magnetic 
impurity leads to a bound state (3q inside the superconducting energy gap. There is some critical strength h c ~ ep 
of the exchange coupling h that separates two different ground states of the system denoted by ip if h < h c and tp' if 
h > h c . The bound state /3q corresponds to a localized quasiparticle with spin "up" 6 . Since the total electronic spin 
in the state ip is zero one says that the continuum localizes a spin "up" . The energy needed to create a quasiparticle 
excitation decreases when increasing h. At h = h c the state tp becomes unstable against a spontaneous creation of an 
excitation with spin "up" and the transition to the state tp' occurs. In this state the electronic spin at the impurity 
site is now equal to —1/2. All the works considering this problem focused the attention on the subgap structure of the 
spectrum and did not addressed the problem of the screening of the magnetic moment by the continuum spectrum. 
This is of no surprise because a sufficiently large magnetic moment of the impurity (S 3> 1) cannot be screened by 
the quasiparticles. 

D. Spin-orbit interaction and its effect on the proximity effect 

In this section we discuss the influence of the spin-orbit (SO) interaction on the proximity effect. Although in 
general its characteristic energy scale is much smaller than the exchange energy h, it can be comparable with the 
superconducting gap A and therefore this effect can be very important. Since the SO scattering leads to a mixing of 
the spin channels, we expect that it will affect not only the singlet component of the condensate but also the triplet 
one in the ferromagnet. 

In conventional superconductors the SO interaction does not affect thermodynamic properties. However, a non- 
vanishing magnetic susceptibility at zero temperature (Knight shift) obs erved in small superconducting samples and 
films can be explained only if the SO interaction is taken into account ijAbrikosov and Got 'kovl 119621) . In the F/S 
structures considered here the exchange field h breaks the time-reversal symmetry in analogy to the external magnetic 
field in the Knight shift problem. Therefore the SO interaction in the superconductor is expected to influence the 
inverse proximity effect studied in this Chapter. 

In this Section we will generalize the analysis of the long-range proximity effect and the inverse proximity effect 
presented above taking th e SO interaction i nto a ccount. The quasiclassical equations in the prese nce of the SO 
intera ction were derived bv lAlexander et "all l|l985|) and used for the first time for the F/S systems bv lDemler et all 
11199711 . 

The derivation of these equations is presented in the Appendix^ The resulting Usadel equation takes the form 

- iDO r (gd r g) + i (f 3 d t g + d v gf 3 ) + [A,g] + [hS,g] 

+ — [St 3 gt 3 S,g] = 0. (5.28) 

T~s.o. 

All symbols are defined in the Appendix The spin-orbit relaxation time r s . . takes very different values dependin g 
on the material used in the experiments. Some estimates for the values of l//ir s o .can be found in lOh et all fcOOOft . 
For example, for transition metals like Fe one obtains l/hr s , . ~ 10 -2 , while for a typical magnetic rare earth the 
value l//ir s o . ~ 0.3 is more typical. In the latter case the SO interaction should clearly affect the penetration of the 
condensate into the ferromagnet. 

In order to study the influence of the SO interaction on both the long-range and the inverse proximity effect we will 
use Ea. (|5.28|l . We consider first the well known problem of the Knight shift. This example will show the convenience 
of using the quasiclassical approach. 



One assumes that the magnetic impurity has spin up. 
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The Knight shift in superconductors 

Since the pioneering work of lAbrikosov and Gor it is well established that the magnetic susceptibility of 

small superconducting samples is no t zero due to the spin - orbit interaction. This explains the experiments performed 
for the first time many years ago bv lAndroes and Knight! l)l96lh who used the nuclear magnetic resonance technique. 

Let us consider a superconductor in an external magnetic field H . In the Usadel equation, Eq. I|5.28|l . the field H 
plays the role of the exchange energy h. We are interested in the linear response to this filed, i.e. in the magnetic 
susceptibility xs of the superconductor. We assume that the superconductor is homogeneous and therefore we drop 
the gradient term in Eq. H5.28|) : 

-u[f 3 ,g]+i[A,g] + iH[n,g]-(l/r s . .) [Sf 3 gf 3 S,g] = 0, (5.29) 

g 2 = 1 . (5.30) 

The solution of Eq. (|5.29|1 has the form 

g = (gBCS + 930-3) h + (fBCsa-3 + fo) , (5.31) 

where the functions 173 and fo are corrections to the normal gscs an d anomalous fscs Green's functions. In the 
particle- hole space the matrix g has the usual form, i.e. it is expanded in matrices f 3 and if 2- In the spin space the 
triplet component (the g 3 and fo terms) appears due to the magnetic field acting on the spins. Using Eqs. 1)5.2915.31(1 
one can readily obtain 

93 = -'m^rh^j ■ (5 - 32) 



where E u = VA 2 +uj 2 . 

Substituting Eq. I|5.32|l into Eq. 1)5.22)1 we can write the magnetization M as follows 

« a M,- M »( m ^ ff j 4/T| j il (5.33) 

The first term in Eq. (|5.33|) cannot be calculated in the framework of the quasiclassical theory and one should use 
exact Green's functions. It corresponds to the Pauli paramagnetic term given by M — figvH. In the quasiclassical 
approach this term is absent. This term does not depend on temperature on the energy scale of the order of T c and 
originates from a contribution of short distances where the quasiclassical approximation fails. 

This situation is rather typical for the quasiclassical approach and one usually adds to fo rmulae obtained within 
this approach contributions coming from short dis tances or times by hand (see, for e xample , ^^mgnko^md YQlkovl 
(1980l): lKopninl l)200l[k iRammer and Smith! (|l986h ). Eq. lj533j) was first obtained bv lAbrikosov and Gor'kovl l)l962f) 

In the absence of the spin orbit interaction the magnetization at T = is, as expected, equal to zero. However, if 
the SO interaction is finite the spin susceptibility xs does not vanish at T = 0. It is interesting that, as follows from 
Eq. 1)5. 29|) . the singlet component of the condensate is not affected by the SO interaction. The origin of the finite 
susceptibility is the existence of the triplet component fo of the condensate. 

In the S/F structures there is no exchange field in the superconductor and therefore the situation is in princi- 
ple different. However, we have seen that due to the proximity effect the triplet component fo is induced in the 
superconductor. 

From the above analysis one expects that the SO interaction may affect the penetration length of such component 
in the superconductor. In the next sections we consider the influence of the SO the superconducting condensate in 
both the ferromagnet and the superconductor. 



1. Influence of the Spin-Orbit interaction on the long-range Proximity Effect 

Now we consider again the S/F/S/F/S structure of Section llV. Al and assume that the long-range triplet component 
is created, which is possible provided the angle a between the magnetizations differs from and n. In order to 
understand how the SO interaction affects the triplet component it is convenient to linearize Eq. I|5.28|l in the F-layer 
assuming, for example, that the proximity effect is weak. One can easily obtain a linearized equation similar to Eq. 
(|3.15|l for the condensate function /. The solution of this equation is represented again in the form 



f(x) = if 2 ® {fo(x)a + fs(x)a 3 ) + if\ <8> fi{x)ai 



(5.34) 
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The functions fi(x) are given as before by fi(x) = bj exp[njx] but now the new eigenvalues Kj are written as 

(5.35) 
(5.36) 

We see from these equations that both the singlet and triplet components are affected by the spin-orbit interaction 
making the decay of the condensate in the ferromagnet faster. In the limiting case, when 4/r so > h, T c , both the 
components penetrate over the same distance £ s . . = \Jt so Df and therefore the long-range effect is suppressed. In 
this case the characteristic oscillations of the singlet component are destroyed ((Dernier et all Il997). In the more 
interesting case 4/r so ~ T c < h, the singlet component does not change and penetrates over the short distance £f 
. At the same time, the triplet component is more sensitive to the spin-orbit interaction and the penetration length 
equals min(£ so ,£ T )> £ F . 

Therefore, if the spin-orbit interaction is not very strong, the penetration of triplet condensate over the long 
distances discussed in the preceding sections is still possible, although the penetration length is reduced. 



2. Spin-Orbit Interaction and the Inverse Proximity Effect 

Studying a S/F bilayer we have seen that the induced magnetic moment in the superconductor S is related to 
the appearance of the triplet component fa. Moreover, we have shown that this component is affected by the SO 
interaction, while the singlet one f% is not. So, one should expect that the SO interaction may change the scale over 
which the magnetic moment is induced in the superconductor and one can estimate easily this length. 

Assuming that the Green's functions in the superconductor take values close to the bulk values we linearize the 
Usadel equation l|5.28|l in the superconductor. The solution has the same form as before, Eq. (|5.2Ufl . but ks should 
be replaced by 

4 - 4 + 4 D , (5.37) 

where n 2 so — &Ds/t so . Therefore, the length of the penetration of gs3 and, in its turn, of Ms into the S region 
decreases if 4 ~ Cg 2 < 4<r 

In principle, one ca n measure the spatial distribution of the magnetic moment in the S region as it was done by 
iLuetkens etal\ l)2003(l by means of muon spin rotation and get an information about the SO interaction in supercon- 
ductors. As Eq. H5.37fl shows, this would be an alternative method to me asure the strength of the SO interaction in 
superconductors, complementary to the measurement of the Knight shift ijAndroes and Knightl Il96l|l . 



VI. DISCUSSION OF THE RESULTS AND OUTLOOK 

In this review we have discussed new unusual properties of structures consisting of conventional superconductors 
in a contact with ferromagnets. It has been known that such systems might exhibit very interesting properties like a 
non-monotonous reduction of the superconducting temperature as a function of the thickness of the superconductor, 
possibility of a 7r-contact in Josephson junctions with ferromagnetic layers, etc. 

However, as we have seen, everything is even more interesting and some spectacular phenomena are possible that 
even might look at first glance as a paradox. The common feature of the effects discussed in this review is that almost 
all of them originate in situations when the exchange field is not homogeneous. As a consequence of the inhomogeneity, 
the spin structure of the superconducting condensate function becomes very non-trivial and, in particular, the triplet 
components are generated. In the presence of the inhomogeneous exchange field, the total spin of a Cooper pair is 
not necessarily equal to zero and the total spin equal to unity with all projections onto the direction of the exchange 
field is possible. 

We have discussed the main properties of the odd triplet superconductivity in the S/F structures. This super- 
conductivity differs from the well known types of superconductivity: a) singlet superconductivity with the s-wave 
(conventional T c superconductors) and d-wave (high T c superconductors) types of pairing; b) odd in momentum p 
and even in frequency u> triplet superconductivity observed, e.g., in SV2-RWO4. 

The odd triplet superconductivity discussed in this Review has a condensate (Gor'kov) function that is an odd 
function of the Matsubara frequency u> and an even function (in the main approximation) of the momentum p in the 
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diffusive limit. It is insensitive to the scattering on nonmagnetic impurities and therefore may be realized in thin film 
S/F structures where the mean free path is very short. 

For the first time, the condensate function of this type has been suggested bv lBerezinskiil (|l975fl many years ago as 
a possible candidate to describe superfluidity in He 3 . Later, it has been established that the superfluid condensate 
in He 3 had a different structure - it was odd in p and even in u>. In principle, there is an important difference between 
the triplet superconductivity discussed here and that predicted by Berezinskii who assumed that the order parameter 
A was also an odd function of lu. In our case the order parameter A is determined by the singlet, s-wave condensate 
function and has the ordinary BCS structure (i.e., it does not depend on the momentum p and frequency to). On the 
other hand the structure of the triplet condensate function / in the diffusive case considered here coincides with 
that suggested by Berezinskii: it is an odd function of the Matsubara frequency u) and, in the main approximation, is 
constant in the momentum space. The antisymmetric part of / is small compared with the symmetric part, being 
odd in p and even in w. 

The triplet component with the projection of the total spin S z = ±1 penetrates the ferromagnet over a long 
distance of the order of £jv ~ \ZDf/2ttT, which shows that the exchange field does not affect the triplet part of the 
condensate. At the same time, the exchange field suppresses the amplitude of the singlet component at the S/F 
interface that determines the amplitude of the triplet component. The long-range triplet component arises only in 
the case of a nonhomogeneous magnetization. The triplet component appears also in a system with a homogeneous 
magnetization but in this case it corresponds to the projection S z — and penetrates the ferromagnet over a short 
length £p = ^D F /h «£ N . . 

The triplet component exists also in magnetic superconductors l)Bulaevskii et adll985llKuric and Kulidl200lj) with 
a spiral magnetic structure. However, it always coexists with the singlet component and cannot be separated from 
it. In contrast, in the multilayered S/F structures with a nonhomogeneous magnetization and with the thickness of 
the F layers cLf exceeding £p, the Josephson coupling between S layers is realized only through the long-range triplet 
component and this separates the singlet and triplet components from each other. As a result, the "real" odd triplet 
superconductivity may be realized in the transverse direction in such structures. 

Another interesting peculiarity of the S/F structures is the inverse proximity effect, namely, the penetration of the 
magnetic order parameter (spontaneous magnetic moment M) into the superconductor and a spatial variation of the 
magnetization direction in the ferromagnet under the influence of the superconductivity. It turns out that both effects 
are possible. A homogeneous distribution of the magnetization Mp in the S/F bilayer structures may be energetically 
unfavorable in F even in a one-domain case resulting in a nonhomogeneous distribution of in the ferromagnet. 

Moreover, the magnetic moment penetrates the superconductor (induced ferromagnetism) changing sign at the S/F 
interface. Therefore the total magnetic moment of the system is reduced. Under some condition the full spin screening 
of Mp occurs. For example, at zero temperature the itinerant magnetic moment of a ferromagnetic grain embedded 
into a superconductor is completely screened by spins of the Cooper pairs in S. The radius of the screening cloud is 
of the order of the superconducting coherence length £5. If the magnetization vector Mp is oriented in the opposite 
direction to the ferromagnetic exchange field h, the anti-screening is possible. 

As concerns the experimental situation, certainly there are indications in favor of the long-range triplet component, 
although an unambiguous evidence does not exist so far. For example, the resistance of ferromagnetic films or 
wir es in the S/F structures changes on distances that exceed the length of the dec ay of the singlet component 
£h l)Aumentado and Chandrasekharl l200lt iGiroud et all Il998t iPetrashov et all Il999l) . A possible reason for this 
long-range proximity effect in the S/F systems is the long-range penetration of the triplet component. However a 
simpler effect might also be the reason for this long-range proximity effect. It is related to a rearrangement of a 
domain structure in the ferromagnet when the temperature lowers below T c . The Meissner currents induced in the 
superco nductor by a stray m agnetic field affect the domain structure, and the resistance of the ferromagnet may 
change l|Dubo nos et al 1 12002^ . At the same time, the Mei ssner currents shoul d be considerably reduced in an one 
dimensional geometry for the ferromagnet like that used in jGiroud et all Il998j) and the explanation in terms of the 
lo ng range pene t ration of the triplet component are more probable here. 

ISefrioui et al\ l)2003j) also obtained some indications on the existence of a triplet component in a multilayered 
S/F/S/F... structure. The samples used bv lSefrioui et al\ contained the high T c material YBa 2 Cu 3 07 (as a super- 
conductor) and the half-metallic ferromagnet LaQ^CaQ.^MnO^ (as a ferromagnet). They found that superconductivity 
persisted even in the case when the thickness of the F layers (If essentially exceeded £f (<^f ^ lOnm and £f ~ 5nm). 
In a half-metal ferromagnet with spins of free electrons aligned in one direction the singlet Cooper pairs cannot exist. 
Therefore it is reasona ble to assume that the superconducting coupling between neighboring S layers is realized via 
the triplet component ijEschrig et all l2003t IVolkov et all\200$) . 

A red uction of the magnetic moment of the S/F structures due to superconducting correlations has been observed 
already l|Garifullin et ad 12002). This re duction may be caused both by the spin screening of the magnetic moment 
Mp and by the rotation of Mp in space llBereeret et all 120001 l2004a|) . Perhaps, the spin screening can be observed 
directly by probing the spatial distribution of the magnetic field (or magnetic moment M) with the aid of the muon 
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spin rotation technique ((Luetkens et all l2003|) . The variation of the magnetic moment M occurs on a macroscopic 
length £5 and therefore can be detected. 

An evidence in favor of the inverse proximity effect has also been obtained in another experimental work 
l)Stahn et aZJ .l2005). Analyzing data of neutron refiectometry on a multilayered Y Ba^Cu^O'j / ' La2/zCaiisMnOz struc- 
ture, the authors concluded that a magnetic moment was induced in the superconducting Y Ba2Cu^,Oj layers. The 
sign of this induced moment was opposite to the sign of the magnetic moment in the ferromagnetic La2/zCai/j,MnOs, 
layers, which correlates with our prediction. 

In spite of these experimental results that may be considered as, at least preliminary, confirmation of the existence 
of the triplet component in the S/F structures, there is a need in additional experimental studies of the unconventional 
superconductivity discussed in this review. One of the important issues would be to understand whether the long 
range proximity effects already observed experimentally are due to the triplet pairing or to a simple redistribution 
of the domain walls by the Meissner currents. We believe that measurements on thin ferromagnetic wires where the 
Meissner currents are reduced may clarify the situation. 

It is very interesting to distinguish between the two possible inverse proximity effects experimentally. Although both 
the formation of the cryptoferromagnetic state and the induction of the magnetic moments in the superconductors 
are very interes t ing ef fe cts, it is not c lear y et which of these effects causes the magnetization reduction observed by 
lOarifullin eial\ ll2002ft : iMiihge e< al\ lll998t) . 

The enhancement of the Josephson current by the presence of the ferromagnet near the junction is one more 
theoretical prediction that has not been observed yet but, certainly, this effect deserves an attention. An overview 
for experimentalists interested in all these subjects is presented in Appendix [5] where we discuss briefly different 
experiments on S/F structures, focusing our attention on the materials for which, we expect, the main effects discussed 
in this review may be observed. 

In addition, further theoretical investigations are needed. The odd triplet component has been studied mainly 
in the diffusive limit (hr << 1). It would be interesting to investigate the properties of the triplet component for 
an arbitrary impurity concentration (hr ^ 1). No theoretical work on dynamics of magnetic moments in the S/F 
structures has been performed yet, although the triplet component may play a very important role in the dynamics 
of the S/F structures. Transport properties of the S/F structures require also further theoretical considerations. It 
would be useful to study the influence of domain structures on properties of the S/F structures, etc. In other words, 
physics of the proximity effects in the superconductor-ferromagnet structures is evolving into a very popular field of 
research, both experimentally and theoretically. 

The study of the proximity effect in S/F structures may be extended to include ferromagnets in contact with high 
temp erature superconductors. Although some experiments have been done already l)Sefrioui et all l2003HStahn et all 
2005), one can expect much more broad experimental investigations in the future. The modern technique allows the 
preparation of multilayered S/F/S/F.. structures consisting of thin ferromagnetic layers (as La 2 /3Cai/ 3 Mn03) and 
thin layers of high T c superconductor (as Y BC12C1L3O7) with variable thicknesses. It would be very interesting to 
study, both experimentally and theoretically, such a system with non-collinear magnetization orientations. In this 
case d-wave singlet and odd triplet superconductivity should coexist in the system. It is known that many properties 
of the ordinary BCS superconductivity remain unchanged in the high T c superconductors. This means that many 
effects considered in this review can also occur in S/F structures containing high T c materials, but there will certainly 
be differences with respect to the conventional superconductors with the s-pairing. 

We hope that this review will encourage experimentalists and theoreticians to make further investigations in this 
fascinating field of research. 
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APPENDIX A: Basic equations 

Throughout this review we use mainly the well established method of quasiclassical Green's functions. Within 
this method the Gor'kov equations can be dras tically simplified by integrating the Green's func t ion ov er the mo- 
mentum. T his method wa s first introduced bv lLark in and Ovchinnikovl l)l968j) and lEilenbergerl l)l968(l . and then 
extended by lUsad el ( 1970) for a dirty case and by lEliashberd 1119711) for a non-e q uilibrium case . The me t hod o f 
the quasiclassical Green's functions is discussed in many reviews feelzig et al\ l)l999j) : lLarkin and Ovchinnikovl l)l984fk 
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iRammer and S^ithUT986> : Is erene and Reined (^983) and in the book bv lKopnir] l)200l|) . In this Appendix we present 
a brief derivation of equations for the quasiclassical Green's functions and write formulae for the main observable quan- 
tities in terms of these functions. A special attention will be paid to the dependence of these functions on the spin 
variables that play a crucial role in the S/F structures. In particular, we take into account the spin-orbit interaction 
alongside with the exchange interaction in the fcrromagnet. 

We start with a general Hamiltonian describing a conventional BCS-superconductor/ferromagnet structure: 

H = Yj i a i [((4><W + eV ) + u *mp) 6ss ' + Us °- 

- (h.a)] a s , p , - (Aa^at, pl + c.c.) } . (Al) 

The summation is carried out over all momenta (p,p') and spins (s, s') (the notation s, p means inversion of both spin 
and momentum), £ p = p 2 /2m — ep is the kinetic energy counted from the Fermi energy tp, V is a smoothly varying 
electric potential. The superconducting order parameter A must be determined self-consistently. It vanishes in the 
ferromagnetic regions. The potential Ui„ lp = U(p — p') describ es the interaction of the elect rons with nonmagnetic 
impurities, and U s . . describes a possible spin-orbit interaction l|Abrikosov and Gor'kovl Il962lk 

(t) 

Us.o. = V (P Xp>. 

i p F 

Here the summation is performed over all impurities. 

The representation of the Hamiltonian in the form (|A1|) implies that we use the mean-field approximation for the 
superconducting (A) and magnetic (h) order parameter. The exchange field h is parallel to the magnetization M.p 
in F 7 . In strong ferromagnets the magnitude of h is much higher than A and corresponds to an effective magnetic 
field H exc — h/nB of the order 10 6 Oe (where [Ib = g^Bohr, 9 is the ^-factor and [iBohr is the Bohr magneton). 

In order to describe the ferromagnetic region we use a simplified model that catches all physics we are interested 
in. Ferromagnetism in metals is caused by the electron-electron interaction between electrons belonging to different 
bands that can correspond to localized and conducting states. Only the latter participate in the proximity effect. If 
the contribution of free electrons strongly dominates (an itinerant ferromagnet), one has Mp = M e and the exchange 
energy is caused mainly by free electrons. 

If the polarization of the conduction electrons is due to the interaction with localized magnetic moments, the 
Hamiltonian Hp can be written in the form 

Hp = -h X ^ { a i S * °ss«ls'p>} (A2) 
{P, s } 

where S = ^2 a S a S(r — r a ), S a is the spin of a particular ion. A constant h\ is related to h via the equation: 
h = HiUmSq , where Um is the concentr ation of magnetic ions and So is a maximum value of S a (we consider 
these spins as classical vectors; see Ref. l|Gor'kov and Rusinovi llQfil l. In this case the magnetization is a sum: 
M = Mice + M e , and the magnetization M e can be aligned parallel (hi > 0, the ferromagnetic type of the exchange 
field) to M or antiparallel (hi < 0, the antiferromagnetic type of the exchange field). In the following we will assume 
a ferromagnetic exchange interaction (M e and M are oriented in the same direction). In principle, one can add to Eq. 
(|A2I) the term X]{a b} * ^b} describing a direct interaction between localized magnetic moments but in the most 
part of the review this term is not important except Section V.A. whe re the cryptoferromagnetic sta te is discussed. 

Starting from the Hamiltonian i|Al|l and using a standard approach l|Larkin and Ovchinnikovlll984|) . one can derive 
the Eilenberger and Usadel equations. Initially these equations have been derived for 2x2 matrix Green's functions 
g n ,n>, where indices n,n' relate to normal (311,522) and anomalous or condensate (/i2,/2i) Green's functions. These 
functions describe the singlet component. In the case of a non- homogenous magnetization considered in this review 
one has to introduce additional Green's functions depending on spins and describe not only the singlet but also the 
triplet component. These matrices depend not only on n,n' indices but also on the spin indices s, s', and are 4x4 
matrices in the spin and Gor'kov space (sometimes the n, n' space is called the Nambu or Nambu-Gor'kov space). 

In order to define the Green's functions in a customary way it is convenient to write the Hamiltonian (|A1I) in terms 
of new operators c^ s and c nsp that are related to the creation and anhilation operators af and a s by the relation 



7 (we remind that the exchange field h is measured in energy units, see also the Footnote on page 78) 
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(we drop the index p related to the momentum) 

2s, n = 1 



t o (A3) 

These operators (for s = 1) were introduced by Nambu l|NambuL Il96(1) . The new operators allow one to express 
the anomalous averages < Of • a-\ > introduced by Gor'kov as the conventional averages < c\ ■ > and therefore 
the theory of superconductivity can be constructed by analogy with a theory of normal systems. Thus, the index n 
operates in the particle-hole (Numbs-Gor'kov) space, while the index s operates in the spin space. In terms of the c ns 
operators the Hamiltonian can be written in the form 

H = Yl C ns^(nn')(ss')Cn's' , (A4) 

where the summation is performed over all momenta, particle-hole and spin indices. The matrix Ti. is given by 

^ = 2 { ^ p6pp ' + eV "> + Umpi ^3 ® CT + A ® ct 3 - hf 3 S 

+ E#(PXP')4 • (A5) 

The matrices Ti and &i are the Pauli matrices in the particle-hole and spin space respectively; i = 0, 1, 2, 3, where fo 
and Co are the corresponding unit matrices. The matrix vector S is defined as 

S = ((71, 172, T3CT3) , 

and the matrix order parameter equals A — fii?eA — T2-ZmA. Now we can define the matrix Green's functions (in 
the particlc-holc<X>spin space) in the Keldysh representation in a standard way 

G{U,t' k ) = 1 (t c (c ns (U)cl s ,(t' k ))) , (A6) 

where the temporal indices take the values 1 and 2, which correspond to the upper and lower branch of the contour 
C, running from —00 to +00 and back to —00. 

One can introduce a matrix in the Keldysh space of the form 

where the retarded (advanced) Green's functions G{t^t') R ^ are related to the matrices G{ti,t' k ) : G(t,t') R ( A ' = 
G(ti,t[) — G(ti(2), £ 2 (i))- AH these elements are 4x4 matrices. These functions determine thermodynamic properties 

of the system (density of states, the Josephson current etc). The matrix G(t, t') = G{t\, t' 2 ) + Gfe, t[) is related to 
the distribution function and has a nontrivial structure only in a nonequilibrium case. In the equilibrium case it is 
equal to: G(e) K = J d{t - t')G{t - t') K exp(ie(t - t')) = [G(e) R - G(e) ] tanh(e/2T). 

In order to obtain the equat ions for the quasiclassical Green's functions, we follow the procedure introduced by 
lLarkin and Ovchinnikovl |)1984|) . The equation of motion for the Green's functions is 

(id t - H - t imp - E s . .) G = 1 , (A8) 

where 

d 2 - ~ 

H = -f 3 -^ - e F - hf 3 S + A ® <7 3 
2m 

and "Simp and are the self-energies given in the Born approximation by 



d£ P / -r-;G 

^s.o. ~~ N imp'U's.o. {^-*) s.o. i 

(G) s .o. = »J d ZpJ Jnxn')SGS(nxn'). 



(A9) 
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Here Ni mp is the impurity concentration, v is the density of states at the Fermi level and n is a unit vector parallel 
to the momentum. 

Next step is to subtract from Ea. (jA8|) . multiplied by f 3 from the left, its conjugate equation multiplied by t 3 from 
the right. Then one has to go from the variables (r, r') to ((r + r')/2, r — r') and to perform a Fourier transformation 
with respect to the relative coordinate. By making use of the fact that the Green's functions are peaked at the Fermi 
surface, one can integrate the resulting equation over £ p , and finally one obtains 

f 3 d t g + d t > gf 3 + v P V<? - i [hS,g] - i [A,g] 

+ ^ [<<?), Si + ^—[r 3 (9)s.o.T 3 ,9] =0 (A10) 

where A = 73A and the quasiclassical Green's functions g{ti, t' k ) are defined as 

#(p P ,r) = ^(f 3 <g><7 ) J d£ p G(ti,t' k ;p,r) , (All) 

and vf is the Fermi velocity. The scattering times appearing in Eq. (IA10|) are defined as 

r- 1 = 2<KvN imp ut mp (A12) 

1 1 f dQ 2 „ 

T s.o. = T^vN lmp J — < G sin 9 (A13) 

Eq. (|A10|) is a generalization of an equation derived bv lLarkin and Ovchinnikovl l|l968l) . and lEilenbergerl {l968) 
for a general noneq uilibrium case. This generalizatio n (in the absence of spin - depen dent interactions) has been done 
bv lEliashberdlli"97i and Larkin and Ovchinnikov l)Larkin and Ovchinnikovl, 1 1984() . A solution for Eq. l|A10|) is not 
unique. The proper solutions must obey the so called normalization condition 

J (dei/27r)^(pF,r;e,ei).p(p F ,r;ei,e') = 1 (A14) 

Generalization for the case of exchange and spin-orbit interaction was presented in ijBereeret et all l200df) and 
(Bcre eret et all feOQlo?) . The solution for Ea. (|A10|l can be obtained in some limiting cases, for example, in a homo- 
geneous case. However finding its solution for nonhomogeneous structures with an arbitrary impurity concentration 
may be a quite difficult task. Further simplifications can be made in the case of a dirty superconductor when the 
energy r _1 related to the elastic scattering by nonmagnetic impurities is larger than all other energies involved in 
the problem, and the mean free path I is smaller than all characteristic lengths (except the Fermi wave length that is 
set in the quasiclassical theory to zero) . In this case one can expand the solution of Eq. (|A10I) in terms of spherical 
harmonics and retain only the first two of them, i.e. 

<7(PF,r;) = g s (r) + (p F /p F )g a (r) , (A15) 

where <7 s (r) is a matrix that depends only on coordinates. The second term is the antisymmetric part (the first 
Legendre polynomial) that determines the current. It is assumed that the second term is smaller than the first one. 
The parameter 1/xq determines it's smallness, where I is the mean free path and Xq is a characteristic length of the 
problem. In S/F structures xq w y/Dp/h is the shortest length because usually ft, > A. In the limit 1/xq « 1, that 
is, if the product hr is small, one can express g a (r) from Ea. (|A10|) in terms of <? s (r) 



g a (r;e,e') = -Z&(r;e, ei)Vg s (r;ei, e'), (A16) 

When obtaining Eq. i|A16|) . we used the relations 

&(r;e,ei) o &(r;ei,e') = 1 , (A17) 
g as (r;e,ei) o g s {r;ex,e') + g s (r;e, ei) o g a (r;ei, e') = 

The symbolically written products in Eas. (|A16|) . i|A17| ) imply an integration over the internal energy t\ as it is shown 
in 

The equation for the isotropic component of the Green's function after averaging over the direction of the Fermi 
velocity vp reads 

- iDV{g\7g) + i {f 3 d t g + ft,$f&) + [A, g] + [hS, g) 

[Sf 3 gf 3 S,g] =0, (A18) 



% 

7~s.o. 
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where D is the diffusion coefficient. 

If we take the elements (If) or (22) of the supermatrix g, we obtain the Usadel equation for the retarded and 
advanced Green's functions g RlyA \t,t') generalized for the case of the exchange field acting on the spins of electrons. 
In this review we are interested mainly in stationary processes, when the matrices g R ^ A \t, t') depend only on the time 
difference (t — t'). Performing the Fourier transformation g R ^ A \e) = J d(t — t')g R ( A \t — t') exp(ze(i — t')), we obtain 
for g R ( A \e) the following equation (we drop the indices R(A)) 



Dd x (gd x g) +ie[T 3 a ,g]+ ih{[%a 3 ,g] cos a(x) + [T a 2 ,g}sma(x)} + i [A, g] H [Sf 3 <7f 3 S, <?] =0. (A19) 

Ts.o. 

It is assumed here that h has the components /i(0, sin a, cos a). This equation was first obtained bv lUsadel il97(l> 
and it is known as the Usadel equation. An inclusion of the exchange and spin-orbit interaction was made in 
l|Alexander et all 119851: iDemler et aZLll997|) . 

Eq. JAT8)| can be solved analytically in many cases and it is used in most of previous sections in order to describe 
different S/F structures. Solutions for the Usadel equation must obey the normalization condition 

3(p F ,r;e).g(p F ,r;e) = 1 (A20) 

The Usadel equation is complemented by the boundary cond itions presented bv lKuprianov and Lukichevl l)l988|) on 
the basis of the Zaitsev's boundary conditions llZaitsevl I l984l). Various aspects of the bou ndary conditions have been 
discussed bv lKopu et all l)2004j) : Eambert et all l)l997|) : iNazarovl l)l999j) : IXia et all l)2002() . In the absence of spin-flip 
processes at the interface they take the form: 

gAgi = — [51,52] , (A21) 

27a 

where 71 = -R^cti, o\ is the conductivity of the conductor 1 and Ri, is the interface resistance per unit area, the 
x-coordinate is assumed to be normal to the plane of the interface. 

The boundary condition (| A2 1|> implies that we accept the simplest model of the S/F interface which is used in 
most papers on S/F structures. We assume that the interface separates two dirty regions: a singlet superconductor 
and a ferromagnet. The superconductor and the ferromagnet are described in the mean field approximation with 
different order parameters: the off-diagonal order parameter A in the superconductor (in the weak coupling limit) 
and the exchange field h in the ferromagnet acting on the spins of free electrons. No spin-flip scattering processes 
are assumed at the S/F interface. A gener alization of the boundary conditio n s we use to the case of a spin-act i ve S/F 
interface was c arried out i n the papers bv lEschrid l)2000j) : iFogelstroml l|2000j) : IKopu et all (|2004|k iMillis et all i|l98Sft . 

Eqs. I|A18[) and l|A21|l together with the self-consistency equation that determines the superconducting order 
parameter A, constitute a complete set of equations from which one can obtain the Green's functions. 

The Usadel equation can be solved in some particular cases. We often use the linearized Usadel equation. In order 
to obtain the linearized Usadel equation we represent the Green's functions g in the superconductor in the form 

3(p F ,r;w) = gscsiu) + Sgs + Sf s , (A22) 

where ()bcs(u) = T 3 g B cs(^) + ihfBCS, 9bcs(^) = (!w/A)/ BCS , f B cs = A/Wuj 2 + A 2 . We have written the 
matrix g in the so called Matsubara representation. This means that a substitution e iu; (u> = ttT(2ji + l),n — 

0, ±1,±2, ) is done and g(uj) coincides with g R (e) for positive u> and with g (e) for negative u>. The linearized 

Usadel equation has the form 



in the S region and 



dlJfs ~ K 2 8 5fs = 2i(SA/D s )g% cs (A23) 



dlJf - K 2 Jf + in 2 h {[d- 3 ,Sf] + cosa±T 3 [d- 2 ,Sf]-sma} =0 (A24) 



in the F region. Her e re| = 2E U /D S> nl = 2\uj\/D f , n\ = h sgnu/D F and [A, B]± = AB ± BA, SA = ir 2 a 3 5A. The 
signs ± in Eq. lA24p correspond to the right and left layer respectively. 

The boundary conditions for 5fs and Sf F = f (in zero-order approximation 5f F = 0) are obtained from Ea. (|A2ip . 
They have the form 

d x 5fs = {^hs)[9Bcs 5 f ~ 9BcsfBCS^39F3 ~ 9bcsIs] (A25) 
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dJ F = (lh F )[9BcsSf - fs] (A26) 

where jf,s = RbCF,s- 

If the Green's functions are known, one can calculate macroscopic quan tities such as the current, magnetic moment 
etc. For example, the current is given bv lLarkin and Ovchinnikovl l)l984j) 

I 8 = (L y L z /l6)a F Tr(T 3 a ) J de{g s dg s /dx) 12 (A27) 

where L y ^ z are the widths of the films in y and z direction (the current flows in the transverse ir-direction) and 
subscript (12) shows that one has to take the Keldysh component of the supermatrix g s dg s /dx. A variation of the 
magnetic moment due to proximity effect is determined by formulae 

SM Z = yL B v{l/2)iirT^Tr{T 3 (g> a 3 5g) (A28) 



5M x , y = ^(I^mT^) IY(t (8 a ia 8g) (A29) 

where v is the density-of-states at the Fermi level in the normal state and fis — gfJ-Bohr is an effective Bohr magneton. 

Finally, it is important to make remarks concerning the notations used in this review. In most works where the 
S/F structures with homogeneous magnetization are studied, the Green's function g is a 2 x 2 matrix with the usual 
normal and Gor'kov's components. Of course, this simplification can be made provided the magnetizations of the F 
layers involved in the problem are aligned in one direction. However, this simple form leads to erroneous results 
if the magnetizations are arbitrarily oriented with respect to each other. The 4x4 form of the Green's function is 
unavoidable if one studies structures with a no n-homogeneo us magnetization. Of course, the c- operators in Eq. (|A3|) 
can be defined in different ways. For example. iMaki (J.1969) introduced a spinor representation of the field operators, 
which is equivalent to letting in Eq. IIA3I) the spin index of the operator a unchanged when n = 2. This notation was 
used in later works ({Alexander et q/|ll985tlDemler et q/.Ul99?l e.g) in which the Green's functions have a 2 x 2 block 
matrix form. The diagonal blocks represent the normal Green's functions, while the off-diagonal blocks represent the 
anomalous one. With this notation the matrix, Eq. (|A5|) . changes its form. For example, the term containing A is 
proportional to iai and not to 173. The choice of the notation depends on the problem to solve. In order to study 
the triplet superconductivity induced in S/F systems and to see explicitly the three projections {S z = 0, ±1) of the 
conde nsate function it is more co nvenient to use the operators defined in Eci . (|A3(I (see for example iFominov et al\ 
<|2003ft and lBergeret et al\ l)2001cj) V 



APPENDIX B: Future direction of the experimental research 

As we have seen throughout the paper there is a great number of experiments on S/F structures. The variety of 
superconducting and ferromagnetic materials is very large. In this section we review briefly some of these experiments. 
We will not dwell on specific fabrication techniques but rather focus on the discussion: which pairs of material (S and 
F) are more appropriate for the observation of the effects studied in this review. 

First experiments on S/F structures used stron g ferromagnets (lar ge exchange fields) as Fe, Ni, Co or Gd and 
conventional superconductors like Nb, Pb, V, etc. ijHauser et aZi ll9631. In these experiments the dependence of the 
superconducting transition temperature on the thicknesses of the S and F layers has been measured. In other words 
the suppression of the superconductivity due to the strong exchange field of the ferromagnet was analyzed. It is clear 
that for such strong ferromagnets the spin splitting is large and therefore a mismatch in electronic parameters of the 
S and F region s is lar ge. This leads to a low interface transparency and a weak proximity effect. This was confirmed 
bv lAarts et all l|l997l) in experiments on V/Vi- x Fe x multilayers. Varying the concentration of Fe in the VFe alloys 
they could change the values of the exchange field and indirectly the transparency of the interface. Such systems 
consisting of a conventional superconductor and a ferromagnetic alloy, both with similar band structure (in the above 
experiment the mismatch was < 5%) , are good candidates for observing the effects discussed in sections llV. Al . IV~Bl 
andlvTTl 

Weak ferromagnets have been used in the last years in many experiments on S/F structur es. Before we turn ou r 
attention to ferromagnets with small exchange fields it is worth mentioning the experiment bv lRusanov et all l)2004f) . 
They analyzed the so-called spin switch effect. In particular they studied the transport properties of Permalloy (Py) /Nb 
bilayers. They observed an enhancement of superconductivity in the resistive transition in the field range where the 
magnetization of the Py switches and many domains are present. Interesting for us is that Py shows a well defined 
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magnetization switching at low fields and therefore it could be used in order to detect the long-range triplet component 
that appears when the magnetization of the ferromagnet is not homogeneous (see section llll.Cfl . Finally, a magnetic 
configuration analy si s of the strong-ferroma gnetic structures used in transport experiments as those performed by 
iGiroud et all lll998ft : iPetrashov et all §999) may also serve to confirm the predictions of section UlI.CI As it was 
discussed before the increase in the conductance of the ferromagnet for temperatures below the superconducting T c 
may be explained assuming a long-range proximity effect. 

The proximity effect in S/F is stronger if one uses dilute ferromagnetic alloys. Thus, such materials are the 
best candidates in order to observed most of t he effects discussed in this review. The idea of using ferromagnetic 
alloys with small exchanges field was used bv iRvazanov et all <)200l[) . They were the first in observing the sign 
reversal of the critical current in a S/F/S Josephson junction. Nb was used as supercond uctor while C u n.4RN in.52 
alloy as a ferromagnet (exchange field ~ 25K). (Later on similar results were obtaine d bvlKon tos et all l)2002|) on 
Nb/Al/Al 2 3 /PdNi/Nb structures). The CuNi alloy was also used in the experiment bv lGu et all l)2002bj) on F/S/F 
structures. In this experiment the authors determined the dependence of the superconducting transition temperature 
on the relative magnetization-orientation of the two F layers. In order to get different alignments between the two 
CuNi layers an exchange-biased spin-valve stack of CuNi/Nb/CuNi/FesnMnso was employed. With a small magnetic 
field the authors could switch the magnetization direction of the free NiCu layer. This technique could be very useful 
in order to observe the Josephson coupling via the triplet component a s described in section llV.AI 

Finally, it is worth mentioning the experiment by Stahn et all l)2005|) on YBa 2 Cu307/La2 /3 Ca 1/3 Mn03. Using the 
neutron reflectometry technique they observed a induced magnetic moment in the superconductor. Although the 
materials employed in thi s experiment canno t be quantitatively described with the methods presented in this review 
(the ferromagnet used bv IStahn et all 1120051) is a half-metal with a exchange field comparable to the Fermi-energy 
and the superconductor is unconventional), the experimental technique may be used in other experiments in order to 
detected the induced magnetization predicted in sections IV.BI and IV. Cl 



APPENDIX: List of Symbols and Abbreviations 



s 


superconductor 


N 


nonmagnetic normal metal 


F 


ferromagnetic metal 


I 


insulator 


LRTC 


long-range triplet component 


ft, I = 1,2,3 


pauli matrices in particle-hole space 


&i, i = 1,2,3 


Pauli matrices in spin space 


to, <5o 


unit matrices. 


D 


diffusion coefficient 


V 


density of states 


w = nT(2n+ 1) 


Matsubara Frequency 


e 


real frequency (energy) 


9bcs 


quasiclassical normal Green's function for a bulk superconductor 


Jbcs 


quasiclassical anomalous Green's function for a bulk superconductor 


T c 


superconducting critical temperature 


Ic 


Josephson critical current 


R b 


interface resistance per unit area 


f-bN = DN/2Rb<TNdN 


minigap induced in a normal metal 


TS,F 


conductivity in the normal state 


1S,F 


Rb&S.F 


7 


ratio opjos 


J 


magnetic coupling between localized magnetic moments. 


h 


exchange field acting on the spin of conducting electrons 




characteristic penetration length of the condensate into a dirty normal 
metal 




characteristic penetration length of the condensate into a dirty ferro- 
magnet 

superconducting coherence length for a dirty superconductor 
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